Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,543,442 C→G 7.7% L29L (CTG→CTC) ‡ cysP ← thiosulfate/sulfate ABC transporter periplasmic binding protein CysP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,543,4420CG 7.7% 212.4 / 7.6 65L29L (CTG→CTC) ‡cysPthiosulfate/sulfate ABC transporter periplasmic binding protein CysP
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base G (3/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

CCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGG  >  NC_000913/2543294‑2543591
                                                                                                                                                    |                                                                                                                                                     
ccTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTca                                                                                                                                                      <  1:694416/150‑1 (MQ=255)
ccTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTca                                                                                                                                                      >  2:652624/1‑150 (MQ=255)
   gTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTcagca                                                                                                                                                   >  1:797878/1‑150 (MQ=255)
   gTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTcagca                                                                                                                                                   <  2:633449/150‑1 (MQ=255)
     aaaaTCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGt                                                                                                                                                 >  1:1795431/1‑150 (MQ=255)
         tCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCg                                                                                                                                             <  1:1819025/150‑1 (MQ=255)
           gCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGtt                                                                                                                                           >  1:1342239/1‑150 (MQ=255)
           gCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGCTGAACAGTTcttat                                                                                                                                           <  2:292232/150‑5 (MQ=255)
           gCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGCTGAACAGTTcttat                                                                                                                                           <  1:498208/150‑5 (MQ=255)
              aGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGcc                                                                                                                                        <  2:1795605/150‑1 (MQ=255)
              aGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGcc                                                                                                                                        <  1:1450325/150‑1 (MQ=255)
                 gcCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGt                                                                                                                                     >  2:709806/1‑150 (MQ=255)
                  ccTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTa                                                                                                                                    <  1:916225/150‑1 (MQ=255)
                     gTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACAt                                                                                                                                 >  1:1843739/1‑150 (MQ=255)
                      ttttGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATg                                                                                                                                >  2:1167538/1‑150 (MQ=255)
                               cccGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCcag                                                                                                                       <  2:797878/150‑1 (MQ=255)
                                  ggCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCcagcag                                                                                                                    >  2:1376907/1‑150 (MQ=255)
                                   gCAGGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCcagcagc                                                                                                                   <  1:1576218/150‑1 (MQ=255)
                                      tGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCaga                                                                                                                <  1:1441994/150‑1 (MQ=255)
                                       gagaTTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCagag                                                                                                               <  1:1585680/150‑1 (MQ=255)
                                          aTTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAg                                                                                                            <  2:1487587/150‑1 (MQ=255)
                                          aTTGTTTAATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACCTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGCCCCGCCAGCAGCAGAGaac                                                                                                            >  2:596459/1‑149 (MQ=255)
                                           ttGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGc                                                                                                           <  2:116677/150‑1 (MQ=255)
                                             gTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGa                                                                                                         >  1:1845106/1‑150 (MQ=255)
                                               tttATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGAcc                                                                                                       <  1:13671/150‑1 (MQ=255)
                                               tttATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGAcc                                                                                                       <  2:584057/150‑1 (MQ=255)
                                                  aTCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCCGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAgc                                                                                                    >  2:780662/1‑150 (MQ=255)
                                                  aTCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTACGTTGCCTGTACATGGCCCCCCAGCAGCAGAGAAGCGACCAgc                                                                                                    >  2:63526/1‑150 (MQ=255)
                                                    cGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAgcgc                                                                                                  >  1:1988351/1‑150 (MQ=255)
                                                    cGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAgcgc                                                                                                  >  1:1360190/1‑150 (MQ=255)
                                                           ttGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCgagtgag                                                                                           <  1:6102/150‑1 (MQ=255)
                                                           ttGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCgagtgag                                                                                           <  2:18270/150‑1 (MQ=255)
                                                                 ccgccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGCCCAGCGCGAGTGAGTTCttt                                                                                     >  1:138894/1‑150 (MQ=255)
                                                                    ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGACCTGTTCAGCAGTTCCGTTGCCTGTACATGGCCGGCCAGCAGCAGAGAAGCGACCAGCGCGAGGGAGTTCTTTTTc                                                                                  <  2:1021706/150‑1 (MQ=255)
                                                                    ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                  <  2:1845106/150‑1 (MQ=255)
                                                                    ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                  <  1:1376907/150‑1 (MQ=255)
                                                                    ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                  <  1:446182/150‑1 (MQ=255)
                                                                    ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                  >  1:1798715/1‑150 (MQ=255)
                                                                    ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                  <  1:172426/150‑1 (MQ=255)
                                                                                gCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACg                                                                      <  1:652624/150‑1 (MQ=255)
                                                                                gCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACg                                                                      >  2:1088579/1‑150 (MQ=255)
                                                                                      tgttgCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCAtt                                                                <  1:1776814/150‑1 (MQ=255)
                                                                                       gttgCTCACTCGGCGGATTCAGGGCGGCAAACAGCTCGGGGGAGACGTCATAAGAAGTGTTCAGCAGTTCCGTTGCCTGTTCATGGCCGGCGAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCGATTg                                                               <  2:1652571/150‑1 (MQ=255)
                                                                                           cTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCa                                                           >  1:208667/1‑150 (MQ=255)
                                                                                             cAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCAcc                                                         >  2:1924813/1‑150 (MQ=255)
                                                                                             cAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCAcc                                                         >  2:72510/1‑150 (MQ=255)
                                                                                               aaCGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCt                                                       <  2:766894/150‑1 (MQ=255)
                                                                                                    cggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATaa                                                  >  2:728273/1‑150 (MQ=255)
                                                                                                    cggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATaa                                                  <  1:1088579/150‑1 (MQ=255)
                                                                                                     ggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATaaa                                                 <  2:1174485/150‑1 (MQ=255)
                                                                                                          cAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTaa                                            >  2:1746359/1‑150 (MQ=255)
                                                                                                                  cAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTc                                    >  1:775196/1‑150 (MQ=255)
                                                                                                                  cAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTc                                    >  2:1793061/1‑150 (MQ=255)
                                                                                                                   aaaCAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCGGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCt                                   <  2:1342239/150‑1 (MQ=255)
                                                                                                                   aaaCAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCt                                   <  2:1798715/150‑1 (MQ=255)
                                                                                                                         cTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGc                             >  1:167824/1‑150 (MQ=255)
                                                                                                                           cgcgGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCca                           >  1:1242789/1‑150 (MQ=255)
                                                                                                                             cgGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCcatc                         <  2:1843739/150‑1 (MQ=255)
                                                                                                                                      tCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTtata                <  1:1038314/150‑1 (MQ=255)
                                                                                                                                       caacggAACTGCTGAACAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTCATAAATTTAATGACTTTCTAATAGCCATCATATTTATaa               >  1:1653688/6‑150 (MQ=255)
                                                                                                                                        aacggAACTGCTGAACAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAAc              >  1:782079/5‑150 (MQ=255)
                                                                                                                                        aTAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAAc              <  2:167824/150‑1 (MQ=255)
                                                                                                                                         tAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACg             >  1:162645/1‑150 (MQ=255)
                                                                                                                                         acggAACTGCTGAACAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACg             >  1:113986/4‑150 (MQ=255)
                                                                                                                                          aaGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACgg            >  1:1975431/1‑150 (MQ=255)
                                                                                                                                               cTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTg       >  1:559045/1‑150 (MQ=255)
                                                                                                                                                tGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGa      <  1:60193/150‑1 (MQ=255)
                                                                                                                                                 gTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGaa     >  1:1057409/1‑150 (MQ=255)
                                                                                                                                                   tCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAg   <  2:1651531/150‑1 (MQ=255)
                                                                                                                                                    cagcagTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAgg  >  1:1649123/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGG  >  NC_000913/2543294‑2543591

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