Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,587,755 C→A 5.4% T183N (ACC→AAC)  catA‑II → catechol 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,587,7550CA 5.4% 458.8 / 8.7 139T183N (ACC→AAC) catA‑IIcatechol 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (68/63);  new base A (5/3);  total (73/66)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

GGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAG  >  NC_002947/3587606‑3587901
                                                                                                                                                     |                                                                                                                                                  
ggTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTAc                                                                                                                                                    <  1:1836774/150‑1 (MQ=255)
  tGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCg                                                                                                                                                  >  1:2387360/1‑150 (MQ=255)
    aTGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGAt                                                                                                                                                <  2:450143/150‑1 (MQ=255)
      gTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGc                                                                                                                                              <  2:2364523/150‑1 (MQ=255)
        tCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGGTAACGATGcgc                                                                                                                                            <  1:1097645/150‑3 (MQ=255)
          cTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGACGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAg                                                                                                                                          >  2:1411559/1‑150 (MQ=255)
           tGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAgg                                                                                                                                         >  1:1865510/1‑150 (MQ=255)
           tGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCCTTACCGATGCCCAgg                                                                                                                                         >  2:554578/1‑150 (MQ=255)
             gAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAgggg                                                                                                                                       >  2:2206045/1‑150 (MQ=255)
              aaGGCCAGGTGCTGGACCCGCCCGGACGCGCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGgc                                                                                                                                      <  1:731687/150‑1 (MQ=255)
              aaGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGgc                                                                                                                                      >  2:68944/1‑150 (MQ=255)
                   cAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGGTAACGATGcgcc                                                                                                                                           <  1:308188/140‑4 (MQ=255)
                   cAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGGTAACGATGcgcc                                                                                                                                           >  2:308188/1‑137 (MQ=255)
                      gTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACcgc                                                                                                                              <  2:2093235/150‑1 (MQ=255)
                          tGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACcgcgcgc                                                                                                                          >  1:1808209/1‑150 (MQ=255)
                               ccGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCATACACACTCAATCGGACCATCGCTAATGACGCCCACGGGCGCTACCTCGCGCGCTcc                                                                                                                     >  2:872900/1‑150 (MQ=255)
                                 gCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGcccggt                                                                                                                                         <  1:2269100/128‑4 (MQ=255)
                                 gCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAgg                                                                                                                                         >  2:2269100/1‑128 (MQ=255)
                                     ggACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAACACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTg                                                                                                               <  2:1740163/150‑1 (MQ=255)
                                          cccGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATc                                                                                                          <  2:1808209/150‑1 (MQ=255)
                                           ccGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCg                                                                                                         >  1:1669536/1‑150 (MQ=255)
                                           ccGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCg                                                                                                         <  2:1275869/150‑1 (MQ=255)
                                           ccGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCg                                                                                                         >  1:1099413/1‑150 (MQ=255)
                                            cGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGa                                                                                                                                                 <  2:2386650/109‑1 (MQ=255)
                                             gCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCggg                                                                                                       <  1:1159368/150‑1 (MQ=255)
                                                 ccGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTAt                                                                                                   >  2:1671623/1‑150 (MQ=255)
                                                     gTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCt                                                                                               <  1:1392643/150‑1 (MQ=255)
                                                        ccACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCCACGCCAGTCGGCGTCCACCATGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCg                                                                                            >  1:1498068/1‑150 (MQ=255)
                                                         cACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTGGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGa                                                                                           >  1:561988/1‑150 (MQ=255)
                                                         cACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGa                                                                                           >  2:1334197/1‑150 (MQ=255)
                                                         cACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGa                                                                                           >  1:1796696/1‑150 (MQ=255)
                                                          aCGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGc                                                                                                                                              <  1:2104952/98‑1 (MQ=255)
                                                          aCGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGc                                                                                                                                              >  2:2104952/1‑98 (MQ=255)
                                                                   cTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAgggg                                                                                 <  2:1886373/150‑1 (MQ=255)
                                                                        gCATGGGAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAAc                                                                            >  2:2213193/1‑150 (MQ=255)
                                                                        gCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAgg                                                                                                                                         >  1:2334375/1‑89 (MQ=255)
                                                                        gCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAgg                                                                                                                                         <  2:2334375/89‑1 (MQ=255)
                                                                        gCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAAc                                                                            >  1:1401680/1‑150 (MQ=255)
                                                                         cATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAAcc                                                                           >  2:875936/1‑150 (MQ=255)
                                                                            gCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGa                                                                                                                                                 >  2:766888/1‑77 (MQ=255)
                                                                            gCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGa                                                                                                                                                 <  1:766888/77‑1 (MQ=255)
                                                                              cAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCaa                                                                             <  1:1297630/143‑1 (MQ=255)
                                                                                  aCCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGTGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGc                                                                  <  2:293520/150‑1 (MQ=255)
                                                                                  aCCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGc                                                                  <  2:1865510/150‑1 (MQ=255)
                                                                                  aCCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGc                                                                  <  2:1099413/150‑1 (MQ=255)
                                                                                    ccGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGc                                                                                >  1:555769/1‑134 (MQ=255)
                                                                                    ccGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                >  1:786111/1‑150 (MQ=255)
                                                                                    ccGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                >  2:582031/1‑150 (MQ=255)
                                                                                    ccGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                <  2:786111/150‑1 (MQ=255)
                                                                                    ccGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                >  1:280086/1‑150 (MQ=255)
                                                                                          tACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACct                                                          >  1:1887552/1‑150 (MQ=255)
                                                                                           aCCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                <  1:1743214/143‑1 (MQ=255)
                                                                                           aCCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                >  2:1743214/1‑143 (MQ=255)
                                                                                             cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTAt                                                                                                   >  2:526801/1‑106 (MQ=255)
                                                                                             cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTAt                                                                                                   <  1:526801/106‑1 (MQ=255)
                                                                                             cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                >  2:966933/1‑141 (MQ=255)
                                                                                             cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACctgct                                                       >  2:1504918/1‑150 (MQ=255)
                                                                                             cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACctgct                                                       <  2:280086/150‑1 (MQ=255)
                                                                                              tACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACctgctg                                                      >  2:1248988/1‑150 (MQ=255)
                                                                                                cTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTAc                                                                                                                                                    >  1:659042/1‑54 (MQ=255)
                                                                                                cTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTAc                                                                                                                                                    <  2:659042/54‑1 (MQ=255)
                                                                                                 tCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTAt                                                                                                   <  1:485779/102‑1 (MQ=255)
                                                                                                 tCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGc                                                   >  1:850690/1‑150 (MQ=255)
                                                                                                   gTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGc                                                                                                      >  2:1672264/1‑97 (MQ=255)
                                                                                                   gTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGc                                                                                                      <  1:1672264/97‑1 (MQ=255)
                                                                                                   gTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCg                                                 >  2:1066360/1‑150 (MQ=255)
                                                                                                   gTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCg                                                 >  1:1319410/1‑150 (MQ=255)
                                                                                                   gTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCg                                                 >  1:1277187/1‑150 (MQ=255)
                                                                                                     tcttcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTc                                               <  1:68944/150‑1 (MQ=255)
                                                                                                      cttcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTggccatgac                                              >  2:1698800/1‑143 (MQ=255)
                                                                                                       ttcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCg                                                                                                         <  1:2243130/90‑1 (MQ=255)
                                                                                                       ttcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCAt                                             >  1:2346569/1‑150 (MQ=255)
                                                                                                       ttcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCAt                                             >  2:1180194/1‑150 (MQ=255)
                                                                                                        tcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATg                                            <  1:582031/150‑1 (MQ=255)
                                                                                                          gACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGc                                          >  1:1107757/1‑150 (MQ=255)
                                                                                                          gACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGc                                          >  2:954954/1‑150 (MQ=255)
                                                                                                            ccAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCa                                        <  1:875936/150‑1 (MQ=255)
                                                                                                              aaaGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCt                                                                                               >  2:1397794/1‑93 (MQ=255)
                                                                                                              aaaGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAggg                                                                                  >  1:799757/1‑106 (MQ=255)
                                                                                                              aaaGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAgc                                      <  1:1248988/150‑1 (MQ=255)
                                                                                                                 gCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGcc                                   >  1:1239527/1‑150 (MQ=255)
                                                                                                                  ccAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGccc                                  >  1:1711951/1‑150 (MQ=255)
                                                                                                                     gTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTcc                                                                                                                     <  2:180804/64‑1 (MQ=255)
                                                                                                                     gTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGgc                               >  2:533/1‑150 (MQ=255)
                                                                                                                      tCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGgcg                              >  2:348303/1‑150 (MQ=255)
                                                                                                                                  cTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGccc                                  >  1:681937/1‑134 (MQ=255)
                                                                                                                                  cTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGccc                                  <  2:681937/134‑1 (MQ=255)
                                                                                                                                  cTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTg                        <  1:335722/144‑1 (MQ=255)
                                                                                                                                  cTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTg                        >  2:335722/1‑144 (MQ=255)
                                                                                                                                   tGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCActtct                 <  1:1504918/150‑1 (MQ=255)
                                                                                                                                      ccctgggcATCGGTAACGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTa              >  2:1486418/7‑150 (MQ=255)
                                                                                                                                        cGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTAtc            <  2:561988/150‑1 (MQ=255)
                                                                                                                                        cGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTAtc            <  2:2256881/150‑1 (MQ=255)
                                                                                                                                          gggcATCGGTAACGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTAtctc          >  1:1356847/3‑150 (MQ=255)
                                                                                                                                             cATCGGTAACGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTgg                                                     >  2:417239/1‑104 (MQ=255)
                                                                                                                                             cATCGGTAACGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTgg                                                     <  1:417239/104‑1 (MQ=255)
                                                                                                                                                 gTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCg         >  2:1150167/1‑144 (MQ=255)
                                                                                                                                                 gTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCg         <  1:1150167/144‑1 (MQ=255)
                                                                                                                                                  ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAAt                                                                    >  2:1842664/1‑84 (MQ=255)
                                                                                                                                                  ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAg  <  2:1107864/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAG  >  NC_002947/3587606‑3587901

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