Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,984,496 A→C 18.4% V504G (GTG→GGG)  aslA ← putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,984,4960AC18.4% 133.1 / 4.1 62V504G (GTG→GGG) aslAputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base A (33/15);  new base C (4/10);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 1.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGC‑AGATCTCAGC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTT  >  NC_000913/3984281‑3984719
                                                                                                                                                                                                                              |                                                                                                                                                                                                                               
gttttggtggggtgggaaaGGTGCGTCCGCATGGCATCCGGCAAGGGC‑AGTTCCCAGT‑GATGGCTCCGGTTGAGTCAGTTTAAATC‑TGGGGGGGTGGTGG‑ATATTTTTCCAGGATCGCCA‑TATCCGCGTGCATTTCGGTCTG‑T‑AGCGGTACCCCCATCGGAATAGGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCccg                                                                                                                                                                                                                     <  1:434809/211‑1 (MQ=255)
  ttCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGC‑AGATCTCAGC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCccg                                                                                                                                                                                                                     >  2:588138/1‑226 (MQ=255)
  ttCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGC‑AGATCTCAGC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGGGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATcccc                                                                                                                                                                                                                     >  2:370066/1‑225 (MQ=255)
  ttCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGC‑AGATCTCAGC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGCGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGCGTCGCTTTCCTGCGGACCGGTGTAGAGGTTAAACACCGACGAACccgcc                                                                                                                                                                                                                   >  1:313587/1‑228 (MQ=255)
       tagcgtggaaaaggagggcaaGCATCGCAT‑CGGGAAAGGCGAGTGCCGAGC‑GATTGCGCCGGGTTAGTAAAATTTATTC‑GGCGCGCGGGGGGG‑TTATTTTTTCAGTACCCCCA‑TATCCGCGGGCATTTCGGTCTG‑T‑AGCGGTACCCCCACCGGAATAGGGGGCACGCCGATGGAGTCGCTTTCCTGAGGATCGGTGTAGAGGTTAAACACCGACGATACCGCCGTTTg                                                                                                                                                                                                              <  1:370066/208‑1 (MQ=255)
            tGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGC‑AGATCTCAGC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGCCGCACGCCGAGGGAGTGGCTTTCCTGCGGCTCGGGGTAAAGGTGAAAACCCGACGAGCCCGCCGtttacctac                                                                                                                                                                                                         >  1:225682/1‑222 (MQ=255)
            tGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGC‑AGATCTCAGC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGTGG‑ATATTTTCTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCCCACGCCGATGGAGTCGCTTTCCTGCGGTACGGCGTCAAGGGTAAAACCCGAGAAGCCCGCCGTTTGCAttc                                                                                                                                                                                                         >  1:589503/1‑227 (MQ=255)
               aTAAGATGCGTCAGCATCGCATCCGGCAAAGGC‑AGATCTCAGC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGGGTAGGGGTTAAAACCCGACGGACCCGCCCTTTTCATAACt                                                                                                                                                                                                       >  1:338954/1‑227 (MQ=255)
               aTAAGATGCGTCAGCATCGCATCCGGCAAAGGC‑AGATCTCAGC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGGGTAGAGGTTAAAACCCGCCGATCCCGCCGTTTGCATTCCtg                                                                                                                                                                                                      >  2:316955/1‑228 (MQ=255)
                    aTGCGTCAGCATCGCATCCGGCAAAGGC‑AGATCTCAGC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAAAACCGCCGATCCCGCCGTTTGCATTACTGTGcc                                                                                                                                                                                                  >  2:573965/1‑227 (MQ=255)
                        gtaagctccgactcgggaaaaGGA‑AGTTCTCGGG‑GATGGCGCTGGCTTAGCCAGTTTTAATC‑TGGGTGGGTGGTGG‑ATGTTTTTGCAGGATCTCCA‑TATCGGCGTGCATTTCGGTCTT‑G‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGa                                                                                                                                                                                             <  1:588138/210‑1 (MQ=255)
                                cgcttcCGGCAAAGGG‑ATCTCGCAGG‑GATGGCGCCGGTTTTTTCAGA‑TTAATCGTGCGCGGGTGGTGG‑ATTTTTTTTCAGCGTCTCAT‑TAGACGGGTGCGTTTAGTGTTG‑A‑GGCGTTACCCCCAGTGGAATAGGGCGCACGCCGAGGGGGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCAcc                                                                                                                                                                                      <  1:38036/227‑1 (MQ=255)
                                  cactcGGGAAAGGG‑GGATCTCCGGGGATGGCGCCGGTTTGGTCAGTTTTAATC‑TGGGGGGGGGGTGG‑ATATTTTTTCAGGATCCCAA‑AATCCGCGTGCATTTGGGTCTG‑G‑GGCGGTACCCCCATCGGAATATGGCCCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAACCCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCAcc                                                                                                                                                                                      <  1:573965/222‑1 (MQ=255)
                                                      tccggggttggggcCGGTTTGGTCAGATTTAACC‑TGGGGGCGTGGTGG‑TTTTTTTTTAGGGATCCCCA‑TTTCCGGGTGTCTTTGGGT‑TG‑TGAGCGGTACCCCCACGGGAATATGGCGCACGCCGATGGAGTCGCTTTCCGGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGtat                                                                                                                                                                  <  1:330256/215‑1 (MQ=255)
                                                         gC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTTCATTACTGTGCCGGTGAATCCACCCGGATATCCCCCCTGGGCAAAAGCg                                                                                                                                                             >  1:651733/1‑228 (MQ=255)
                                                             ttgggtccGGCTTAGTCAGATTTAATC‑TGGCAGCGTGGTGG‑ATATTTTTTCAGGATCCCCA‑TATACGCGTGAATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTa                                                                                                                                                           <  2:517224/221‑1 (MQ=255)
                                                               aGCCCCGGTTTAGTCAGTTTTAATC‑GGCGCGGGTGGTGG‑ATATTTTTTCGGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAg                                                                                                                                                         <  2:338954/227‑1 (MQ=255)
                                                                gccgcgttttgggcataTTAAATC‑TGCGGGCGTGGTGG‑ATTTTTTTACAGGTTCCCCA‑TATACCGGTGCTTGGTGTTCTG‑G‑AGGGGCACCCCCCTCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTaa                                                                                                                                                          <  1:532575/209‑1 (MQ=255)
                                                                 cgcCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGGGAATCCACCCTGATTTCCGCCCTGGGCTTAAGCGTAAGGt                                                                                                                                                       >  2:451300/1‑227 (MQ=255)
                                                                  tcgcgtttAGTCAAATTTAATC‑CGCGCGCGTGGTGG‑ATATTTTTTGAGGTTCTCCA‑TATCGGGGTGCATTGCGGTCTT‑T‑AGCGGTCCACCCAGCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTg                                                                                                                                                     <  2:188857/222‑1 (MQ=255)
                                                                  gcCGGTTTAGTCAGTTTTATTC‑GGGCGGCTTGGTGG‑ATATTTTTTC‑GGGTCTCCATTGTGCGGGAGCTTTTGGGTCTG‑T‑AGCGGTACACCCATCGGAATGTGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTg                                                                                                                                                     <  1:354588/228‑1 (MQ=255)
                                                                     taaaggacaaGATTTTATC‑TGGGGGGGGTGTGGAATATTTTTCCAGGTTCCCCA‑TATGCGCGGGCTTTTC‑GTGTGTT‑AGGGTTACCGCCTTGGAAATGTGGCGCACCCCGATGGGGTCGCTTTCCGGCGGATCGGTGTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCt                                                                                                                                                   <  1:116357/219‑1 (MQ=255)
                                                                                aTTTAATC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCCCTGGGGATAAGCGTAAGGTGGCTAAATCAGGACgg                                                                                                                                       >  1:316447/1‑227 (MQ=255)
                                                                                     aTC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGGTAATCCACCCTGCTATCCGCTCTGGGAATAAGGGTCAGGTTGCTGAATCAGGCCGTGaac                                                                                                                                   >  2:295345/1‑225 (MQ=255)
                                                                                           cgcgcgTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTCCTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGaa                                                                                                                             >  2:619751/1‑228 (MQ=255)
                                                                                           cgcgcgTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGc                                                                                                                                                    >  2:116357/1‑205 (MQ=255)
                                                                                           cgcgcgTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGAAATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGa                                                                                                                              >  2:169562/1‑227 (MQ=255)
                                                                                           cgcgcgTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGGGAATCCACCCTGATATCCGCTCTTGGTATAAGCGTAAGGTTGCTGAATACGGACGTGATACTTGaa                                                                                                                             >  2:247010/1‑228 (MQ=255)
                                                                                           cgcgcgTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGAGTAGAGGGTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGAGAATCCACCCTGATATCCGCTCTGGGTATACGCGTAAGGTTGCTGAATCAGAACGGGCAACTTGa                                                                                                                              >  2:399315/1‑227 (MQ=255)
                                                                                               cgTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGAATAAGCGTAAGGTTGCTGAATCAGGACGTGATCCTTGAACTct                                                                                                                         >  2:576399/1‑227 (MQ=255)
                                                                                                          atTTTTTCAGGATCCCCA‑TATGCGCTTGCATTTCGGTCTG‑T‑AGCGGTACACCCATGGGAATATGGCGCAGGCCGATGGAGTCGCTTTCCGGGGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGc                                                                                                                <  2:492098/227‑1 (MQ=255)
                                                                                                                   aggatccccattagccggTGCATTTCGGTCTG‑G‑AGCTGTACCCCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCACCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACagcagc                                                                                                        <  2:612982/210‑1 (MQ=255)
                                                                                                                         cccatgtccggggtcATTTCGGTCTG‑T‑ACGGGTCCCCCCAGCGGAATAGGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGGTTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTa                                                                                                 <  2:315041/214‑1 (MQ=255)
                                                                                                                             tataCGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGGATCAGGACGTGATACTTGAACTCATCCAAACGCACCGAAGCGAGTTTCCCGt                                                                                             >  1:511317/1‑228 (MQ=255)
                                                                                                                                   cgTGCATTTCGGTCTG‑T‑AGCGGTACACCCCTCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCGGATATCCGCTCTGGGTCTCAGCGGAAGGGTTCTGAATCAGGACATGATACTTGAACTCATACATACGCCACGCGGCAAGGTTACCGGTGAGGa                                                                                       >  1:148073/1‑228 (MQ=255)
                                                                                                                                    tgtctctcttggggttgtgggggTACCCCAATGGGAATGTGGGGCACGCCGGTGGGGTCCGTTGCCTGCGGTGCGGGGTAGAGTTTAACCCCCGCCGTTCCCGCCGTTTGCAGTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAgg                                                                                        <  1:545055/207‑1 (MQ=255)
                                                                                                                                     tGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGaa                                                                                                                                                >  2:20539/1‑169 (MQ=255)
                                                                                                                                     tGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGCCGTGATACTTGAACTCATCCATACGAACAGCCGCGAGTTTACCGTTGAGGaaa                                                                                     >  1:315041/1‑227 (MQ=255)
                                                                                                                                     tGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATCAGCACAGCAGCGAGTTTACCGTTGAGGaa                                                                                      >  2:326664/1‑227 (MQ=255)
                                                                                                                                     tGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAACCGTAAGGTTGCTGAACCAGGACGTGATACTTGAACTCATCCATACGCCAAGCAGGAAGTTTACAGTTGAGGaa                                                                                      >  1:334773/1‑227 (MQ=255)
                                                                                                                                     tGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCCCTCTGGGTATAAGCGTAAGGTTGCTGAATTAGGACGGGATACTTGAACTCATCCATCCGCACAGCACCGCGTTTACCGTTGAGGag                                                                                      >  1:646201/1‑226 (MQ=255)
                                                                                                                                       cATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTATAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGACTCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTACCCGTTGGGGAagta                                                                                   >  2:182698/1‑228 (MQ=255)
                                                                                                                                            cggtcgg‑t‑AGCGGTACCCCCATGGGAATATGGCGCACGCCGATGGAGTCGCTTTCCGGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCt                                                                              <  2:562756/228‑1 (MQ=255)
                                                                                                                                                    t‑AGGGGTACACCAATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCt                                                                        <  2:316447/227‑1 (MQ=255)
                                                                                                                                                    ggtgggttccccctatgggAATGTGGCGCACCCCGATGGGGTCCCTTTCCTGCGGATCGGTGTAGAGGTTAAACCCCGCCGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCt                                                                        <  1:169562/211‑1 (MQ=255)
                                                                                                                                                     agggggTCCACCCATCGGAATAGGGCGAACGCCTATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCt                                                                        <  2:252872/223‑1 (MQ=255)
                                                                                                                                                     aggcgttccccccATCGGAATTTGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCt                                                                        <  2:334773/217‑1 (MQ=255)
                                                                                                                                                          gTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGt                                                                  >  2:166274/1‑228 (MQ=255)
                                                                                                                                                                              cgcACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGAAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTccc                                              >  2:2903/1‑228 (MQ=255)
                                                                                                                                                                                 aCGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGt                                                                                             >  1:74840/1‑178 (MQ=255)
                                                                                                                                                                                 aCGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGt                                                                                             <  2:74840/178‑1 (MQ=255)
                                                                                                                                                                                            gTCGCTTTCCCGCGGAGGGGTTTAGGGGTTAACACCCGACGTCCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGt                                <  1:166274/228‑1 (MQ=255)
                                                                                                                                                                                               gCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGc                                                                                                                                                                           >  2:619878/1‑86 (MQ=255)
                                                                                                                                                                                                cTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGct                                                                                                                                                                          <  1:619878/86‑1 (MQ=255)
                                                                                                                                                                                                           aTCGGTGTGGGGGTTAACCACCGGCGGTCCCGCCGTTTGCATTACTGTGCCGGTAAATCCCACCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATc                 <  1:295345/228‑1 (MQ=255)
                                                                                                                                                                                                                     agGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGTCAAAggtg       >  1:156481/1‑228 (MQ=255)
                                                                                                                                                                                                                       gTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAggtggc     >  2:112881/1‑227 (MQ=255)
                                                                                                                                                                                                                          aaaCACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCCGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGCAGTAGAGCTCGCCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCCTCAAAAAAGGGGGtttt  >  1:15246/1‑228 (MQ=255)
                                                                                                                                                                                                                          aaaCACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGtttt  >  2:44510/1‑228 (MQ=255)
                                                                                                                                                                                                                          aaaCACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGATGGtttt  >  2:248207/1‑228 (MQ=255)
                                                                                                                                                                                                                          aaaCACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGAGGTCTGGTCCAACCCATCGATAAAGTGGGttt   >  2:591596/1‑227 (MQ=255)
                                                                                                                                                                                                                          aaaCACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATAAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTTGTGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGAGTGCTGGTCCACACCATCGATAAAGGTGGtttt  >  1:562160/1‑228 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                               
CGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGC‑AGATCTCAGC‑GATAGCGCCGGCTTAGTCAGATTTAATC‑TGCGCGCGTGGTGG‑ATATTTTTTCAGGATCTCCA‑TATACGCGTGCATTTCGGTCTG‑T‑AGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTT  >  NC_000913/3984281‑3984719

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