Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,543,444 G→A 7.7% L29L (CTG→TTG) ‡ 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,4440GA 7.7% 195.1 / 5.0 65L29L (CTG→TTG) ‡cysPthiosulfate/sulfate ABC transporter periplasmic binding protein CysP
Reads supporting (aligned to +/- strand):  ref base G (28/32);  new base A (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 = 8.49e-01

GTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAG  >  NC_000913/2543297‑2543593
                                                                                                                                                   |                                                                                                                                                     
gTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTcagca                                                                                                                                                     <  2:633449/150‑1 (MQ=255)
gTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTcagca                                                                                                                                                     >  1:797878/1‑150 (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                                                                                             <  2:18270/150‑1 (MQ=255)
                                                        ttGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCgagtgag                                                                                             <  1:6102/150‑1 (MQ=255)
                                                              ccgccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGCCCAGCGCGAGTGAGTTCttt                                                                                       >  1:138894/1‑150 (MQ=255)
                                                                 ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGACCTGTTCAGCAGTTCCGTTGCCTGTACATGGCCGGCCAGCAGCAGAGAAGCGACCAGCGCGAGGGAGTTCTTTTTc                                                                                    <  2:1021706/150‑1 (MQ=255)
                                                                 ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                    <  1:446182/150‑1 (MQ=255)
                                                                 ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                    <  2:1845106/150‑1 (MQ=255)
                                                                 ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                    <  1:172426/150‑1 (MQ=255)
                                                                 ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                    <  1:1376907/150‑1 (MQ=255)
                                                                 ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                    >  1:1798715/1‑150 (MQ=255)
                                                                             gCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACg                                                                        >  2:1088579/1‑150 (MQ=255)
                                                                             gCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACg                                                                        <  1:652624/150‑1 (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                                      >  2:1793061/1‑150 (MQ=255)
                                                                                                               cAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTc                                      >  1:775196/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)
                                                                                                                                                  agcagTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGa   >  1:1207153/1‑150 (MQ=255)
                                                                                                                                                   gcagTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAg  >  1:1915299/1‑150 (MQ=255)
                                                                                                                                                   gcagTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAg  >  1:865621/1‑150 (MQ=255)
                                                                                                                                                   gcagTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAg  <  2:104929/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAG  >  NC_000913/2543297‑2543593

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: