Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,911,819 T→G 8.4% intergenic (‑294/+20) pstS ← / ← glmS phosphate ABC transporter periplasmic binding protein/L‑glutamine:D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,911,8190TG 8.4% 147.9 / 5.7 51intergenic (‑294/+20)pstS/glmSphosphate ABC transporter periplasmic binding protein/L‑glutamine:D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base G (3/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

ATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATA  >  NC_000913/3911524‑3912118
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                           
ccgatctttcCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGc                                                                                                                                                                                                                                                                                                        <  1:504246/293‑1 (MQ=255)
aTAATGTGTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGc                                                                                                                                                                                                                                                                                                        <  1:273775/301‑1 (MQ=255)
aTAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGACACGATTACGc                                                                                                                                                                                                                                                                                                        >  2:466940/1‑301 (MQ=255)
   aTGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAgg                                                                                                                                                                                                                                                                                                     <  1:556053/301‑1 (MQ=255)
     gTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACACAggga                                                                                                                                                                                                                                                                                                   >  2:459163/1‑300 (MQ=255)
      tctcCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGCGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCACCGCTTACGCAgggac                                                                                                                                                                                                                                                                                                  >  1:239295/1‑299 (MQ=255)
       ctcCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAgatcgg                                                                                                                                                                                                                                                                                                 >  2:504918/1‑296 (MQ=255)
            tggatgaTTTATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCTGCGTAAGCGGGGCAtttttc                                                                                                                                                                                                                                                                                            <  2:277789/300‑6 (MQ=255)
              gaggaTTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAAAAGGAAGAAAAATGCACCGATTACGCAGGGCATCCAttat                                                                                                                                                                                                                                                                                          >  1:228072/1‑299 (MQ=255)
                        aaaGCATTGTTTGTTGGCGACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAAcc                                                                                                                                                                                                                                                                                <  1:57210/301‑1 (MQ=255)
                             tttgtttgtTTGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAAAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAAc                                                                                                                                                                                                                                                                           <  2:324084/300‑1 (MQ=255)
                             tttgtttgtTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAAc                                                                                                                                                                                                                                                                           <  2:228315/300‑1 (MQ=255)
                                        gCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCa                                                                                                                                                                                                                                                                <  2:363118/301‑1 (MQ=255)
                                                 aGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAAcctggcaaaatcggttacggt                                                                                                                                                                                                                                                       <  2:9231/301‑20 (MQ=255)
                                                        ccctacacgacgctcttccgatcttGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTc                                                                                                                                                                                                                                                <  1:155473/277‑1 (MQ=255)
                                                             cAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACAGAtcggaagagcacacgtctgaactccagta                                                                                                                                                                                                                                           >  2:633780/1‑273 (MQ=255)
                                                                                tGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCagatcggaagagcacacgtctgaa                                                                                                                                                                                                                        >  2:155624/1‑278 (MQ=255)
                                                                                  tAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAgcgcac                                                                                                                                                                                                                      >  1:422599/1‑299 (MQ=255)
                                                                                     ggAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGAGACAt                                                                                                                                                                                                                   >  1:503574/1‑301 (MQ=255)
                                                                                                                  aaaGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGgtg                                                                                                                                                                                      <  2:423162/301‑1 (MQ=255)
                                                                                                                                taTTTAAATAATCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGc                                                                                                                                                                        <  2:283192/300‑1 (MQ=255)
                                                                                                                                    tAAATAAGCGTTGAGATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCAGTCTGTGACAGGGAAAAAGTAGCGGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACATGAAGAAAAATGCCCCGCGTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAGCGGTGTAGAAGATCGGTGCAAtg                                                                                                                                                                    <  1:458599/301‑2 (MQ=255)
                                                                                                                                              ttGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCca                                                                                                                                                          <  2:504246/301‑1 (MQ=255)
                                                                                                                                                                      ttAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGCATTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTACCATCACCTCTTACACATGCGGAATATCgatgacgacaat                                                                                                                                 >  1:523599/1‑294 (MQ=255)
                                                                                                                                                                         ataaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGAAGATGTGAATGt                                                                                                                               >  1:277113/1‑301 (MQ=255)
                                                                                                                                                                           aaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGtaa                                                                                                                             >  1:149597/1‑299 (MQ=255)
                                                                                                                                                                            aCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTAt                                                                                                                            >  2:560302/1‑301 (MQ=255)
                                                                                                                                                                            aCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGACTCCCGATGCTGTCAATGTTAt                                                                                                                            >  1:524641/1‑301 (MQ=255)
                                                                                                                                                                              gAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATGTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCg                                                                                                                          <  2:239562/301‑1 (MQ=255)
                                                                                                                                                                               aaGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGc                                                                                                                         <  1:440760/301‑1 (MQ=255)
                                                                                                                                                                                  agaTGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTAc                                                                                                                      <  2:123557/301‑1 (MQ=255)
                                                                                                                                                                                                    cATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCt                                                                                                    >  2:71137/1‑301 (MQ=255)
                                                                                                                                                                                                                        gTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACaaa                                                                                >  2:43818/1‑300 (MQ=255)
                                                                                                                                                                                                                         tAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACt                                                                               <  2:525313/301‑1 (MQ=255)
                                                                                                                                                                                                                              gAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCAACCTGATCGGCAAAAACATACAACTAAccc                                                                          >  1:188637/1‑300 (MQ=255)
                                                                                                                                                                                                                               aaGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGc                                                                         <  2:277416/301‑1 (MQ=255)
                                                                                                                                                                                                                                     aCGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGACCGGCAAAAACATCCAACGCACCCCCAcgcc                                                                   >  1:453216/1‑300 (MQ=255)
                                                                                                                                                                                                                                              cacaTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAActt                                                          <  2:188879/301‑1 (MQ=255)
                                                                                                                                                                                                                                                caTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGCGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAAACATACAACTGACAGCCACGCGCCCGAActtct                                                        >  1:168121/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         ggCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACACGCATCCTGATCGGCGAAGACATACAACTAACCGCCACGCGCCCGAACTTCTCCAATGTTg                                               >  1:481780/1‑301 (MQ=255)
                                                                                                                                                                                                                                                            atgATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCGTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGAt                                            <  2:368869/299‑1 (MQ=255)
                                                                                                                                                                                                                                                               atggtggATGAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCTTCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTc                                         <  1:466376/295‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTcc                                >  2:190797/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                         aaaaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGCTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGAACGGCGAAGACATACAACTGACCTACACTCGCGCGAACTTCTTCAATGGTGGATTTCAGTTTTTccc                               >  1:626836/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                          aaaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCCTCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATCCAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTcaac                              >  1:82888/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                    aGGAAGAAAAATGCCCCGCTTACGCAGGGCATTCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgttg                    <  1:190554/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    aGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgttg                    <  2:168318/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         aatggaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTAGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATACGTTTTTCgg               >  1:521692/6‑301 (MQ=255)
                                                                                                                                                                                                                                                                                            ggaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACTGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTAGTTGTTCGGTGc            >  2:456521/3‑301 (MQ=255)
                                                                                                                                                                                                                                                                                            ggaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGc            >  2:248577/3‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                      cTTACTCAGGGCATCCATTTATTACTCAACCGTATCCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgttgttgcaccgaacaacgaat  <  1:71047/301‑17 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                           
ATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATA  >  NC_000913/3911524‑3912118

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