Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,489 C→A 5.7% V142V (GTG→GTT amaB ← L‑piperidine‑6‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,4890CA 5.7% 585.8 / 16.7 175V142V (GTG→GTTamaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (76/89);  new base A (5/5);  total (81/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

CTTTTTCGAACAACGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTT  >  NC_002947/6005341‑6005638
                                                                                                                                                    |                                                                                                                                                     
cTTTTTCGAACAACGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACca                                                                                                                                                      <  2:2998541/150‑1 (MQ=255)
 tttttCGAACAACGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACcac                                                                                                                                                     >  1:5073871/1‑150 (MQ=255)
  ttttCGAACAACGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACcaca                                                                                                                                                    >  1:125621/1‑150 (MQ=255)
    ttCGAACAACGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACAcc                                                                                                                                                  <  1:3637225/150‑1 (MQ=255)
    ttCGAACAACGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACAcc                                                                                                                                                  >  1:1430125/1‑150 (MQ=255)
    ttCGAACAACGCCTGGCAGGCCAGGGCGGTCAGCGAACACAACACCCACCCAATCCAAACCACCGACTAGCCCGCCACCACGGCCAGGCCGGTGTTCCAGGCCCACACGCCAACCGGGAAGTTCAAGGCGCTGATCACGCCGACCACAcc                                                                                                                                                  >  2:5256829/1‑150 (MQ=255)
     tCGAACAACGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACAccc                                                                                                                                                 >  1:4299628/1‑150 (MQ=255)
      cGAACAACGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCa                                                                                                                                                <  2:4642442/150‑1 (MQ=255)
        aacaacGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGc                                                                                                                                              <  1:4025822/150‑1 (MQ=255)
         acaacGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCg                                                                                                                                             >  1:4392960/1‑150 (MQ=255)
            cgcccTGGCCTGCCAGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTg                                                                                                                                         >  1:653552/4‑150 (MQ=255)
                cTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCa                                                                                                                                      <  1:1527988/150‑1 (MQ=255)
                 tGGCAGGCCAGGGCGGTCAGCGGAGTCGTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAg                                                                                                                                     <  1:4872922/150‑1 (MQ=255)
                  ggCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAgg                                                                                                                                    <  1:526359/150‑1 (MQ=255)
                   gCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGt                                                                                                                                   >  1:4897659/1‑150 (MQ=255)
                   gCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGt                                                                                                                                   >  2:3975365/1‑150 (MQ=255)
                      ctgccaggCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTGTc                                                                                                                                >  2:418855/7‑150 (MQ=255)
                       gCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATCCCACCGAGTTGCCGGCCACCCGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCa                                                                                                                               <  1:3322048/150‑1 (MQ=255)
                        ccAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCAc                                                                                                                              >  2:2796757/1‑150 (MQ=255)
                        ccAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCAc                                                                                                                              <  1:4488194/150‑1 (MQ=255)
                        ccAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCAc                                                                                                                              <  1:3253516/150‑1 (MQ=255)
                         cAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACg                                                                                                                             <  1:5227802/150‑1 (MQ=255)
                            ggCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCAt                                                                                                                          >  2:4934489/1‑150 (MQ=255)
                                    gCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACACATGTGGTGGc                                                                                                                  >  1:4156800/1‑150 (MQ=255)
                                            ttctCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGgcgc                                                                                                          >  1:4830637/1‑150 (MQ=255)
                                            ttctCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGgcgc                                                                                                          >  1:4254694/1‑150 (MQ=255)
                                             tctcGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCt                                                                                                         <  2:4659952/150‑1 (MQ=255)
                                             tctcGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCt                                                                                                         <  2:1106446/150‑1 (MQ=255)
                                                 ggACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCACCCGGGAAGGTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCAGGTGGTGGCCCGGGCGCTCGGa                                                                                                     >  2:4423869/1‑150 (MQ=255)
                                                 ggACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGa                                                                                                     >  1:2208579/1‑150 (MQ=255)
                                                  gACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGaa                                                                                                    >  2:3317191/1‑150 (MQ=255)
                                                         tCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATgg                                                                                             >  1:1373024/1‑150 (MQ=255)
                                                         tCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATgg                                                                                             <  2:1355822/150‑1 (MQ=255)
                                                           cATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTc                                                                                           <  1:2332331/150‑1 (MQ=255)
                                                             tACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAg                                                                                         >  1:986141/1‑150 (MQ=255)
                                                                caccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCa                                                                                                                               <  1:1678246/109‑1 (MQ=255)
                                                                caccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCa                                                                                                                               >  2:1678246/1‑109 (MQ=255)
                                                                caccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGcc                                                                                      >  2:3720556/1‑150 (MQ=255)
                                                                caccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGcc                                                                                      >  1:5309993/1‑150 (MQ=255)
                                                                 accGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCCACCGGGCAGGTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCg                                                                                     >  2:1570652/1‑150 (MQ=255)
                                                                  ccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGt                                                                                    <  1:4983923/150‑1 (MQ=255)
                                                                  ccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGGAGAAGGCGCTGATCACGCCGACCACACCCAGCGGGGGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGt                                                                                    >  2:4955195/1‑150 (MQ=255)
                                                                       ttGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGc                                                                               >  1:3641628/1‑150 (MQ=255)
                                                                          ccGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTgg                                                                            >  1:2363999/1‑150 (MQ=255)
                                                                           cGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGgc                                                                           >  1:4576967/1‑150 (MQ=255)
                                                                               caccaGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGa                                                                       >  2:4382082/1‑150 (MQ=255)
                                                                                 ccagggccagggCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACa                                                                     >  1:89871/1‑150 (MQ=255)
                                                                                    gggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGc                                                                  <  2:2174363/150‑1 (MQ=255)
                                                                                     ggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGtggtgg                                                                                                                   >  2:4891520/1‑100 (MQ=255)
                                                                                       ccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCaa                                                               <  2:2431702/150‑1 (MQ=255)
                                                                                         agggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCa                                                                                                                               <  1:1243206/84‑1 (MQ=255)
                                                                                         agggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCa                                                                                                                               >  2:1243206/1‑84 (MQ=255)
                                                                                             cGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGa                                                         <  2:2098174/150‑1 (MQ=255)
                                                                                              ggTGTTCCAGGCCCACACGGCAACCGGGACGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGag                                                        >  2:5249261/1‑149 (MQ=255)
                                                                                                  ttCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCg                                                    >  2:3548848/1‑150 (MQ=255)
                                                                                                  ttCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCg                                                    <  1:2715102/150‑1 (MQ=255)
                                                                                                    ccAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCa                                                  <  1:4846146/150‑1 (MQ=255)
                                                                                                    ccAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACACATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCa                                                  <  2:4156800/150‑1 (MQ=255)
                                                                                                       ggCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGAt                                               <  2:4847522/150‑1 (MQ=255)
                                                                                                          ccACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTc                                            >  2:3030600/1‑150 (MQ=255)
                                                                                                           cacaCGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCg                                           <  1:1728577/150‑1 (MQ=255)
                                                                                                           cacaCGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCg                                           <  2:2667226/150‑1 (MQ=255)
                                                                                                              ccGGCACCCGGGAAGTTGAAGGCGCTGATCACGCCGACCACCCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCAGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATc                                        <  1:4261977/149‑1 (MQ=255)
                                                                                                                   aaCCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAgg                                                                                        >  2:607938/1‑97 (MQ=255)
                                                                                                                     ccGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCCGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTccc                                 <  1:1650470/150‑2 (MQ=255)
                                                                                                                     ccGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCt                                 <  1:3380345/150‑1 (MQ=255)
                                                                                                                     ccGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCt                                 <  2:772182/150‑1 (MQ=255)
                                                                                                                     ccGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCt                                 <  1:2382510/150‑1 (MQ=255)
                                                                                                                      cGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTg                                <  2:2760025/150‑1 (MQ=255)
                                                                                                                       gggCAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGc                               >  2:3243942/1‑150 (MQ=255)
                                                                                                                       gggAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGc                               >  1:4220838/1‑150 (MQ=255)
                                                                                                                       gggAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGc                               <  2:5387079/150‑1 (MQ=255)
                                                                                                                       gggAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGc                               >  2:4041062/1‑150 (MQ=255)
                                                                                                                        ggAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGaa                                                        >  1:2754415/1‑124 (MQ=255)
                                                                                                                        ggAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCa                              <  1:4858192/150‑1 (MQ=255)
                                                                                                                        ggAAGTTGAAGGCCGAGATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTgta                              >  2:4628706/1‑148 (MQ=255)
                                                                                                                           aGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACtt                           <  1:2751480/150‑1 (MQ=255)
                                                                                                                               gAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCAcc                       <  2:5427751/150‑1 (MQ=255)
                                                                                                                                aaGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCAccc                      <  1:2994507/150‑1 (MQ=255)
                                                                                                                                  ggCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATc                                        >  2:4605290/1‑130 (MQ=255)
                                                                                                                                    ccgagATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGc                  >  2:463935/5‑150 (MQ=255)
                                                                                                                                     gcTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGcc                 >  2:896222/1‑150 (MQ=255)
                                                                                                                                     gcTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGcc                 >  1:1591096/1‑150 (MQ=255)
                                                                                                                                      gagATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCt                <  2:2002137/148‑1 (MQ=255)
                                                                                                                                        gATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTc              >  1:5267297/1‑150 (MQ=255)
                                                                                                                                        gATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTc              <  1:314620/150‑1 (MQ=255)
                                                                                                                                        gATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTc              <  1:4621637/150‑1 (MQ=255)
                                                                                                                                          tCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCt            >  2:2313092/1‑150 (MQ=255)
                                                                                                                                            aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg          >  1:2461401/1‑150 (MQ=255)
                                                                                                                                            aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg          <  1:2796757/150‑1 (MQ=255)
                                                                                                                                            aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg          >  1:4722022/1‑150 (MQ=255)
                                                                                                                                            aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg          <  2:4851441/150‑1 (MQ=255)
                                                                                                                                            aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg          <  2:2208579/150‑1 (MQ=255)
                                                                                                                                              gCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGt        <  2:5267297/150‑1 (MQ=255)
                                                                                                                                               ccGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTGGGTg       <  1:3866624/150‑1 (MQ=255)
                                                                                                                                                 gACCACACCCAGCGGGTGCCAGGTTTCACGCAAGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACGTCACCCAGGAAGGCCTGGGTGaa     <  2:125621/150‑2 (MQ=255)
                                                                                                                                                  aCCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATc    <  2:1264043/150‑1 (MQ=255)
                                                                                                                                                   ccACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCt   >  2:2761846/1‑150 (MQ=255)
                                                                                                                                                    cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt  <  1:418855/150‑1 (MQ=255)
                                                                                                                                                    cacaCCCAGCGGGTGCCAGGGTTCACGCAGGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACCGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTGCCGGGGTGATCtt  >  2:5219754/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTTTTCGAACAACGCCTGGCAGGCCAGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTT  >  NC_002947/6005341‑6005638

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