Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,781 Δ1 bp 6.0% intergenic (+551/+367) PP_1865 → / ← aroH transposase/3‑deoxy‑7‑phosphoheptulonate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,087,7810G. 6.0% 296.3 / 13.4 100intergenic (+551/+367)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base G (52/42);  new base . (3/3);  total (55/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 101 total aligned reads displayed.

GCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACAT  >  NC_002947/2087634‑2087930
                                                                                                                                                   |                                                                                                                                                     
gctcttgctcttGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCA                                                                                                                                                     >  1:1419987‑M1/1‑150 (MQ=255)
  tcttgctcttgcTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGC                                                                                                                                                   >  2:86895‑M1/1‑150 (MQ=255)
    ttgctcttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGC                                                                                                                                                 <  1:1776133‑M1/150‑1 (MQ=255)
         cttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGAGGGGCATGGATGCCCGTCAAGCGCTAGGCCCCATCCTGGGGCC‑‑CAGCGCttgacgg                                                                                                                                          <  1:1594386‑M1/150‑8 (MQ=255)
            gCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCctcc                                                                                                                                         <  1:1515602‑M2/150‑5 (MQ=255)
             ctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTccc                                                                                                                                        >  2:1727666/1‑150 (MQ=32)
             ctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTccc                                                                                                                                        >  1:748989/1‑150 (MQ=32)
              tttGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCg                                                                                                                                       >  1:448690/1‑150 (MQ=32)
               ttGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCgg                                                                                                                                      <  1:882617/150‑1 (MQ=32)
                tGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCggg                                                                                                                                     >  2:1412000/1‑150 (MQ=32)
                tGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCggg                                                                                                                                     >  1:202209/1‑150 (MQ=32)
                  gCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGCAGCTGTCTGCAGCGCGGTCCTCCCGGGAg                                                                                                                                   >  1:345/1‑150 (MQ=11)
                    ttttGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCa                                                                                                                                 >  2:1665668/1‑150 (MQ=32)
                    ttttGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCA                                                                                                                                                     <  2:256239‑M1/130‑1 (MQ=255)
                    ttttGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCA                                                                                                                                                     >  1:256239‑M1/1‑130 (MQ=255)
                        gCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGcc                                                                                                                             <  1:1412000/150‑1 (MQ=32)
                         cTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGGCATGGGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGccc                                                                                                                            <  2:1441728/150‑1 (MQ=12)
                          ttCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCg                                                                                                                           >  2:725745/1‑150 (MQ=32)
                          ttCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCg                                                                                                                           >  1:1068724/1‑150 (MQ=32)
                          ttCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCg                                                                                                                           >  1:1292192/1‑150 (MQ=32)
                             tAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCctccgggcttgctcccgggag                                                                                                                        >  2:1503723‑M2/1‑129 (MQ=255)
                             tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCctcccgggcttgctcccggga                                                                                                                        >  2:1230762‑M2/1‑129 (MQ=255)
                               aGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTg                                                                                                                      >  1:1882852/1‑150 (MQ=32)
                                   cgcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAggg                                                                                                                  <  2:773148/150‑1 (MQ=32)
                                    gcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAgggg                                                                                                                 >  1:264501/1‑150 (MQ=32)
                                      ggTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAc                                                                                                               <  2:433329/150‑1 (MQ=32)
                                       gTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCA                                                                                                                                                     <  2:716666‑M1/111‑1 (MQ=255)
                                       gTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCA                                                                                                                                                     >  1:716666‑M1/1‑111 (MQ=255)
                                        tAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAccc                                                                                                             >  1:756613/1‑150 (MQ=32)
                                           ttCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGATGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCgg                                                                                                          >  1:1168955/1‑150 (MQ=17)
                                              aGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagc                                                                                                       >  2:1661824/1‑150 (MQ=32)
                                                gCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGGcccctacgcta                                                                                                     <  1:1562261/150‑8 (MQ=255)
                                                 cGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgta                                                                                                    <  1:1129235/150‑1 (MQ=32)
                                                                   cAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGgatgat                                                                                  <  2:1417730/150‑1 (MQ=11)
                                                                    aggaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGGTCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGatgatg                                                                                 <  1:1085993/150‑1 (MQ=14)
                                                                    aggaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGatgatg                                                                                 >  1:28192/1‑150 (MQ=31)
                                                                      gaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATggg                                                                               <  1:1665668/150‑1 (MQ=32)
                                                                       aggCCGACCGCAGGGCTTGCGGAGGGAGGTGACGGGAATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGTAGTATGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGa                                                                              <  1:242363/150‑1 (MQ=11)
                                                                         gCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGc                                                                            >  2:668800/1‑150 (MQ=32)
                                                                          ccGAGCGCATGGCTTGCGGAGGGAGGTGACGGGTATGCATGCCCGCCAAGCGCTGGGCCCTACGATGGGCCCTGCAGCGCGGTGCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGCGGCGTAGCGTTGGGGCCGGATGATGGGAGcc                                                                           <  2:1040952/150‑1 (MQ=11)
                                                                            gAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGa                                                                         <  2:1882852/150‑1 (MQ=32)
                                                                                   gggCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGtttt                                                                  <  1:725745/150‑1 (MQ=35)
                                                                                    ggCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTg                                                                 <  2:1515602/150‑1 (MQ=35)
                                                                                     gCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGCGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTgg                                                                >  2:162061/1‑150 (MQ=25)
                                                                                     gCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTgg                                                                <  2:1468069/150‑1 (MQ=34)
                                                                                       ttGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGtt                                                              >  1:1537132/1‑150 (MQ=34)
                                                                                               ggaggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGc                                                      <  2:1292192/150‑1 (MQ=35)
                                                                                               ggaggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGc                                                      <  2:345/150‑1 (MQ=35)
                                                                                                aaggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTACGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGcc                                                     >  2:1328488/2‑150 (MQ=21)
                                                                                                 aggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGc                                                      <  2:1073218/148‑1 (MQ=35)
                                                                                                 aggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGc                                                      >  1:1073218/1‑148 (MQ=34)
                                                                                                  cttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGa                                                                              <  2:228655/121‑1 (MQ=255)
                                                                                                  cttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGa                                                                              >  1:228655/3‑123 (MQ=255)
                                                                                                    tGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGAGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGa                                                 <  2:650312/150‑1 (MQ=25)
                                                                                                     gACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGAc                                                <  1:1661824/150‑1 (MQ=35)
                                                                                                     gACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGAc                                                >  1:396675/1‑150 (MQ=255)
                                                                                                     gACGGGCATCCATGCCAGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGAGGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGAc                                                >  1:305357/1‑150 (MQ=255)
                                                                                                        gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCctccgggcttgctcccgggagga                                                                                                                      <  2:642486‑M2/77‑24 (MQ=255)
                                                                                                        gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCctccgggcttgctcccgggagga                                                                                                                      >  1:642486‑M2/1‑54 (MQ=255)
                                                                                                        gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCa                                              <  1:763937/149‑1 (MQ=35)
                                                                                                        gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCa                                              >  2:763937/1‑149 (MQ=34)
                                                                                                          gCATGGATGCCCGTCAAGCGCTGGGACGAAGAATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGACGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAg                                           >  1:544141/1‑150 (MQ=11)
                                                                                                           cATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGt                                          >  2:783188/1‑150 (MQ=35)
                                                                                                             tGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACACTTTTGGTTGCTCTGGTCCCGACCATCGTaa                                        >  2:417607/1‑150 (MQ=11)
                                                                                                              ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAc                                       <  1:417607/150‑1 (MQ=35)
                                                                                                               gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTgg                                                                >  1:1052955/1‑124 (MQ=35)
                                                                                                               gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTgg                                                                <  2:1052955/124‑1 (MQ=37)
                                                                                                                 tGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGTAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccg                                    >  2:1126887/1‑150 (MQ=21)
                                                                                                                    ccGTCAAGCGCTTGGCCCCAGGATGGGACCT‑‑‑GCGCttgacgggcatccatgcccgtcacctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctggg                                                                 <  2:1462953‑M1/115‑81 (MQ=255)
                                                                                                                    ccGTCAAGCGCTTGGCCCCAGGATGGGACCT‑‑‑GCGCttgacgggcatccatgcccgtcacctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctggg                                                                 >  1:1462953‑M1/1‑35 (MQ=255)
                                                                                                                    ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGtttt                                                                  >  1:1897509/1‑117 (MQ=37)
                                                                                                                    ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGtttt                                                                  <  2:1897509/117‑1 (MQ=37)
                                                                                                                    ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                 >  2:1084826/1‑150 (MQ=35)
                                                                                                                    ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                 >  1:957632/1‑150 (MQ=34)
                                                                                                                     cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatctgcgtcgcctgaactaccgcgcacttagaagcaaaa                                          >  2:1161616‑M1/1‑37 (MQ=255)
                                                                                                                     cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatctgcgtcgcctgaactaccgcgcacttagaagcaaaa                                          <  1:1161616‑M1/140‑104 (MQ=255)
                                                                                                                       tCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCttgacgggcatccatgcccgtcacctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaa                              >  2:517944‑M1/1‑35 (MQ=255)
                                                                                                                       tCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaa                                                                                <  1:1242336‑M1/100‑66 (MQ=255)
                                                                                                                       tCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaa                                                                                >  2:1242336‑M1/1‑35 (MQ=255)
                                                                                                                         aaGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCGCttgacgggcatccatgcccgtcacctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgag                         >  1:1408427‑M1/1‑30 (MQ=255)
                                                                                                                          aGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAggg                           <  2:28192/150‑1 (MQ=35)
                                                                                                                           gcgcGGGGCCCCAGGGGGGGGCTGGCAGCGCGGCCCGCCGGGGAGCAGGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTACCCCGCCGGAAGGGc                          <  2:1419987/150‑1 (MQ=11)
                                                                                                                           gcgcGGGGCCCCAGGAGGGGGCCTGCAGCGCGGTCCGCCCGGGAGGCGGCCCTGAGTGAGGGGGCCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAGCCCGCCGGAAGGGc                          <  2:544141/150‑1 (MQ=11)
                                                                                                                               tGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGaa                      >  1:692038/1‑150 (MQ=35)
                                                                                                                               tGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGAAGAAAGGGCAGaa                      >  2:1564160/1‑150 (MQ=255)
                                                                                                                               tGGGCCCCAGGATGGGGCC‑‑CAGCGCttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatctgcgtcgcctgaactaccgc                                                          >  2:476280‑M1/1‑25 (MQ=255)
                                                                                                                               tGGGCCCCAGGATGGGGCC‑‑CAGCGCttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatctgcgtcgcctgaactaccgc                                                          <  1:476280‑M1/112‑88 (MQ=255)
                                                                                                                                   ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTCACCCGCCGGAAGGGCGGAAAGGt                  <  2:202209/150‑1 (MQ=25)
                                                                                                                                   ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGt                  >  1:773412/1‑150 (MQ=35)
                                                                                                                                    cccAGGATGGGGCCTGCAGCGCGGTCctccgggcttgctcccgggag                                                                                                                        <  1:1107624‑M2/47‑22 (MQ=255)
                                                                                                                                    cccAGGATGGGGCCTGCAGCGCGGTCctccgggcttgctcccgggag                                                                                                                        >  2:1107624‑M2/1‑26 (MQ=255)
                                                                                                                                    cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                 >  2:1200995/1‑150 (MQ=35)
                                                                                                                                     ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGAGGTGAGGGGCCCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                <  1:1381940/150‑1 (MQ=255)
                                                                                                                                      cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGaa               >  2:472485/1‑150 (MQ=35)
                                                                                                                                      cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGaa               <  1:1727666/150‑1 (MQ=35)
                                                                                                                                        ggATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGaa               >  1:333540/1‑148 (MQ=35)
                                                                                                                                        ggATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGaa               <  2:333540/148‑1 (MQ=35)
                                                                                                                                                ccTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTa     <  2:1168955/150‑1 (MQ=35)
                                                                                                                                                  tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAATAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACa   >  2:741421/1‑150 (MQ=38)
                                                                                                                                                   gCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACAt  <  1:668800/150‑1 (MQ=38)
                                                                                                                                                   |                                                                                                                                                     
GCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACAT  >  NC_002947/2087634‑2087930

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