Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,587,752 T→G 5.2% V182G (GTT→GGT)  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,7520TG 5.2% 436.1 / 11.2 136V182G (GTT→GGT) catA‑IIcatechol 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (67/61);  new base G (4/3);  total (72/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 145 total aligned reads displayed.

CACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAG  >  NC_002947/3587603‑3587901
                                                                                                                                                     |                                                                                                                                                     
 aCGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGtt                                                                                                                                                      >  1:2093235/1‑150 (MQ=255)
 aCGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGtt                                                                                                                                                      >  2:266993/1‑150 (MQ=255)
   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)
                      cAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGGTAACGATGcgcc                                                                                                                                           >  2:308188/1‑137 (MQ=255)
                         gTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACcgc                                                                                                                              <  2:2093235/150‑1 (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                                                                                                         <  2:1275869/150‑1 (MQ=255)
                                              ccGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCg                                                                                                         >  1:1099413/1‑150 (MQ=255)
                                              ccGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCg                                                                                                         >  1:1669536/1‑150 (MQ=255)
                                                gCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCggg                                                                                                       <  1:1159368/150‑1 (MQ=255)
                                                   gCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTa                                                                                                    <  1:836560/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                                                                                           >  1:685494/1‑150 (MQ=255)
                                                            cACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGa                                                                                           >  2:1183794/1‑150 (MQ=255)
                                                            cACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGa                                                                                           >  1:1796696/1‑150 (MQ=255)
                                                             aCGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGc                                                                                                                                              >  2:2104952/1‑98 (MQ=255)
                                                              aggTCGACCGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTAt                                                                                                   >  2:1688695/2‑139 (MQ=255)
                                                              aggTCGACCGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTAt                                                                                                   <  1:1688695/138‑1 (MQ=255)
                                                                       tgtgGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGc                                                                                <  2:412870/150‑1 (MQ=255)
                                                                           gCATGGGAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAAc                                                                            >  2:2213193/1‑150 (MQ=255)
                                                                           gCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAgg                                                                                                                                         >  1:2334375/1‑89 (MQ=255)
                                                                           gCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAAc                                                                            >  1:1401680/1‑150 (MQ=255)
                                                                            cATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAAcc                                                                           >  2:875936/1‑150 (MQ=255)
                                                                               gCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGa                                                                                                                                                 <  1:766888/77‑1 (MQ=255)
                                                                               gCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGa                                                                                                                                                 >  2:766888/1‑77 (MQ=255)
                                                                               gCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAg                                                                        <  2:685494/150‑1 (MQ=255)
                                                                                 cAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCaa                                                                             <  1:1297630/143‑1 (MQ=255)
                                                                                 cAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCaa                                                                             >  2:1297630/1‑143 (MQ=255)
                                                                                     aCCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGTGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGc                                                                  <  2:293520/150‑1 (MQ=255)
                                                                                     aCCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGc                                                                  <  2:1865510/150‑1 (MQ=255)
                                                                                      cccGTGGTACCTACTCGTTCTTCGACCAAAGCCAGCCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGcc                                                                 <  1:1334197/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                                                                >  1:280086/1‑150 (MQ=255)
                                                                                       ccGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                <  2:786111/150‑1 (MQ=255)
                                                                                           ggTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGAc                                                            <  1:1536363/150‑1 (MQ=255)
                                                                                              aCCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                >  2:1743214/1‑143 (MQ=255)
                                                                                              aCCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACctg                                                         <  2:2387360/150‑1 (MQ=255)
                                                                                                cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGt                                                                                                                <  1:206123/93‑1 (MQ=255)
                                                                                                cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGt                                                                                                                >  2:206123/1‑93 (MQ=255)
                                                                                                cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTAt                                                                                                   <  1:526801/106‑1 (MQ=255)
                                                                                                cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                >  2:966933/1‑141 (MQ=255)
                                                                                                cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCt                                                                <  1:966933/141‑1 (MQ=255)
                                                                                                cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACctgct                                                       >  2:1504918/1‑150 (MQ=255)
                                                                                                   cTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTAc                                                                                                                                                    >  1:659042/1‑54 (MQ=255)
                                                                                                    tCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTAt                                                                                                   >  2:485779/1‑102 (MQ=255)
                                                                                                    tCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTAt                                                                                                   <  1:485779/102‑1 (MQ=255)
                                                                                                    tCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGc                                                   >  1:850690/1‑150 (MQ=255)
                                                                                                      gTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGc                                                                                                      >  2:1672264/1‑97 (MQ=255)
                                                                                                      gTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCg                                                 >  2:1066360/1‑150 (MQ=255)
                                                                                                        tcttcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTc                                               <  1:68944/150‑1 (MQ=255)
                                                                                                          ttcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCg                                                                                                         <  1:2243130/90‑1 (MQ=255)
                                                                                                          ttcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCAt                                             <  2:2346569/150‑1 (MQ=255)
                                                                                                          ttcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCAt                                             >  2:1180194/1‑150 (MQ=255)
                                                                                                          ttcGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCAt                                             >  1:2346569/1‑150 (MQ=255)
                                                                                                             gACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGc                                          >  2:954954/1‑150 (MQ=255)
                                                                                                             gACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGc                                          >  1:1107757/1‑150 (MQ=255)
                                                                                                                 aaaGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCt                                                                                               >  2:1397794/1‑93 (MQ=255)
                                                                                                                 aaaGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAgc                                      <  1:1248988/150‑1 (MQ=255)
                                                                                                                 aaaGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAg                                       >  2:305329/1‑149 (MQ=255)
                                                                                                                 aaaGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAg                                       <  1:305329/149‑1 (MQ=255)
                                                                                                                    gCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGcc                                   >  1:1239527/1‑150 (MQ=255)
                                                                                                                     ccAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGccc                                  >  1:1711951/1‑150 (MQ=255)
                                                                                                                        gTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTcc                                                                                                                     <  2:180804/64‑1 (MQ=255)
                                                                                                                        gTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTcc                                                                                                                     >  1:180804/1‑64 (MQ=255)
                                                                                                                        gTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGgc                               >  2:533/1‑150 (MQ=255)
                                                                                                                         tCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGgcg                              >  2:348303/1‑150 (MQ=255)
                                                                                                                               tACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTg                        >  1:1283877/1‑150 (MQ=255)
                                                                                                                                     cTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGccc                                  <  2:681937/134‑1 (MQ=255)
                                                                                                                                     cTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTg                        <  1:335722/144‑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                                                     <  1:417239/104‑1 (MQ=255)
                                                                                                                                                    gTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCg         <  1:1150167/144‑1 (MQ=255)
                                                                                                                                                    gTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCg         >  2:1150167/1‑144 (MQ=255)
                                                                                                                                                     ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAAt                                                                    <  1:1842664/84‑1 (MQ=255)
                                                                                                                                                     ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAAt                                                                    <  1:643818/84‑1 (MQ=255)
                                                                                                                                                     ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAAt                                                                    >  2:643818/1‑84 (MQ=255)
                                                                                                                                                     ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAAt                                                                    >  2:1842664/1‑84 (MQ=255)
                                                                                                                                                     ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAg  <  2:1107757/150‑1 (MQ=255)
                                                                                                                                                     ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAg  <  2:1107864/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAG  >  NC_002947/3587603‑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: