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

CGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTT  >  NC_000913/3984281‑3984719
                                                                                                                                                                                                                              |                                                                                                                                                                                                                               
gttttggtggggtgggaaaGGTGCGTCCGCATGGCATCCGGCAAGGGCAGTTCCCAGTGATGGCTCCGGTTGAGTCAGTTTAAATCTGGGGGGGTGGTGGATATTTTTCCAGGATCGCCATATCCGCGTGCATTTCGGTCTGTAGCGGTACCCCCATCGGAATAGGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCccg                                                                                                                                                                                                                     <  1:434809/211‑1 (MQ=255)
  ttCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCccg                                                                                                                                                                                                                     >  2:588138/1‑226 (MQ=255)
  ttCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGGGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATcccc                                                                                                                                                                                                                     >  2:370066/1‑225 (MQ=255)
  ttCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGCGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGCGTCGCTTTCCTGCGGACCGGTGTAGAGGTTAAACACCGACGAACccgcc                                                                                                                                                                                                                   >  1:313587/1‑228 (MQ=255)
       tagcgtggaaaaggagggcaaGCATCGCATCGGGAAAGGCGAGTGCCGAGCGATTGCGCCGGGTTAGTAAAATTTATTCGGCGCGCGGGGGGGTTATTTTTTCAGTACCCCCATATCCGCGGGCATTTCGGTCTGTAGCGGTACCCCCACCGGAATAGGGGGCACGCCGATGGAGTCGCTTTCCTGAGGATCGGTGTAGAGGTTAAACACCGACGATACCGCCGTTTg                                                                                                                                                                                                              <  1:370066/208‑1 (MQ=255)
            tGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGCCGCACGCCGAGGGAGTGGCTTTCCTGCGGCTCGGGGTAAAGGTGAAAACCCGACGAGCCCGCCGtttacctac                                                                                                                                                                                                         >  1:225682/1‑222 (MQ=255)
            tGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTCTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCCCACGCCGATGGAGTCGCTTTCCTGCGGTACGGCGTCAAGGGTAAAACCCGAGAAGCCCGCCGTTTGCAttc                                                                                                                                                                                                         >  1:589503/1‑227 (MQ=255)
               aTAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGGGTAGGGGTTAAAACCCGACGGACCCGCCCTTTTCATAACt                                                                                                                                                                                                       >  1:338954/1‑227 (MQ=255)
               aTAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGGGTAGAGGTTAAAACCCGCCGATCCCGCCGTTTGCATTCCtg                                                                                                                                                                                                      >  2:316955/1‑228 (MQ=255)
                    aTGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAAAACCGCCGATCCCGCCGTTTGCATTACTGTGcc                                                                                                                                                                                                  >  2:573965/1‑227 (MQ=255)
                        gtaagctccgactcgggaaaaGGAAGTTCTCGGGGATGGCGCTGGCTTAGCCAGTTTTAATCTGGGTGGGTGGTGGATGTTTTTGCAGGATCTCCATATCGGCGTGCATTTCGGTCTTGAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGa                                                                                                                                                                                             <  1:588138/210‑1 (MQ=255)
                                cgcttcCGGCAAAGGGATCTCGCAGGGATGGCGCCGGTTTTTTCAGATTAATCGTGCGCGGGTGGTGGATTTTTTTTCAGCGTCTCATTAGACGGGTGCGTTTAGTGTTGAGGCGTTACCCCCAGTGGAATAGGGCGCACGCCGAGGGGGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCAcc                                                                                                                                                                                      <  1:38036/227‑1 (MQ=255)
                                  cactcGGGAAAGGGGGATCTCCGGGGATGGCGCCGGTTTGGTCAGTTTTAATCTGGGGGGGGGGTGGATATTTTTTCAGGATCCCAAAATCCGCGTGCATTTGGGTCTGGGGCGGTACCCCCATCGGAATATGGCCCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAACCCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCAcc                                                                                                                                                                                      <  1:573965/222‑1 (MQ=255)
                                                      tccggggttggggcCGGTTTGGTCAGATTTAACCTGGGGGCGTGGTGGTTTTTTTTTAGGGATCCCCATTTCCGGGTGTCTTTGGGTTGTGAGCGGTACCCCCACGGGAATATGGCGCACGCCGATGGAGTCGCTTTCCGGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGtat                                                                                                                                                                  <  1:330256/215‑1 (MQ=255)
                                                         gCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTTCATTACTGTGCCGGTGAATCCACCCGGATATCCCCCCTGGGCAAAAGCg                                                                                                                                                             >  1:651733/1‑228 (MQ=255)
                                                             ttgggtccGGCTTAGTCAGATTTAATCTGGCAGCGTGGTGGATATTTTTTCAGGATCCCCATATACGCGTGAATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTa                                                                                                                                                           <  2:517224/221‑1 (MQ=255)
                                                               aGCCCCGGTTTAGTCAGTTTTAATCGGCGCGGGTGGTGGATATTTTTTCGGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAg                                                                                                                                                         <  2:338954/227‑1 (MQ=255)
                                                                gccgcgttttgggcataTTAAATCTGCGGGCGTGGTGGATTTTTTTACAGGTTCCCCATATACCGGTGCTTGGTGTTCTGGAGGGGCACCCCCCTCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTaa                                                                                                                                                          <  1:532575/209‑1 (MQ=255)
                                                                 cgcCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGGGAATCCACCCTGATTTCCGCCCTGGGCTTAAGCGTAAGGt                                                                                                                                                       >  2:451300/1‑227 (MQ=255)
                                                                  tcgcgtttAGTCAAATTTAATCCGCGCGCGTGGTGGATATTTTTTGAGGTTCTCCATATCGGGGTGCATTGCGGTCTTTAGCGGTCCACCCAGCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTg                                                                                                                                                     <  2:188857/222‑1 (MQ=255)
                                                                  gcCGGTTTAGTCAGTTTTATTCGGGCGGCTTGGTGGATATTTTTTCGGGTCTCCATTGTGCGGGAGCTTTTGGGTCTGTAGCGGTACACCCATCGGAATGTGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTg                                                                                                                                                     <  1:354588/228‑1 (MQ=255)
                                                                     taaaggacaaGATTTTATCTGGGGGGGGTGTGGAATATTTTTCCAGGTTCCCCATATGCGCGGGCTTTTCGTGTGTTAGGGTTACCGCCTTGGAAATGTGGCGCACCCCGATGGGGTCGCTTTCCGGCGGATCGGTGTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCt                                                                                                                                                   <  1:116357/219‑1 (MQ=255)
                                                                                aTTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCCCTGGGGATAAGCGTAAGGTGGCTAAATCAGGACgg                                                                                                                                       >  1:316447/1‑227 (MQ=255)
                                                                                     aTCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGGTAATCCACCCTGCTATCCGCTCTGGGAATAAGGGTCAGGTTGCTGAATCAGGCCGTGaac                                                                                                                                   >  2:295345/1‑225 (MQ=255)
                                                                                           cgcgcgTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTCCTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGaa                                                                                                                             >  2:619751/1‑228 (MQ=255)
                                                                                           cgcgcgTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGc                                                                                                                                                    >  2:116357/1‑205 (MQ=255)
                                                                                           cgcgcgTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGAAATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGa                                                                                                                              >  2:169562/1‑227 (MQ=255)
                                                                                           cgcgcgTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGGGAATCCACCCTGATATCCGCTCTTGGTATAAGCGTAAGGTTGCTGAATACGGACGTGATACTTGaa                                                                                                                             >  2:247010/1‑228 (MQ=255)
                                                                                           cgcgcgTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGAGTAGAGGGTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGAGAATCCACCCTGATATCCGCTCTGGGTATACGCGTAAGGTTGCTGAATCAGAACGGGCAACTTGa                                                                                                                              >  2:399315/1‑227 (MQ=255)
                                                                                               cgTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGAATAAGCGTAAGGTTGCTGAATCAGGACGTGATCCTTGAACTct                                                                                                                         >  2:576399/1‑227 (MQ=255)
                                                                                                          atTTTTTCAGGATCCCCATATGCGCTTGCATTTCGGTCTGTAGCGGTACACCCATGGGAATATGGCGCAGGCCGATGGAGTCGCTTTCCGGGGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGc                                                                                                                <  2:492098/227‑1 (MQ=255)
                                                                                                                   aggatccccattagccggTGCATTTCGGTCTGGAGCTGTACCCCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCACCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACagcagc                                                                                                        <  2:612982/210‑1 (MQ=255)
                                                                                                                         cccatgtccggggtcATTTCGGTCTGTACGGGTCCCCCCAGCGGAATAGGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGGTTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTa                                                                                                 <  2:315041/214‑1 (MQ=255)
                                                                                                                             tataCGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGGATCAGGACGTGATACTTGAACTCATCCAAACGCACCGAAGCGAGTTTCCCGt                                                                                             >  1:511317/1‑228 (MQ=255)
                                                                                                                                   cgTGCATTTCGGTCTGTAGCGGTACACCCCTCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCGGATATCCGCTCTGGGTCTCAGCGGAAGGGTTCTGAATCAGGACATGATACTTGAACTCATACATACGCCACGCGGCAAGGTTACCGGTGAGGa                                                                                       >  1:148073/1‑228 (MQ=255)
                                                                                                                                    tgtctctcttggggttgtgggggTACCCCAATGGGAATGTGGGGCACGCCGGTGGGGTCCGTTGCCTGCGGTGCGGGGTAGAGTTTAACCCCCGCCGTTCCCGCCGTTTGCAGTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAgg                                                                                        <  1:545055/207‑1 (MQ=255)
                                                                                                                                     tGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGaa                                                                                                                                                >  2:20539/1‑169 (MQ=255)
                                                                                                                                     tGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGCCGTGATACTTGAACTCATCCATACGAACAGCCGCGAGTTTACCGTTGAGGaaa                                                                                     >  1:315041/1‑227 (MQ=255)
                                                                                                                                     tGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATCAGCACAGCAGCGAGTTTACCGTTGAGGaa                                                                                      >  2:326664/1‑227 (MQ=255)
                                                                                                                                     tGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAACCGTAAGGTTGCTGAACCAGGACGTGATACTTGAACTCATCCATACGCCAAGCAGGAAGTTTACAGTTGAGGaa                                                                                      >  1:334773/1‑227 (MQ=255)
                                                                                                                                     tGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCCCTCTGGGTATAAGCGTAAGGTTGCTGAATTAGGACGGGATACTTGAACTCATCCATCCGCACAGCACCGCGTTTACCGTTGAGGag                                                                                      >  1:646201/1‑226 (MQ=255)
                                                                                                                                       cATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTATAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGACTCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTACCCGTTGGGGAagta                                                                                   >  2:182698/1‑228 (MQ=255)
                                                                                                                                            cggtcggtAGCGGTACCCCCATGGGAATATGGCGCACGCCGATGGAGTCGCTTTCCGGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCt                                                                              <  2:562756/228‑1 (MQ=255)
                                                                                                                                                    tAGGGGTACACCAATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCt                                                                        <  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)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                               
CGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTT  >  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: