Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,446,775	C→T	A233A (GCG→GCA)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,446,775	0	C	T	97.5%	125.8 / ‑4.4	40	A233A (GCG→GCA)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (18/21);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

ATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGA  >  NC_000913/2446530‑2447016                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    atcatcTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCTGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATg                                                                                                                                                                                                                                                >  2:143707/1‑251 (MQ=255)    agctcAACTTCCTCGCCAAAGCGGTTAATGGTACGCCCCGGCACGGTAATGCTGGAGGTCGGTTGCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGc                                                                                                                                                                                                                                             <  2:154650/249‑1 (MQ=255)       tcAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGAt                                                                                                                                                                                                                                          >  1:235657/1‑251 (MQ=255)         aaCTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGt                                                                                                                                                                                                                                        >  1:45622/1‑251 (MQ=255)                    aaaaGGGGTTATTGGTGCGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATAGGGGCAATGATGGGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAgg                                                                                                                                                                                                                             <  2:54816/250‑1 (MQ=255)                          ggTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTc                                                                                                                                                                                                                        >  2:143191/1‑250 (MQ=255)                                         ccGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCa                                                                                                                                                                                                        <  1:57517/251‑1 (MQ=255)                                             cACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTc                                                                                                                                                                                                    <  1:236113/251‑1 (MQ=255)                                                      gCTGGAGGTCGGTTTCAGCGCCATATGGGCAAGGATTTGCTGCCCGCTGCTGATGCCGCCGAGAATGCCTCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGaa                                                                                                                                                                                           <  1:165077/251‑1 (MQ=255)                                                                     aagcgcgttgggggCAATGATTTGTGGCCCCCTGAGGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCTGTTCTGGCTGCCGCGCAGCGCCACCACGGCAAAGCCGTCGCCAATTTCCACGCCTTTCTCCGC‑GTTGATGCTCATCAGTGCATGGGCGATTTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGcaa‑caa                                                                                                                                                                           <  2:142320/240‑1 (MQ=255)                                                                      agagcCATATGGGCAATGATTTGTTGCCTGCTGCTGTGTCCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCTCCAAGTCCGGCAGGAACGCCACTGGCaa‑caac                                                                                                                                                                          <  2:72144/248‑1 (MQ=255)                                                                      aGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCCGGAACGCCACTGGcaa‑caa                                                                                                                                                                           >  1:85477/1‑250 (MQ=255)                                                                              atGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGCACGCCACTGGCAA‑CAACGGTGACtt                                                                                                                                                                  >  1:12814/1‑251 (MQ=255)                                                                                 ggCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGCGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCCAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGGGACttt                                                                                                                                                                 >  1:30927/1‑249 (MQ=255)                                                                                       gATTTGCGGCCCGCTGCTGATGCCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATGTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACGTGCAA‑CAACGGTGACTTTAGCGCCGa                                                                                                                                                         <  2:179198/251‑1 (MQ=255)                                                                                               gcCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGt                                                                                                                                                   <  1:136572/249‑1 (MQ=255)                                                                                                   gctgctGATACCGCCGAGCATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCCAAGACCGGCTCCCCAAGTCCGGGAGGAAAGCCCCCGGCAAACAA‑GGTGGCTCTAGCGCCGATCGAGTCGcccc                                                                                                                                             >  2:189824/1‑250 (MQ=255)                                                                                                                            cccgcATGGTTGCTCTGGAAACCGTCTTGGGTGATTGCATCGCCGTTCTGGCTGCCGCGCAGCGCCCCCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTTATGCTCATCAGCGCAGGGGCGATGTCAGCATCCAGGCGTTCAAAGACCGGCTCGCCAAGTCCGTCAGGAACGCCACTGGCAC‑CAATGGGGACTTTAGCGCCGATGGAGTCGCCCTCGTTTTT‑CAGCGCATGCATCAACTc                                                                                                                   <  2:162425/251‑1 (MQ=255)                                                                                                                                  tGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCCAGTCCGGCAGGAAAGCCACCGGCCA‑CAAAGGGGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑TAgccccacactccactcgtctt                                                                                                              >  2:147448/1‑231 (MQ=255)                                                                                                                                     ttGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGACTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCCAGTTCGGCAGGAACGCCACTGGCAA‑CAACGGCGACTTTAGCGCCGGTGGGGTGCCCCTCTTTTTTTCA‑CGCACGCCTACACCCGTCTAAcgc                                                                                                          >  2:126435/1‑251 (MQ=255)                                                                                                                                          ctGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTc                                                                                                       <  1:65425/249‑1 (MQ=255)                                                                                                                                              aaaCCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGAtttt                                                                                                 <  2:235657/251‑1 (MQ=255)                                                                                                                                                          gTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCgg                                                                                     <  2:12814/251‑1 (MQ=255)                                                                                                                                                           tGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATg                                                                                                                                                       <  1:260593/185‑1 (MQ=255)                                                                                                                                                                                       cAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGGGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCg                                                        >  1:23617/1‑251 (MQ=255)                                                                                                                                                                                           gcCACCACGTCAAAGCCGTAGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCCAGTCCGGCAGGAAACCCACTGGCAA‑CAACGGGGACTTTTGCGCCTGTGGAGTCCCCCTCTTTTTT‑CAGCTCACGCATCAACTCGTATAACGCGTCGATTTTTTCTGGGACTGGGCAAAAAAACGGAtttttctttaccgccgacc                                                     >  2:225475/1‑235 (MQ=255)                                                                                                                                                                                              cccacGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTc                                                 <  2:85477/250‑1 (MQ=255)                                                                                                                                                                                              cccacGTCAAAGCCCTCGCCAATTCCCACGCCTTTCACCGC‑TTTGATGCCCACCAGTGCATGGGCGATGCCCGCATCCAGGCGGTCAAAGTCCGGCTCGCCAAGCCCGGGAGGAACCCCACTGGCAA‑AAACGGTGCCTTTAGCGCCCATGAAGTCGCCCCCTTTGTT‑CAGCGCACGCACCAACTCGTCTATCCCTTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCACGCTCTGCAACCAGTc                                                 <  2:30927/250‑1 (MQ=255)                                                                                                                                                                                               ccacGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCt                                                <  1:126435/251‑1 (MQ=255)                                                                                                                                                                                                cacGTCAAAGCCGTCGCCAATTTCCACGCCTTTCAC‑GC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCttt                                              >  1:64865/1‑251 (MQ=255)                                                                                                                                                                                                         gCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATcc                                       >  2:156474/1‑250 (MQ=255)                                                                                                                                                                                                            gTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCg                                   <  1:189824/251‑1 (MQ=255)                                                                                                                                                                                                               gCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGCTGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCCCTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTTTCTGGGTCCGGGCAAAAAAACGGATTTTTCTCGCCCTGCCACCACTCTTTGATATCCCGCGGa                                 >  2:21353/1‑250 (MQ=255)                                                                                                                                                                                                                        cACGCCTTTCGCCGCTGTT‑ATTCTCATCATTGCCTGGGCGATGTCAGCATCCAGTCGGTCAAAGACCGGCTCGCCAAGTCCGGCATGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGTAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCa                       <  1:10977/251‑1 (MQ=255)                                                                                                                                                                                                                         aCGCCTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACCGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGAACTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCAt                      >  2:103755/1‑251 (MQ=255)                                                                                                                                                                                                                            ccTTTCGCCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGAGTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTg                   <  2:159123/251‑1 (MQ=255)                                                                                                                                                                                                                             cTTTCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTAAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGAGAATTTAGCGCAGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATGAACTCGTCTAACGCGTAGATTTTGTCGGGGTCCGGGCAAAAAAAAGGGTTTTGCTCGACCTTAGACCAGTCTTTGATATCCAGCGGAATGTAGCACATCTgg                  >  1:35949/1‑251 (MQ=255)                                                                                                                                                                                                                                tCACCGC‑GTTGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTggggc               >  1:138960/1‑249 (MQ=255)                                                                                                                                                                                                                                       cggtGATGCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAg         <  2:23617/248‑1 (MQ=255)                                                                                                                                                                                                                                              gCTCATCAGTGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGa  <  1:192880/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGC‑GTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAA‑CAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTT‑CAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGA  >  NC_000913/2446530‑2447016

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