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.T95.9% 181.6 / 0.5 49intergenic (+43/‑76)phoA/psiFalkaline phosphatase/phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base T (29/18);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCT  >  NC_000913/402959‑403453
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
gCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATGGCGGCGTATGGCCCGCATGCCGCGAATGTTGTTAGACTGACCGACCAGACCGAGCTCTTCTACACCATGAAATCCGCTCTGGTGCTGAAATAAAACCTCCCCCGGCAGTGAATTTTCGCTGCCGGGTGGttttttttgct                                                                                                                                                                                                                                                       <  1:605056/251‑5 (MQ=255)
             cccAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAAcc                                                                                                                                                                                                                                             >  1:673777/1‑248 (MQ=255)
              ccAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGAc                                                                                                                                                                                                                                         >  1:2394/1‑251 (MQ=255)
              ccAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGAc                                                                                                                                                                                                                                         >  1:3223/1‑251 (MQ=255)
                   cgcTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAAt                                                                                                                                                                                                                                    >  1:480090/1‑251 (MQ=255)
                       aaaTACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGc                                                                                                                                                                                                                                 >  1:16979/1‑250 (MQ=255)
                       aaaTACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGc                                                                                                                                                                                                                                 >  1:693248/1‑250 (MQ=255)
                                                gtgatgAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGt                                                                                                                                                                                                       <  1:382736/251‑1 (MQ=255)
                                                            gggAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGa                                                                                                                                                                                             >  1:199772/1‑249 (MQ=255)
                                                                   ccGAAGAGGATTCACAAGAACATACCGGCAGGCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATg                                                                                                                                                                                    <  2:477088/251‑1 (MQ=255)
                                                                      aaGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTg                                                                                                                                                                                 <  1:37544/251‑1 (MQ=255)
                                                                                 caAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCTCGCTTGCCGCCAATGTTGTTGGACTGACCGACCCGACCGATCTCTTCCACACCATGACAGCCGCTCTGGGGCTGACATAAAACCGCGCCAGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCTGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACCTGAAGAGGGCTGGTGCTATGAAAAt                                                                                                                                                                       >  1:892492/1‑250 (MQ=255)
                                                                                      aCATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGgcgc                                                                                                                                                                                                                    >  1:212260/1‑200 (MQ=255)
                                                                                      aCATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGgcgc                                                                                                                                                                                                                    <  2:212260/200‑1 (MQ=255)
                                                                                        aTACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTAc                                                                                                                                                               >  1:365098/1‑251 (MQ=255)
                                                                                           ccGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCCGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCCTACGCTCATGGTTTAAACATGAAGAGGGATGGTGCTATGAAATTAACATTACTgg                                                                                                                                                            >  2:621892/1‑251 (MQ=255)
                                                                                                aggcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAAGGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGGGAATTTTCGCGGCCGGGTGGTTGTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAATAGGGATGGTGGTATGAAAATAACATTACTGGTTAcc                                                                                                                                                       <  2:569050/248‑1 (MQ=255)
                                                                                                       tGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCttt                                                                                                                                                >  1:73379/1‑251 (MQ=255)
                                                                                                            aTTGCGGCGTATGGCCCGCATGCCGACAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGt                                                                                                                                           >  1:785254/1‑251 (MQ=255)
                                                                                                                 ggcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGt                                                                                                                                                                                  <  2:14878/207‑1 (MQ=255)
                                                                                                                 ggcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGt                                                                                                                                                                                  >  1:14878/1‑207 (MQ=255)
                                                                                                                  gcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGtt                                                                                                                                     >  1:327920/1‑251 (MQ=255)
                                                                                                                           ccGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGTGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTGTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACTTTCATGGTTAAAACATGAAGAGGGATGGTGCTATGTAAAGAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAacc                                                                                                                            <  1:585642/251‑1 (MQ=255)
                                                                                                                               aTGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCg                                                                                                                        >  1:252787/1‑251 (MQ=255)
                                                                                                                                     ccAATGTTGTTGGACTGCCCGACCAGACCGATCTCTTCTACACCATGAAAGCAGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATCACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGgcg                                                                                                                  >  2:285295/1‑251 (MQ=255)
                                                                                                                                                  aCTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCTCCCTGCATTTAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACTTGAATAGGTATTGTGTTATGAAAATAACATTACTGGTTACCTTGGTTTTCTGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAAc                                                                                                     <  1:96707/251‑1 (MQ=255)
                                                                                                                                                     gaccgaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACt                                                                                                    <  2:480090/249‑1 (MQ=255)
                                                                                                                                                                  gATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCacaaca                                                                                      >  2:630384/1‑250 (MQ=255)
                                                                                                                                                                   aTCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCacaacaac                                                                                    >  2:608883/1‑251 (MQ=255)
                                                                                                                                                                    tctcTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCacaacaac                                                                                    >  1:270336/1‑250 (MQ=255)
                                                                                                                                                                      tcttctACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCacaacaaccgc                                                                                 >  2:874360/1‑248 (MQ=255)
                                                                                                                                                                       cttcaccacCAGGAAAGCCGCCCTGGGGCTAAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCacaacaac                                                                                    <  2:2394/241‑1 (MQ=255)
                                                                                                                                                                            acgccATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAAGTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTGAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCTTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTAt                                                                             <  1:621892/246‑1 (MQ=255)
                                                                                                                                                                                cATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcctc                                                                       <  2:3223/251‑1 (MQ=255)
                                                                                                                                                                                   gaaagcctgtctgtgtgggAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCAAACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCGTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCTGCGCGGCGGAGAGAACGTTAACCCCACGTCAGCAGCGTATGACCTCCTg                                                                    <  2:892492/233‑1 (MQ=255)
                                                                                                                                                                                       gCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGCGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAAt                                                                >  1:570319/1‑251 (MQ=255)
                                                                                                                                                                                             ctGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATcagc                                                            <  2:673777/249‑1 (MQ=255)
                                                                                                                                                                                             ctGGGGATGAAATAAAACCGCGCACGGCAGAGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAACCATGAAGCGGGATGGTGCTATGAAAATAACATTACTGGTTACATTGCTTTTCGGTCTGGTTTTTTTCACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATcagc                                                            >  2:375886/1‑249 (MQ=255)
                                                                                                                                                                                                          aaaaCCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGaa                                                                                                      <  1:491209/194‑1 (MQ=255)
                                                                                                                                                                                                          aaaaCCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGaa                                                                                                      >  2:491209/1‑194 (MQ=255)
                                                                                                                                                                                                                   cccGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGAc                                                     >  2:148453/1‑234 (MQ=255)
                                                                                                                                                                                                                   cccGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGAc                                                     <  1:148453/234‑1 (MQ=255)
                                                                                                                                                                                                                    ccGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTAAGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAgg                                   <  1:608883/251‑1 (MQ=255)
                                                                                                                                                                                                                      ggCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCACCCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGAt                                                                                                                                                                                     >  1:587934/1‑103 (MQ=255)
                                                                                                                                                                                                                      ggCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCATATCACGGGCGCATACGCTCATTGTTAAAACATGAAGAGGGAt                                                                                                                                                                                     <  2:587934/103‑1 (MQ=255)
                                                                                                                                                                                                                            gAATTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTc                           <  2:693248/251‑1 (MQ=255)
                                                                                                                                                                                                                             aaTTTTCGCTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCg                          <  2:252787/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGCCGGGTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCt                  >  1:915020/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gggTGGTTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGCAAGGGGATGCTCGTAAGACCTACATg             >  1:589256/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ggtggtTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATg             >  1:217935/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    ttttGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGcc   >  1:806584/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      ttGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCt  <  2:785254/250‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
GCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCT  >  NC_000913/402959‑403453

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: