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.T97.7% 160.7 / ‑2.8 44intergenic (+43/‑76)phoA/psiFalkaline phosphatase/phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (26/17);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TT‑GCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATT  >  NC_000913/402960‑403449
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
cgccGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGttttt‑tttgctg                                                                                                                                                                                                                                                  <  1:44628/249‑6 (MQ=255)
  ccGGGCCGCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGGTTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGttttt‑tttgctg                                                                                                                                                                                                                                                  <  2:172511/249‑6 (MQ=255)
          cACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCa                                                                                                                                                                                                                                        <  1:11305/251‑1 (MQ=255)
             ccAGGCGCTAAATCCCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGAc                                                                                                                                                                                                                                     >  2:428766/1‑251 (MQ=255)
             ccAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATTGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGAc                                                                                                                                                                                                                                     <  1:576925/251‑1 (MQ=255)
               aGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGAGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTATGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACtt                                                                                                                                                                                                                                   >  1:530739/1‑251 (MQ=255)
                          aCCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGc                                                                                                                                                                                                                             >  1:51224/1‑246 (MQ=255)
                          aCCAAAGATGATGCAGTGATCGTGATGACTCACGGGAACTACGAAGAGGACACACAAGAACATACCGGCAATTAGTTGCCAATCGCCGCGTACGACCCCAATGTCGCCAATCTTGTTGGACTGACCGACCAGACCGATCTATTCTACACCATGAAAGCCGCTCTGGGGCTGAAATCAAACCGCGACCGGCAGTGAATTTTTGCTGCAGGGTGGTTTTT‑TTTGCTGTTAGAACCCTGACTTAATGGCAGAtt                                                                                                                                                                                                                        >  2:383492/1‑250 (MQ=255)
                                    gcgcAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCAt                                                                                                                                                                                                              <  2:108867/251‑1 (MQ=255)
                                        aGTGATGGTGATCAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGc                                                                                                                                                                                                          >  2:205119/1‑251 (MQ=255)
                                            atggtgatgAGTTACGGGAACTCCGACGAGGATTCTAAAGAACATACCGGCAGTCAGTTGCGTATTCCGGCGTATGGCACGCCTGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCGCTGGGGATGAAATAACACCGCGCCCGGCCGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTATTGGCAGATCACGGGCGCCTTCGCTCAt                                                                                                                                                                                                      >  2:383201/1‑251 (MQ=255)
                                                        tACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATg                                                                                                                                                                                          >  2:93913/1‑251 (MQ=255)
                                                                    gAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCAGCCAATGATGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGACATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGAAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATggg                                                                                                                                                                              >  1:553743/1‑250 (MQ=255)
                                                                       gagGATTCAGTAGAACATACCGGCAGGCAGTTGCTTCTTGCGACGTATGGCCCGTATGCCGCCAATGTTGTTGGGCTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGTATTTTCGCTTCCGGGGGGTTTTTGTT‑GCTGTTAGCAGCCAGACTTAAGGGCGGATCACGGGCGCATACGCTGATGATTAAAACATGAAGAGGGATGGTGCt                                                                                                                                                                           <  1:316704/251‑1 (MQ=255)
                                                                       gagGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCt                                                                                                                                                                           <  2:51224/251‑1 (MQ=255)
                                                                         ggATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTAGGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTAAACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCAGGCAGTGAATTTTCGCTGCCGGGGGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGGTCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                                                                                                                                                         >  2:102223/1‑251 (MQ=255)
                                                                          gATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATg                                                                                                                                                                        >  2:438751/1‑251 (MQ=255)
                                                                                caAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAAt                                                                                                                                                                   >  1:246152/1‑250 (MQ=255)
                                                                                 aaGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATa                                                                                                                                                                  >  1:478549/1‑250 (MQ=255)
                                                                                         aCCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTg                                                                                                                                                         >  1:305685/1‑251 (MQ=255)
                                                                                            ggCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGtt                                                                                                                                                      >  2:495064/1‑251 (MQ=255)
                                                                                                   agtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTCTGAAAATAACATTACTGGTTACCTTGc                                                                                                                                               >  2:459885/1‑251 (MQ=255)
                                                                                                                ggcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGt                                                                                                                                  >  1:332849/1‑251 (MQ=255)
                                                                                                                   gTATGGCCCGCATGCCTCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACGTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGtttt                                                                                                                               <  1:421240/251‑1 (MQ=255)
                                                                                                                    tATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttt                                                                                                                              >  1:322209/1‑251 (MQ=255)
                                                                                                                                gccgccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTc                                                                                                                  <  2:141512/251‑1 (MQ=255)
                                                                                                                                       aTGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTg                                                                                                           <  1:205119/251‑1 (MQ=255)
                                                                                                                                               ggACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGTCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTTCTTTTTGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCgagaga                                                                                                   <  1:351065/251‑1 (MQ=255)
                                                                                                                                                    gaccgaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACttt                                                                                              <  1:244042/251‑1 (MQ=255)
                                                                                                                                                      ccgaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTa                                                                                             >  1:127751/1‑250 (MQ=255)
                                                                                                                                                         accAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCAGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGACGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACAACCGTCGGCGCTGCCGAGAGAACTTTAAccc                                                                                         >  1:531041/1‑251 (MQ=255)
                                                                                                                                                              aCCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCcaca                                                                                     >  1:204501/1‑250 (MQ=255)
                                                                                                                                                                    ctctTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCacaacaac                                                                                >  1:609455/1‑249 (MQ=255)
                                                                                                                                                                      cttctACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATTAAAATAACATTACTGGTTACCTGGCTTTTCGGTCTGGTGTGTTTAACCTCCGTCGGCGCTGCCGAGATAACGTTAACCCCacaacaac                                                                                <  2:360320/247‑1 (MQ=255)
                                                                                                                                                                                            ctGGGGCTGAAATAAAACCGCGCCCGGCAGTGCATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCCCAACAACAGCGTATGCCCTCCTGTAATcagc                                                        >  2:269567/1‑249 (MQ=255)
                                                                                                                                                                                              ggggCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGc                                                    <  1:372853/251‑1 (MQ=255)
                                                                                                                                                                                                 gCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTCGCAACCAGACTTAATGGCAGCTCACGGGCGCATACGCTCCTGGGTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACACCCACAGCGTCTGCCCTCCTGTAATCAGCAGGCgcc                                                 >  2:263596/1‑249 (MQ=255)
                                                                                                                                                                                                    gAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGgc                                              <  2:428539/251‑1 (MQ=255)
                                                                                                                                                                                                       aTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCa                                           <  1:300705/251‑1 (MQ=255)
                                                                                                                                                                                                         aaaaCCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAgg                                         >  2:156486/1‑251 (MQ=255)
                                                                                                                                                                                                                 gcCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGaa                                  >  1:405838/1‑250 (MQ=255)
                                                                                                                                                                                                                      gCAGTGAAGTGTCGCTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCATCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGa                            <  1:383201/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCt              >  2:18530/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gCCGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATTGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTTCCGAGAGAACTTTAACCCCGCAACAACAGCGTATGACCTCCTGTAATCAGCAGTCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTa             <  1:267031/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACAt          <  1:18530/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGGTGGTTTTT‑TTTGCTGTTAGCAACCAGACGTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACAt          <  2:449454/251‑1 (MQ=255)
                                                                                                                                                                                                                                                tttt‑tt‑GCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGAtt  <  2:246152/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
CTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTT‑TT‑GCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATT  >  NC_000913/402960‑403449

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