Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,153,573	0	T	G	41.2%	105.1 / 28.0	68	Y92D (TAC→GAC)	pabC	4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base T (4/36);  new base G (28/0);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 3.61e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  NC_000913/1153429‑1153722                                                                                                                                                 |                                                                                                                                                        cAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAGA‑Gc                                                                                                                                                   >  2:853181/1‑150 (MQ=255)    cGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAA‑Gc‑a                                                                                                                                                 >  1:1303563/1‑148 (MQ=255)        tACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACA‑Gc‑a                                                                                                                                                 >  2:589273/1‑144 (MQ=255)         aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAAc                                                                                                                                         <  1:1315694/151‑1 (MQ=255)          cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACa                                                                                                                                        <  2:1441445/151‑1 (MQ=255)          cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACa                                                                                                                                        <  1:1636988/151‑1 (MQ=255)                gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGAc                                                                                                                                  <  2:1147597/151‑1 (MQ=255)                   tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCA‑GCAACAAGGAAAACCGGACCg                                                                                                                                >  2:1434514/1‑151 (MQ=255)                   tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGA‑GC‑ACAATGAACacacgggag                                                                                                                                >  2:306324/1‑142 (MQ=255)                       aGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAAc                                                                                                                           <  2:132499/151‑1 (MQ=255)                          ggTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAGA‑GC‑ACATTTAAAAACGGGCCGGCaa                                                                                                                            >  1:806943/1‑147 (MQ=255)                                 gatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATtct                                                                                                                  <  1:244874/150‑1 (MQ=255)                                  atTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATtctct                                                                                                                <  1:259692/151‑1 (MQ=255)                                      ccTGTGACTTCTGGCCTCACCTTGAACAAGAGGTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAGCAGCGGACCGGCAACGCGTATTCTCTCCGt                                                                                                            <  2:264228/151‑1 (MQ=255)                                                 tGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTTTCCGTTACGGCtt                                                                                                    <  2:267341/148‑1 (MQ=255)                                                   gCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCt                                                                                               <  2:702032/151‑1 (MQ=255)                                                     cTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACA‑AC‑ACATTTAAAAGCGGACCGGCCACGCGGGGTCTCTCCCTTTCGGCTTAt                                                                                                  >  1:1576536/1‑146 (MQ=255)                                                       cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGca                                                                                            <  1:521868/150‑1 (MQ=255)                                                        aGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGcaca                                                                                          <  2:1137394/151‑1 (MQ=255)                                                         gCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACA‑GC‑ACAATTAACAGCGGACCCGCGACACGGGTTTTCTCCGTTTCCGCCTATCCTGcc                                                                                            >  1:1310751/1‑147 (MQ=255)                                                          cTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACA‑GC‑ACAAAGAAAAGCGGGCCCGGCACGGGGGTTCTCTCCCTTCACGCTTATCCTGgacaat                                                                                        >  2:1574529/1‑149 (MQ=255)                                                            tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACAtt                                                                                        <  1:535241/149‑1 (MQ=255)                                                            tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACAtt                                                                                        <  1:409412/149‑1 (MQ=255)                                                             gAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACg                                                                                     <  1:1170001/151‑1 (MQ=255)                                                             gAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCA‑GG‑ACACTGGAAAGGGGGCCCGCAACCCGGGGTTTCTCCGTTAAGGCTTATCCCGccc                                                                                           >  2:529653/1‑143 (MQ=255)                                                                cAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGAGCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGa                                                                                    <  2:1409945/149‑1 (MQ=255)                                                                 aaGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCg                                                                                 <  1:1190266/151‑1 (MQ=255)                                                                 aaGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGGAGGA‑ACATTG‑ACAAAGGACCGGCGCCGCGGGGTGTCTCCCGTACGGCTTTTTCTGCCCACTACGAc                                                                                   >  2:886165/1‑149 (MQ=255)                                                                   gagaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCA‑GC‑ACATTGAACAACGGACCGGCAACGCGGGTTTTCTCCCTTTCGGCTTTTCCTGCACATTACGAcccgt                                                                               >  1:1562867/1‑148 (MQ=255)                                                                   gagaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGAA‑GC‑ACAATGAAAAACGGGCCCGCAACGCGGGGTCTCTCCCTTACGGCTTTTTCTGCACAATaaaacc                                                                                  >  1:1634170/1‑143 (MQ=255)                                                                    agaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACA‑GC‑ACAATTTACAACGGGCCGGCCACGCGGGGTCTCTCCCCTTCGGCCTATCCTGCACCTTACCACCGttt                                                                              >  1:405053/1‑151 (MQ=255)                                                                     gaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACA‑GC‑ACAATTAAAAGCGGGCCGGCCACGCGGGGTCTCTCCCTTACGGCTTATTCTGCACATTACGACCGttt                                                                              >  1:1623465/1‑150 (MQ=255)                                                                      aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGc                                                                            <  2:1195984/151‑1 (MQ=255)                                                                               cTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGCGGGGACA‑GC‑ACAATTAAAAGCGGACCGGCCAAGCGGGTTCTCTCCCTTTCGGCTTATTCCGCACATTACCACCGCTTTCGTAACGAgg                                                                   >  1:161246/1‑151 (MQ=255)                                                                                   cagcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCA‑GC‑ACACTTAAAAAAGGGACGGCAACGCGGAATCTTTCCCCTTCGGCCTATTCTTCACATTAAGACCCCTTGCGGAAAGAgggg                                                                 >  1:1285411/1‑149 (MQ=255)                                                                                       agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACg                                                           <  1:1353484/151‑1 (MQ=255)                                                                                       agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACg                                                           <  2:1467551/151‑1 (MQ=255)                                                                                       agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACg                                                           <  2:751331/151‑1 (MQ=255)                                                                                         aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACA‑GG‑AAAATTAACAACGGGCCGGGCACGCGGGTTCTCTCCCCTTCGGCCTTTCCTGCACATTTCAACCCTTTGCGTAACCAGGGGGttttg                                                           >  2:1582022/1‑146 (MQ=255)                                                                                         aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAA‑AC‑ACATTGAACAGCGGGACCGCAAACCGGGTTTTCTCCCGTACGGCGTTTTCTTCACATTTCGACCGGTTTGGTAACGAGGGGGGTACGtt                                                         >  1:1299158/1‑151 (MQ=255)                                                                                          aCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTg                                                        <  1:1188742/151‑1 (MQ=255)                                                                                           cAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCA‑GC‑ACATTTAACAGCGGACCGGCAACCCGGGTTCTCCCCCTTTCGGCTTATCCTGCACATTACAAACGGTTGCGTAACGAGGGGGTTACGtttg                                                       >  1:1512499/1‑149 (MQ=255)                                                                                               aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAA‑GC‑ACATTTAAAAACGGGCCCGCCACACGGGGTTTCTCCCTTTCGGCGTTTTCTGCACATTACCACCCGTTTTGGAAAGAGGGGGGTACGtt                                                         >  1:405558/1‑145 (MQ=255)                                                                                               aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCA‑AG‑AAATTGAACAACGGACCGGCCACGCGGGGTCTCTCCGGTTCGGCTTATTCTGCCCCTTACGACCGCTTTCGTAAAGAGGGGGTTTCGTTTgcg                                                     >  2:734671/1‑149 (MQ=255)                                                                                                aaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGCGGGGGCC‑GC‑ACATTTTACAGCGGACCGGCCACGCCGGGTCTCTCCCTTTCGGCTTTTCCTTCCCCTTACCACCCTTTTCGTTACGAGGGGGGTTCGTTTgcg                                                     >  2:26919/1‑148 (MQ=255)                                                                                                  aTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAg                                                <  1:970767/151‑1 (MQ=255)                                                                                                       gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCA‑GC‑ACATTGAAAAACGGACCGGCCACGCGGGGTCTCTCCGCTTCGGCCTTTCCTGCACATTAAGACCGGTTGCGTGAAGAGGGGGTTACGTTGGGGCTAAGccccg                                           >  1:222600/1‑149 (MQ=255)                                                                                                       gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACA‑GC‑ACACTTAACAGCGGACCGGCCAAGCGGGGTCTCTCCCCTTCGGCGTTTCCTGCCCATCACCACCCCTTTCGGGACGAGGGGGGTACGTTTGCGCTCAAccc                                             >  1:913964/1‑149 (MQ=255)                                                                                                       gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCgg                                           <  1:449841/151‑1 (MQ=255)                                                                                                        tgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCgg                                           <  1:56004/150‑1 (MQ=255)                                                                                                          cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGCGGGGGCC‑GA‑ACATTCGACACCAGGCCCGCCACGCGGGGTCTCTCCGCTTCGGCCTATCTTCCCCAATCCGACCGTTTGCGCAAAGAGGGGAGTACGTTGGCGGTaaaacccgt                                          >  2:705676/1‑142 (MQ=255)                                                                                                            gAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCC‑GC‑ACAATTAAAAGCGGGCCGGCAACGCGGGGTTTTTCCCTTTCCGCTTATTCTGCCCATTACGACCGCTTGCGTAAAGAGGGGGGTAAGTTGGCGGCAAACCCCGGGCgg                                      >  1:872762/1‑151 (MQ=255)                                                                                                               aGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACA‑GC‑ACAATGAACAGCGGGCCGGCAACGCGGGTTTTCTCCCCTCCGGCTTATTCTGCACATTACGAACGTTTTCGTAACGAGGGGGGTATGTTGGCGCCAAGccccgggggg                                      >  1:1271071/1‑141 (MQ=255)                                                                                                                     gATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgcgc                             <  1:1086809/151‑1 (MQ=255)                                                                                                                        cAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCaa                           <  1:1082841/150‑1 (MQ=255)                                                                                                                            cgcgGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCt                       <  2:1530778/150‑1 (MQ=255)                                                                                                                               ggTAGTGGCGGGCGAGGTTACA‑GC‑ACATTGAACCGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATGACGACCGTTTGCGTAACGGGGGGATTACGTTGGCGGTAAGCCCGGGGCGGCTGGGGCGCAATCCTCa                     <  2:1448605/149‑1 (MQ=255)                                                                                                                                  aGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCtt                 <  1:643512/150‑1 (MQ=255)                                                                                                                                  aGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTg                <  1:734968/151‑1 (MQ=255)                                                                                                                                     ggcgggcgAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCa              <  2:1357058/150‑1 (MQ=255)                                                                                                                                      gcgggcgAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAgg            >  2:137283/1‑151 (MQ=255)                                                                                                                                        gggcgAGGGTACA‑GC‑ACATTGACCAGCGGACCTGCAACGCGGATTCTCTCCGGTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGAGATTACGTTGGCGATAAGCCCGGTGCGGCTGGGGCACAATCCTCATCTTGCAGGTa          <  2:1126220/151‑1 (MQ=255)                                                                                                                                        gggcgAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTa          <  1:1507099/151‑1 (MQ=255)                                                                                                                                         ggcgAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTAt         <  2:666583/151‑1 (MQ=255)                                                                                                                                            gAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa      <  1:1205048/151‑1 (MQ=255)                                                                                                                                            gAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGCGCGGCTGGGGCGCAATCCTCATCTTGCCGGTATTa       >  2:1486679/1‑150 (MQ=255)                                                                                                                                               ggTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTACACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAc    >  2:452435/1‑150 (MQ=255)                                                                                                                                                gTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAt  >  1:798205/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                        CAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACA‑GC‑ACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  NC_000913/1153429‑1153722

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