Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 403,205 (T)7→8 intergenic (+43/‑76) phoA → / → psiF alkaline phosphatase/phosphate starvation‑inducible protein PsiF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913403,2051.T88.6% 107.0 / 4.2 35intergenic (+43/‑76)phoA/psiFalkaline phosphatase/phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (15/16);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑T‑GCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGC  >  NC_000913/402958‑403451
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
aGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCGGGGTGG‑ttttt‑t‑ttgc                                                                                                                                                                                                                                                      <  1:353517/251‑4 (MQ=255)
            cACCCAGGCGCTAAAGACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAATAGGATTCACAAGAACATACCGGCAGTCATTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCa                                                                                                                                                                                                                                          <  2:178444/251‑1 (MQ=255)
                 aGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACtt                                                                                                                                                                                                                                     >  1:132493/1‑251 (MQ=255)
                         aaTACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGc                                                                                                                                                                                                                               >  2:33975/1‑249 (MQ=255)
                                                 gtgatgAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGt                                                                                                                                                                                                     <  1:115172/251‑1 (MQ=255)
                                                           aCGGGCACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCCCGGGCGCATACGCTCATGGGTAAAACATg                                                                                                                                                                                            >  2:81542/1‑250 (MQ=255)
                                                             gggAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATCCGCTCATGGTTAAAACCTGa                                                                                                                                                                                           >  2:25271/1‑249 (MQ=255)
                                                               gAACTCAGCAGAGGATTCACAAGCACATACAAGCACTCAGTTGCGTCTTGCACCGTATGGACCGCATCCCCCCATTGTTGTTGGACTGCCAGACAAGACAGACCTCTTCTCAACCCTGCAAGCCGCTCTGGGGCTGAAATCACACCGCGCCCCGCATGGCATTTTCGCTGCCGGGTGG‑TTTTT‑TCT‑GCTGTTAGCAACCAGACTTAATGGCAGATCACCGGTGCATACGCTCATGGTTAAAACATGAag                                                                                                                                                                                         >  2:267734/1‑249 (MQ=255)
                                                                    ccGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGAt                                                                                                                                                                                   <  1:343754/250‑1 (MQ=255)
                                                                           ggATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGGTGAAATAAAACCTCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGGTAGCAACCAGACTTAATGGCAGATCACGGGCGGATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                                                                                                                                                           <  1:251350/251‑1 (MQ=255)
                                                                           ggATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                                                                                                                                                           <  1:171880/251‑1 (MQ=255)
                                                                           ggATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGGTAAAACATGAATAGGGTTGGTGCTAt                                                                                                                                                                           >  2:10090/1‑251 (MQ=255)
                                                                                      aaCATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGCTGGTGCTATGAAAATAACAt                                                                                                                                                                >  2:46739/1‑251 (MQ=255)
                                                                                            ccGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTgg                                                                                                                                                          <  1:33975/251‑1 (MQ=255)
                                                                                                    cagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTg                                                                                                                                                  >  1:306770/1‑251 (MQ=255)
                                                                                                       ttGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCtt                                                                                                                                               >  1:185823/1‑251 (MQ=255)
                                                                                                                    cgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttt                                                                                                                                  <  2:88946/251‑1 (MQ=255)
                                                                                                                           cccGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAac                                                                                                                           <  1:114499/251‑1 (MQ=255)
                                                                                                                                  gcctccaatggtggtagccTGACCGACCAGACCGAGCCCGTCTACACCATGAAATCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGATTTTT‑T‑T‑GCGGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATTAAGAGGGATGGTGCTATGAAAATAACATGACTGGTTACCTGGCTTTTCGTTCTGGTTTCTTTAACCACCGTc                                                                                                                    <  2:263956/233‑1 (MQ=255)
                                                                                                                                  gccgccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGt                                                                                                                     <  1:75794/250‑1 (MQ=255)
                                                                                                                                     gccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATGACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGgc                                                                                                                 <  1:198953/251‑1 (MQ=255)
                                                                                                                                     gccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACAATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGgc                                                                                                                 >  1:145252/1‑251 (MQ=255)
                                                                                                                                     gccAATGTTGTTGGACTGACCGACAAGACCGATCTCTTCTACACCGGGAAAGCCGCTCTGGGGCGGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCGATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGgc                                                                                                                 <  2:148779/251‑1 (MQ=255)
                                                                                                                                           gttgttGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGcc                                                                                                           >  1:82370/1‑251 (MQ=255)
                                                                                                                                                   aCTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTGACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGGCGGCGCTGCCGAGAGAAc                                                                                                   <  2:404414/251‑1 (MQ=255)
                                                                                                                                                                                ccATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTACCCCCACAACACCAGCGTCTGAcct                                                                      >  2:145361/1‑251 (MQ=255)
                                                                                                                                                                                 cATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcctc                                                                     >  2:203251/1‑251 (MQ=255)
                                                                                                                                                                                 cATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcctc                                                                     <  1:73381/251‑1 (MQ=255)
                                                                                                                                                                                              ctGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGGCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATcagc                                                          <  2:340303/249‑1 (MQ=255)
                                                                                                                                                                                                ggggCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTTGT‑T‑GCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCGCCTGTAATcagcag                                                        <  2:161075/249‑1 (MQ=255)
                                                                                                                                                                                                  ggCTGAAATAATACCGCGCGCGGCAACGACTATTCGCTGGCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCCCGGGAGCATACGCTCCTGGTTAAAAAATGAAGAGGGATGGAGCGATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAACAAGCAGGCg                                                     >  2:77318/1‑250 (MQ=255)
                                                                                                                                                                                                       aaaTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATGACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACATCAGCGTATGACCTCCTGTAATCAGCAGGCGACGgcg                                               <  1:249484/251‑1 (MQ=255)
                                                                                                                                                                                                             aaCCGCGCCCGGCAGTGAATTTTCGATGCCGGGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGCCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGgc                                                                                                                 >  2:220078/1‑179 (MQ=255)
                                                                                                                                                                                                               ccGCCCCCGGCAGTGAATTTTCGCTGCCGTGTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTGAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTGTTCGGTCTGGTTTTTTTAACCACCGGCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGtt                                       <  1:395124/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggcggttgag‑tTTTT‑T‑G‑GCGGTTATGAACCAGACTTAATGGCAGATCCCGGGCGCAGCCGCTCATGGTTAAAACACGAAGAGGGATGGTGCTATTAAAATAACATTACTGGTTACCTTGCTTTTCGGTCGGGTTTTTTTACCCACCGTCGGCGCTGCCGAGAAAAATTTAACCCCACAACAACAGCGAATGACCACCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCt                <  2:139433/240‑1 (MQ=255)
                                                                                                                                                                                                                                           gggTGG‑TTTTT‑T‑TTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATg           >  2:332210/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gtggtttttt‑t‑t‑GCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATg           >  1:22707/6‑249 (MQ=255)
                                                                                                                                                                                                                                               ggtttttt‑t‑t‑GCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTg       <  1:326315/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttt‑t‑t‑GCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTAGGACCTCCTGTAATCAGCAGTCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTg   <  2:274694/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     tt‑t‑t‑GCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGc  <  2:338994/250‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
AGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGG‑TTTTT‑T‑T‑GCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGC  >  NC_000913/402958‑403451

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