Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,492,055 A→C 5.4% intergenic (+115/+15) yhdZ → / ← degS amino acid ABC transporter ATP‑binding protein/serine endoprotease DegS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,492,0550AC 5.4% 643.8 / 18.3 201intergenic (+115/+15)yhdZ/degSamino acid ABC transporter ATP‑binding protein/serine endoprotease DegS
Reads supporting (aligned to +/- strand):  ref base A (126/56);  new base C (7/4);  total (139/62)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

AAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTACCTGGTT  >  NC_002947/1491911‑1492204
                                                                                                                                                |                                                                                                                                                     
   ttgttCCTCAGCCCGATCCTGCATTGAGGCTGTTTGGGCTGCGCGGCCGCCCTGCGCGGGCACGCCCGCTCCCAGCGGTTTAGCGCGGGCTTTGGCGCCTGCGTAATCCCTGTGGGCGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGcc                                                                                                                                               <  1:2989134/150‑1 (MQ=255)
       cccccccccAGACCCTCCATTGATGCTGTTGGGGCTGCCCGGCAGCCCTTCGCGGGCACGCCCGCGCCCACCGTTTTAGCGCTGGCTTCGGAGCCGCCGTACTGCCTGGGGGCGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTgg                                                                                                                                           <  1:4616138/143‑1 (MQ=255)
          tCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGtt                                                                                                                                        >  2:2293283/1‑150 (MQ=255)
             gCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGGTTAGCGCAGGCTTTGGAGCCGGCGGAAATCCGGTGGGAGCGGGCGTGCCCCCGCAGGGCCGCCACGCGGCCCGGGGGTTAt                                                                                                                                     >  1:5378492/1‑150 (MQ=25)
              ccAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTCGCGCTGGCGTTGGAGCCTGCGTCGTTCCTGTGGGCGCGGGCGTGCGCGCGAAGGGCCGCACCGCGGCGCTGGGGTTCtt                                                                                                                                    >  2:4255351/1‑150 (MQ=38)
                      tGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCGGGGTTTATTTCTCTTcc                                                                                                                            >  2:5399040/1‑150 (MQ=255)
                                                     ccTTGGCGGGCCCGCCCGCCCCCAGTGTTTTCGCGCGGGCTTTGGCGCTTGCTTATTTCCGGTGGGTGGGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCa                                                                                             <  1:1426236/150‑1 (MQ=255)
                                                            GGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCGGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCAc                                                                                      <  2:6911319‑M1/150‑1 (MQ=255)
                                                            GGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCGTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCGGGGTTTATTTCTCTTCCTGGCTGGCCGCCGGCGCCGGGGGCGGCCGCAGCCCCAc                                                                                      <  2:4285659‑M1/150‑1 (MQ=255)
                                                                acgcccgctcccacaggaattaatcaggctccaaagacagagctaaaactgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTATCTTACTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCCGCCCCACCACg                                                                                  >  1:950342‑M2/78‑150 (MQ=255)
                                                                  GCCCGCTCCCACTGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGc                                                                                >  2:5846015‑M1/1‑150 (MQ=255)
                                                                  GCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGc                                                                                <  1:3135580‑M1/150‑1 (MQ=255)
                                                                     cgctcccacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTAttt                                                                                                                                   <  2:4426272‑M2/24‑1 (MQ=255)
                                                                     cgctcccacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTAttt                                                                                                                                   >  1:4426157‑M2/73‑96 (MQ=255)
                                                                     CGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGAt                                                                             <  1:3511294‑M1/150‑1 (MQ=255)
                                                                       ctcccacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcTGCCGCGCAGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCa                                                                           >  2:2292957‑M2/71‑150 (MQ=255)
                                                                       ctcccacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcTGCCGCGCAGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCa                                                                           >  2:6044033‑M2/71‑150 (MQ=255)
                                                                       CTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCa                                                                           >  2:5315371‑M1/1‑150 (MQ=255)
                                                                       CTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCa                                                                           >  2:5315598‑M1/1‑150 (MQ=255)
                                                                        tcccacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCACCGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAg                                                                          >  2:1138115‑M2/70‑150 (MQ=255)
                                                                        tcccacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAg                                                                          >  2:5712464‑M2/70‑150 (MQ=255)
                                                                        tcccacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAg                                                                          >  1:5036895‑M2/70‑150 (MQ=255)
                                                                        tcccacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAg                                                                          >  2:5133927‑M2/70‑150 (MQ=255)
                                                                        TCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAg                                                                          <  1:4369368‑M1/150‑1 (MQ=255)
                                                                          ccacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCCACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt                                                                        >  1:3296381‑M2/68‑150 (MQ=255)
                                                                          ccacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt                                                                        >  1:5959445‑M2/68‑150 (MQ=255)
                                                                          CCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt                                                                        <  1:3890843‑M1/150‑1 (MQ=255)
                                                                          CCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt                                                                        <  2:2789255‑M1/150‑1 (MQ=255)
                                                                          CCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGAGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt                                                                        <  1:5302823‑M1/150‑1 (MQ=255)
                                                                           cacaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcTGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt                                                                       >  1:471158‑M2/67‑150 (MQ=255)
                                                                           CACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt                                                                       <  1:519872‑M1/150‑1 (MQ=255)
                                                                            acaggaattacgcaggctccaaagccagcgtaattcctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGttt                                                                                                                                       <  2:2184489‑M2/20‑1 (MQ=255)
                                                                            acaggaattacgcaggctccaaagccagcgtaattcctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGttt                                                                                                                                       >  1:2184489‑M2/66‑85 (MQ=255)
                                                                            acaggaattacgcaggctccaaagccagcgtaattcctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTc                                                                      >  1:4164485‑M2/66‑150 (MQ=255)
                                                                            acaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcTGCCGCGCAGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCgg                                                                                                              <  2:5447665‑M2/45‑1 (MQ=255)
                                                                            acaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTc                                                                      >  1:2831478‑M2/66‑150 (MQ=255)
                                                                            acaggaattacgcaggctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTc                                                                      >  1:4534076‑M2/66‑150 (MQ=255)
                                                                            ACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTc                                                                      <  1:2121684‑M1/150‑1 (MQ=255)
                                                                              AGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCaa                                                                    <  1:7303772‑M1/150‑1 (MQ=255)
                                                                               GGTTTTGCGCTGGCTTTGGCGCCTGCTTATTGCCTGTGGGTGGGGGGGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAAc                                                                   >  1:6945757‑M1/1‑150 (MQ=255)
                                                                                  attacgcaggctccaaagccagcgtaattcctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAActgc                                                                >  2:3707444‑M2/60‑150 (MQ=255)
                                                                                       gcTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGccc                                                           >  2:6837052/1‑150 (MQ=255)
                                                                                       gcTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGccc                                                           >  1:6748067/1‑150 (MQ=255)
                                                                                        cTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGATGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCg                                                          >  2:1105338/1‑150 (MQ=255)
                                                                                           gCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTg                                                       >  2:2248186/1‑150 (MQ=255)
                                                                                            ctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcTGCCGCGCAGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTgc                                                      <  1:2292957‑M2/101‑1 (MQ=255)
                                                                                            ctccaaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcTGCCGCGCAGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTgc                                                      <  1:6044033‑M2/101‑1 (MQ=255)
                                                                                             tttGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTgcg                                                     >  1:5122403/1‑150 (MQ=255)
                                                                                             tttGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTgcg                                                     >  2:3579665/1‑150 (MQ=255)
                                                                                             tttGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTgcg                                                     >  2:4413258/1‑150 (MQ=255)
                                                                                             tttGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTgcg                                                     >  1:5639520/1‑150 (MQ=255)
                                                                                             tttGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCGGCCCGTTgcg                                                     >  1:1719101/1‑150 (MQ=255)
                                                                                               caaaGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCa                                                   >  2:1853739/4‑150 (MQ=255)
                                                                                                ggAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCAt                                                  >  1:5628903/1‑150 (MQ=255)
                                                                                                aaagccagcgctaaacctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCAt                                                  >  2:4119789‑M2/46‑150 (MQ=255)
                                                                                                    ccTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCAcc                                              >  1:4974547/1‑150 (MQ=255)
                                                                                                    ccTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCAcc                                              <  1:5052685/150‑1 (MQ=255)
                                                                                                      tGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTc                                            >  2:2010543/1‑150 (MQ=255)
                                                                                                       gCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAg                                                                                            <  1:6498413/101‑1 (MQ=255)
                                                                                                        cGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGa                                          >  2:1345870/1‑150 (MQ=255)
                                                                                                          tatttCCTGGGGGCGCGGGCGGGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATg                                        <  1:1112302/147‑1 (MQ=255)
                                                                                                          tAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATg                                        <  2:6330913/150‑1 (MQ=255)
                                                                                                          tAATTCCTGGGGGAGGGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATg                                        <  1:5471269/150‑1 (MQ=255)
                                                                                                             aacctgtgggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGtt                                     <  1:1894806‑M2/118‑1 (MQ=255)
                                                                                                                 tgtgggagcgggcgtgccctcgaagggcCGCAACGCGGCCGTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTAtct                                 >  2:2442689‑M2/29‑150 (MQ=255)
                                                                                                                     ggggggggcgtgcccgcgacgggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTtctc                             <  1:5619311‑M2/126‑1 (MQ=255)
                                                                                                                     ggagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTtctc                             >  2:1573553‑M2/25‑150 (MQ=255)
                                                                                                                      gagcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGAGGGTTATCTTCTCg                            >  2:6386266‑M2/24‑150 (MQ=255)
                                                                                                                       cgcgggcacgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCGTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCCGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCGTCTCGt                           >  2:7165929‑M2/23‑150 (MQ=255)
                                                                                                                        gcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGtt                          <  1:6108670‑M2/129‑1 (MQ=255)
                                                                                                                        gcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGtt                          <  1:5860728‑M2/129‑1 (MQ=255)
                                                                                                                        gcgggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGtt                          >  1:5685971‑M2/22‑150 (MQ=255)
                                                                                                                          gggcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTgg                        >  2:3494259‑M2/20‑150 (MQ=255)
                                                                                                                           ggcgtgcccgcgaagggcCGCAACGCGGCCCGGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGGTGCGCATCACCTCGATGGTTATCGTCTCGTTggg                       >  1:3488362‑M2/19‑150 (MQ=255)
                                                                                                                            gcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGc                      >  2:695514‑M2/18‑150 (MQ=255)
                                                                                                                            gcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGc                      >  2:2143226‑M2/18‑150 (MQ=255)
                                                                                                                            gcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGc                      <  1:2010543‑M2/133‑1 (MQ=255)
                                                                                                                            gcgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCGCGATGGTTATCTTCACGTTGGGc                      >  2:7443987‑M2/18‑150 (MQ=255)
                                                                                                                             cgtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGATGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCt                     <  1:1105215‑M2/134‑1 (MQ=255)
                                                                                                                              gtgcccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCtt                    >  2:5385700‑M2/16‑150 (MQ=255)
                                                                                                                               tgcccgcgaagggcCGCAACGCGGCCCGGGGTTTAGTTCTCGTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGGTGCGCATCACCGCGATGGGTATCTGCGCGTGGGGCTTg                   >  2:648884‑M2/15‑150 (MQ=255)
                                                                                                                                 taaacccagggcCGCGTTGCGGCCCTTGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGAt                 >  2:7034327‑M2/13‑150 (MQ=255)
                                                                                                                                 cccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGAt                 >  2:3912387‑M2/13‑150 (MQ=255)
                                                                                                                                  ccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATc                <  2:4798915‑M2/139‑1 (MQ=255)
                                                                                                                                  ccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATc                >  2:3177957‑M2/12‑150 (MQ=255)
                                                                                                                                  ccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATc                >  2:6519482‑M2/12‑150 (MQ=255)
                                                                                                                                  ccgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGAGCAGCTTCAACTGCTGCCCGTTGCGCATCCCCTCGATGGTTAGCTTCGCGTGGGGCTGGATc                >  1:4774464‑M2/12‑150 (MQ=255)
                                                                                                                                   cgcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCTTTGGGCTTGATcc               >  1:7409774‑M2/11‑150 (MQ=255)
                                                                                                                                    gcgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCGCGTTGGGCTTGATCcg              >  2:5617884‑M2/10‑150 (MQ=255)
                                                                                                                                     cgaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCcgc             <  2:7303772‑M2/142‑1 (MQ=255)
                                                                                                                                     cgaagggcCGCAACGCCGCCCTGGGTATATATCTCTTCCTGCCTGGCAGCCGGCGCCGGTGGCGGCCGCACCCCCACCTCGGCGATCAGCTACAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGTGCTTGATCcgc             <  2:5044454‑M2/142‑1 (MQ=255)
                                                                                                                                     cccagggcCGCGTTGCGGCCCTTGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTgcgc                                                    >  2:4646022‑M2/9‑111 (MQ=255)
                                                                                                                                      gaagggcCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCcgcg            >  1:2095439‑M2/8‑150 (MQ=255)
                                                                                                                                          ggcCGCAACGCGGCCCTGGGTTGATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTAc        >  2:1010856‑M2/4‑150 (MQ=255)
                                                                                                                                            cCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGc                                                                                           >  1:5440461‑M2/2‑65 (MQ=255)
                                                                                                                                            cCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTACCt      >  2:2723662‑M2/2‑150 (MQ=255)
                                                                                                                                            cCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTACCt      >  2:3584623‑M2/2‑150 (MQ=255)
                                                                                                                                            cCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTACCt      >  2:1305948‑M2/2‑150 (MQ=255)
                                                                                                                                            cCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTACCt      >  2:4909308‑M2/2‑150 (MQ=255)
                                                                                                                                              gCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGTCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTACCTgg    >  1:1833256/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
AAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTACCTGGTT  >  NC_002947/1491911‑1492204

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