Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,182,136 C→T 5.2% R91H (CGT→CAT)  ygiD ← 4,5‑ DOPA‑extradiol‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,182,1360CT 5.2% 313.6 / 4.1 96R91H (CGT→CAT) ygiD4,5‑ DOPA‑extradiol‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (43/48);  new base T (2/3);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

CAGGTTTGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCAC  >  NC_000913/3181990‑3182281
                                                                                                                                                  |                                                                                                                                                 
cAGGTTTGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGc                                                                                                                                                 <  2:169921/149‑1 (MQ=255)
 aGGTTTGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCt                                                                                                                                                <  1:142857/149‑1 (MQ=255)
 aGGTTTGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCt                                                                                                                                                >  8:341783/1‑149 (MQ=255)
  ggTTTGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCtg                                                                                                                                               <  7:342344/149‑1 (MQ=255)
   gTTTGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCtgt                                                                                                                                              >  2:191728/1‑149 (MQ=255)
        atACTGGCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGAGCGGCGCTAACAGCGCAACCAGACGCTGTGCCa                                                                                                                                          <  3:79255/147‑1 (MQ=255)
         tACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAgc                                                                                                                                        >  3:219205/1‑149 (MQ=255)
          aCTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAgcg                                                                                                                                       >  5:245747/1‑149 (MQ=255)
          aCTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAgcg                                                                                                                                       >  1:92419/1‑149 (MQ=255)
           cTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAgcgc                                                                                                                                      <  4:86714/149‑1 (MQ=255)
              tCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAgg                                                                                                                                   <  2:101738/149‑1 (MQ=255)
              tCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAg                                                                                                                                    <  4:123193/148‑1 (MQ=255)
               cGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGc                                                                                                                                  <  1:62912/149‑1 (MQ=255)
                gATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGc                                                                                                                                  >  6:219650/1‑148 (MQ=255)
                    cTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCACCCAGACGCTGTGCCAGCGCAGGCGAAcc                                                                                                                             >  6:277429/1‑149 (MQ=255)
                     tCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCAACGATCCGTGGTCAAAGACCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGATGTGCCAGCGCAGGCGAAacc                                                                                                                            >  2:6971/1‑149 (MQ=255)
                       aaCTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCgg                                                                                                                          <  3:274105/149‑1 (MQ=255)
                       aaCTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCg                                                                                                                           >  7:77421/1‑148 (MQ=255)
                            cACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGc                                                                                                                       <  5:8453/147‑1 (MQ=255)
                            cACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAg                                                                                                                     <  5:219650/149‑1 (MQ=255)
                                aTCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGAt                                                                                                                  >  1:105256/1‑148 (MQ=255)
                                aTCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGAt                                                                                                                  >  5:305284/1‑148 (MQ=255)
                                aTCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGAt                                                                                                                  >  3:104866/1‑148 (MQ=255)
                                    ggATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTTGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAAGCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGc                                                                                                             <  1:6971/149‑1 (MQ=255)
                                    ggATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGc                                                                                                             >  5:81459/1‑149 (MQ=255)
                                    ggATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGc                                                                                                             >  6:71501/1‑149 (MQ=255)
                                    ggATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGc                                                                                                             <  4:118915/149‑1 (MQ=255)
                                      atatCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGATGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGt                                                                                                           <  1:304850/149‑1 (MQ=255)
                                       tatCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTa                                                                                                          <  6:245747/149‑1 (MQ=255)
                                       tatCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTa                                                                                                          <  7:345938/149‑1 (MQ=255)
                                        atCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTAt                                                                                                         >  6:80795/1‑149 (MQ=255)
                                         tCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCa                                                                                                                                          <  8:39323/115‑1 (MQ=255)
                                         tCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATc                                                                                                        >  8:86990/1‑149 (MQ=255)
                                              gTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGTCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACa                                                                                                   <  6:149490/149‑1 (MQ=255)
                                              gTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTAt                                                                                                         <  5:163649/143‑1 (MQ=255)
                                              gTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTAt                                                                                                         >  6:163649/1‑143 (MQ=255)
                                                      aCATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGa                                                                                                                         >  3:276546/1‑119 (MQ=255)
                                                      aCATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGa                                                                                                                         <  4:276546/119‑1 (MQ=255)
                                                      aCATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGc                                                                                           >  1:207536/1‑149 (MQ=255)
                                                         tCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCgg                                                                                         >  6:114578/1‑148 (MQ=255)
                                                             aaTCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACTCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGaa                                                                                      <  5:103868/147‑1 (MQ=255)
                                                             aaTCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGaa                                                                                      <  5:176979/147‑1 (MQ=255)
                                                             aaTCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAg                                                                                     <  5:80795/148‑1 (MQ=255)
                                                              aTCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGCCGCTGTGCCAGCGCAGGCGAACCCGGCGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGcc                                                                                   >  8:74348/1‑149 (MQ=255)
                                                               tCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGATGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGa                                                                                       <  4:305186/144‑1 (MQ=255)
                                                               tCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGATGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGa                                                                                       >  3:305186/1‑144 (MQ=255)
                                                               tCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGcca                                                                                  >  8:210602/1‑149 (MQ=255)
                                                               tCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGcca                                                                                  <  6:305284/149‑1 (MQ=255)
                                                                 aGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGccac                                                                                 >  3:183357/1‑148 (MQ=255)
                                                                  gCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGccacc                                                                                >  8:282504/1‑148 (MQ=255)
                                                                  gCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGccacc                                                                                >  8:130786/1‑148 (MQ=255)
                                                                  gCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGccacc                                                                                >  2:340065/1‑148 (MQ=255)
                                                                  gCACCCCCCAAGAGCCGTGGTCAAAGCCACATGCTTCTTTATCCAGCGTCAAAGGGATAGTCGCTAACAGCTAAACCAGACGCTGTGCCAGAGCAGGCGAACCATCAGAAGGATAATGAGTAACGTAAAACGCCTTCGGTAAGacaac                                                                                >  8:258293/1‑146 (MQ=255)
                                                                    aCGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCaaa                                                                             <  2:99605/149‑1 (MQ=255)
                                                                       ccccccGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCCGCGGCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTc                                                                          <  4:104866/149‑1 (MQ=255)
                                                                        cccACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCa                                                                         <  7:282504/149‑1 (MQ=255)
                                                                        cccACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCa                                                                         <  2:75803/149‑1 (MQ=255)
                                                                         ccACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAgcgc                                                                                               >  1:242709/1‑126 (MQ=255)
                                                                         ccACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAgcgc                                                                                               <  2:242709/126‑1 (MQ=255)
                                                                         ccACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCAt                                                                        <  1:191728/149‑1 (MQ=255)
                                                                          cACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCCCCGGGATCGGCGCTAACAACTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATg                                                                       <  1:277756/149‑1 (MQ=255)
                                                                          cACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCAt                                                                        >  1:258208/1‑148 (MQ=255)
                                                                              aGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATc                                                                   <  1:340065/149‑1 (MQ=255)
                                                                                 cggggTCAAAGCCCCCCCCTGCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCACCCCGCCCCCGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCAGCGGGAAGCCACCACCGTCATGACTCGTg                                                                <  2:92419/146‑1 (MQ=255)
                                                                                   tGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTggg                                                              <  4:183357/149‑1 (MQ=255)
                                                                                            ccccAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGtctc                                                     <  7:341783/149‑1 (MQ=255)
                                                                                            ccccAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGtctc                                                     >  6:141434/1‑149 (MQ=255)
                                                                                             cccAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCATCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGtctc                                                     >  4:9740/1‑148 (MQ=255)
                                                                                               cAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGCCGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAACGTCATGAATCGGGGGCGGCGTCTccct                                                  >  6:121591/1‑147 (MQ=255)
                                                                                                   cttcttTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATcgcg                                              <  6:81459/149‑1 (MQ=255)
                                                                                                   cttcttTATCCAGCGGCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATcgcg                                              <  7:210602/149‑1 (MQ=255)
                                                                                                      cttTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGCCGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGt                                            <  5:8212/148‑1 (MQ=255)
                                                                                                      cttTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCgg                                             <  2:258208/147‑1 (MQ=255)
                                                                                                      cttTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGt                                            <  2:48179/148‑1 (MQ=255)
                                                                                                       tttATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAAccc                                                                                                                            >  3:256012/1‑67 (MQ=255)
                                                                                                       tttATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAAccc                                                                                                                            <  4:256012/67‑1 (MQ=255)
                                                                                                           tCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTc                                       <  5:277429/148‑1 (MQ=255)
                                                                                                             cAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTg                                    >  6:313292/1‑149 (MQ=255)
                                                                                                              aGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGt                                   <  7:86990/149‑1 (MQ=255)
                                                                                                               gCGTCACCGGGATCGGCGCTAACAGCTCAACCAGATGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGtt                                  >  8:26282/1‑149 (MQ=255)
                                                                                                               gCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGCACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTggt                                  >  1:102934/1‑147 (MQ=255)
                                                                                                                     ccGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACgg                            <  8:77421/149‑1 (MQ=255)
                                                                                                                      cGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACggg                           >  3:147472/1‑149 (MQ=255)
                                                                                                                      cGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACggg                           <  4:219205/149‑1 (MQ=255)
                                                                                                                           tCGGCGCTAACAGCTCAACCAGATGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAAc                      <  7:26282/149‑1 (MQ=255)
                                                                                                                              gcgcTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGCACCCGGAGCAGGATAATGCGTATCGTACCGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAg                   >  7:120536/1‑149 (MQ=255)
                                                                                                                               cgcTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGt                  <  2:105256/149‑1 (MQ=255)
                                                                                                                                   aaCAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCa                                                                         <  7:8705/90‑1 (MQ=255)
                                                                                                                                   aaCAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCa                                                                         >  8:8705/1‑90 (MQ=255)
                                                                                                                                   aaCAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGc              <  4:140031/149‑1 (MQ=255)
                                                                                                                                   aaCAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGc              >  4:56121/1‑149 (MQ=255)
                                                                                                                                         tCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAac        >  3:209090/1‑149 (MQ=255)
                                                                                                                                           aaCCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCGTGTTCCACGGGTAAACCAGTGAGCCGAAac        <  6:143635/147‑1 (MQ=255)
                                                                                                                                             ccAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAaccac     >  3:14785/1‑148 (MQ=255)
                                                                                                                                             ccAGACGCTCTCCCAGCGCACGCGAACCCGGAGCAGGATAATTCGTATCGTACAGCGCATGCGGGAAGCCACCAAAGTCATGAATCTTGGGCGGCGTCTCCATCGCGGTCACTCCAGTACCACGGCTAAACCACTGAGCCGAAacaacc    >  6:175624/1‑149 (MQ=255)
                                                                                                                                                gACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTACACGGGTAAACCAGTGAGCCGACaccacccc  >  3:184931/1‑146 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CAGGTTTGCTACTGTCGATACTCAACTGCACCATCGGGATATCAGCGTCAGGATACATCTTAATCAGCACGCCCCACGAGCCGTGGTCAAAGCCCCAGGCTTCTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCAC  >  NC_000913/3181990‑3182281

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑