Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	6,809	T→G	K384T (AAA→ACA)	yaaJ ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	6,809	0	T	G	98.6%	248.3 / NA	72	K384T (AAA→ACA)	yaaJ	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (36/35);  total (36/36)
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.98e-01

ATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATG  >  NC_000913/6578‑7033                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                   atatCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGGGCCGCCGATGACCGCTGCGAAGGTGCAGATCCCCAAACACCAGATCGCTGTAgg                                                                                                                                                                                                                                >  1:138238/1‑236 (MQ=255)       ggATTACGCAATGGATCAAACACGGGCCCCACGCCGAGTTTACGCTGGCGTAGATAAGCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGg‑ttgtt                                                                                                                                                                                                                         <  2:125885/236‑1 (MQ=255)        gATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAAAACCAGATCGCTGTAggggt                                                                                                                                                                                                                             >  1:377544/1‑230 (MQ=255)                  tGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTaaaa                                                                                                                                                                                                              >  1:26510/1‑236 (MQ=255)                       cAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAaagaag                                                                                                                                                                                                         >  1:62060/1‑236 (MQ=255)                       cAAACACAGGTAGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAaagaag                                                                                                                                                                                                         <  2:255456/236‑1 (MQ=255)                         aaCACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGTAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAaagaa                                                                                                                                                                                                          <  1:8110/233‑1 (MQ=255)                              cGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTg                                                                                                                                                                                                   >  1:363704/1‑235 (MQ=255)                                     aCGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCGTTATGATATCTGCCAGTGGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCgg                                                                                                                                                                                            <  1:141021/235‑1 (MQ=255)                                     aCGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCgg                                                                                                                                                                                            <  1:724723/235‑1 (MQ=255)                                       gCCGATTTTACGATGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCa                                                                                                                                                                                          <  1:145460/235‑1 (MQ=255)                                       gCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGTGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCa                                                                                                                                                                                          <  1:183354/235‑1 (MQ=255)                                                 cGCTGGAGTAGATAAGCCCTGGCAATGGTAGGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGtt                                                                                                                                                                               <  2:74149/236‑1 (MQ=255)                                                      gCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGGGCCGCCGATGACCGTTGCGAAGGGGCAGATCCGCAAACACCAGATCGCGGTAGGG‑TTGGTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGGaa                                                                                                                                                                          >  1:741467/1‑236 (MQ=255)                                                         tAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAgaga                                                                                                                                                                                                      <  1:263118/205‑1 (MQ=255)                                                               aTCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCgg                                                                                                                                                                                            <  1:347813/209‑1 (MQ=255)                                                                tCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGtttt                                                                                                                                                                                               >  2:516744/1‑205 (MQ=255)                                                                tCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGtttt                                                                                                                                                                                               <  1:516744/205‑1 (MQ=255)                                                                tCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCt                                                                                                                                                                >  1:669514/1‑236 (MQ=255)                                                                tCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGAACGTTGCGAAGGTGCAGAACCGCAAAAACCAGAACGGTGTAGGG‑GTGTTCAGGGGTAAAAAGAAGAGATAGattt                                                                                                                                                                                               >  2:263118/1‑205 (MQ=255)                                                                    tGGCAATGTTAGGCACTAGAGTCGACATCCGTGAAATGGC‑GTCATAGTTGGTAATAGCCATGCAGGACATTATGATCTATGCCAGTTCCCACATCAACGGAAGGCTTGGCAATGGGCCGCCGATGGCCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCTCTGTAGGG‑TTTTTCAAGCATAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑g                                                                                                                                                           <  2:449091/236‑1 (MQ=255)                                                                     ggCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑gg                                                                                                                                                          >  1:551618/1‑236 (MQ=255)                                                                      gCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGGAATAGCCATGCAGGCCATTATGATATCTGCCAGGTGCCACATCAGCGGAAGGGTTAGCAAGGGGCCGCCGATGAACGTTGCGAAGGGGCAGAACCGCAAAAACCAGAACGCTGTAGGG‑TTGGTCAGGCGTAAAAAGAAGAGATTGGTTTAGGCAAAAAt                                                                                                                                                                                     >  2:347813/1‑209 (MQ=255)                                                                        aaTGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑g                                                                                                                                                           <  2:26510/232‑1 (MQ=255)                                                                          tGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGGCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑GTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAAc                                                                                                                                                     >  1:512861/1‑236 (MQ=255)                                                                             tATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGa                                                                                                                                                  <  2:363704/236‑1 (MQ=255)                                                                                       gggctgtgtcAGTAAAAGGGCGTTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACaa                                                                                                                                        <  2:741467/227‑1 (MQ=255)                                                                                         gCGAGAGCAGTAAAATGGCGGTCA‑AGTTGATAATAGCCAGGCAGGCCATTATGATAGCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGAAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAgg                                                                                                                                      <  2:138238/236‑1 (MQ=255)                                                                                           gagtGCAGTAAAATGGCCGTCA‑AATTGATAAGAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGa                                                                                                                                                  <  2:23172/222‑1 (MQ=255)                                                                                            agagCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAGACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCa                                                                                                                                                                                          <  2:695426/182‑1 (MQ=255)                                                                                            agagCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCa                                                                                                                                                                                          >  1:695426/1‑182 (MQ=255)                                                                                            agagCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTa                                                                                                                                   >  1:449456/1‑236 (MQ=255)                                                                                             gagCAGTAAAATGGCGGTAA‑AATGGGAAATAGCCATGCAGGCCATTAGGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTa                                                                                                                                   <  2:669514/235‑1 (MQ=255)                                                                                                 aGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCTACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGCAGAGATTGTATTCGGCATAAATGTAGTTGGCAAAGATGGAGCTGAAGGGCAAAC‑GAATAAACACAAGGGGAACaaa                                                                                                                              >  2:450458/1‑236 (MQ=255)                                                                                                  gTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACaaa                                                                                                                              <  1:510372/235‑1 (MQ=255)                                                                                                       aTGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAg                                                                                                                         <  2:551618/235‑1 (MQ=255)                                                                                                             gTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAATGGTAACAAACTCAGCAcccc                                                                                                                   >  1:59727/1‑235 (MQ=255)                                                                                                             gTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGGAACAAAATCAGCACCCCa                                                                                                                  >  1:418761/1‑236 (MQ=255)                                                                                                                   aTTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAgg                                                                                                                <  2:449456/233‑1 (MQ=255)                                                                                                                     tGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCTACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAAccc                                                                                                           <  1:450458/236‑1 (MQ=255)                                                                                                                             gCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTc                                                                                                                           <  2:138904/212‑1 (MQ=255)                                                                                                                             gCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTc                                                                                                                           >  1:138904/1‑212 (MQ=255)                                                                                                                                 tGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCg                                                                                                >  2:433230/1‑235 (MQ=255)                                                                                                                                  gCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCaaa                                                                                                                                                      <  1:596350/180‑1 (MQ=255)                                                                                                                                  gCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCaaa                                                                                                                                                      >  2:596350/1‑180 (MQ=255)                                                                                                                                  gCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACaa                                                                                                                                        >  2:454056/1‑194 (MQ=255)                                                                                                                                  gCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACaa                                                                                                                                        <  1:454056/194‑1 (MQ=255)                                                                                                                                  gCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCAACCGc                                                                                               >  1:182875/1‑235 (MQ=255)                                                                                                                                                tatCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAaagaa                                                                                                                                                                                                          <  1:11176/115‑1 (MQ=255)                                                                                                                                                tatCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAaagaa                                                                                                                                                                                                          >  2:11176/1‑115 (MQ=255)                                                                                                                                                tatCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATa                                                                                >  1:623669/1‑236 (MQ=255)                                                                                                                                                  tCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGa                                                                                  >  2:591323/1‑232 (MQ=255)                                                                                                                                                    tGCCAGTCGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCt                                                                            >  1:4925/1‑236 (MQ=255)                                                                                                                                                         gTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAAt                                                                        >  2:495484/1‑235 (MQ=255)                                                                                                                                                          ttGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATa                                                                       >  1:402677/1‑235 (MQ=255)                                                                                                                                                              cacaTCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTc                                                                  <  1:495484/236‑1 (MQ=255)                                                                                                                                                                  tCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGc                                                              <  2:4925/236‑1 (MQ=255)                                                                                                                                                                    aGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCgg                                                            <  1:433230/236‑1 (MQ=255)                                                                                                                                                                          aGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTa                                                      >  1:633705/1‑236 (MQ=255)                                                                                                                                                                           ggCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAAc                                                                                                                                 >  2:245110/1‑160 (MQ=255)                                                                                                                                                                           ggCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAAc                                                                                                                                 <  1:245110/160‑1 (MQ=255)                                                                                                                                                                                gacagggTGCCGACGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGATTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGtt                                                  <  2:62060/229‑1 (MQ=255)                                                                                                                                                                                aGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGtt                                                  <  1:591323/234‑1 (MQ=255)                                                                                                                                                                                aGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTtgt                                                >  1:589299/1‑236 (MQ=255)                                                                                                                                                                                            caggtGACCGTGGCGAAGGTGCAGATCCGAAAACCCAAGTTCGCTGTAGGT‑TTGTTCAGGATTAAAAAGAAGAGTTTTTTTCCGGCATAAATTTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAAAACCCACAAGGATAACAAACTCAGCCCCCCAGGACCCCATTAGCCCCCGCTTCCCCTTCTGGATAAGCTGAATTCCTCCCGGCGGCATGTAGGTTGGGCCGTTCCCCg                                     <  2:512861/231‑1 (MQ=255)                                                                                                                                                                                                              gTGCAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGAAGGTTGTGCCGTTACCCCCCAGTAATATCAGCaagac                  >  1:545499/1‑232 (MQ=255)                                                                                                                                                                                                                caaGATCCGCAAACTCCGGA‑CACGTTGGGGTTTGTGCGGGCATAAAAAGAAGAGATTGTTTGCGACATAAATGGAGTTGAAAACGATGGAGCGGAA‑GCCAAACAAAATACCCAAAAGGGAAAAAAACTACGCACCCCAGGAACCCTTTGGCCCCCGCACCGCCTTCTGGATAAGCTAAATCCCTCCCGGGGGTATGGGGGTTGTGCCTTTACCCGCCAGTAATATCAGCATGgcgc                <  2:182875/234‑1 (MQ=255)                                                                                                                                                                                                                 cAGATCCGCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCCGCGGCATGTAGGTTGTGCCGTTACCCGCCCGCCATCTCAGCATGgcgc                >  1:594331/1‑235 (MQ=255)                                                                                                                                                                                                                        gCACACACAAGAGCGCTGTGGGT‑TTGTTCAGGCTAAAAAAGAAGATTTGTTTTTGGGCATAAATGTAGTTGGAAACGAGGGAGCGGAA‑GGCAAACAAAATACCCCCAAGGGAAACAAACTCACCCCCCCAGGACCCCATTAGCCCCCGCTTCGCCTTCTGGATAAGCTAAATCCCTCCCAGGGGCATGTGGGTTGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCt               <  2:402677/229‑1 (MQ=255)                                                                                                                                                                                                                        gCAAACACCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTa        >  1:417592/1‑236 (MQ=255)                                                                                                                                                                                                                         cAAACACCAGATCGCTGTAGGG‑TTGTTCAGG‑GTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACa      >  1:401735/1‑236 (MQ=255)                                                                                                                                                                                                                              acCAGATCGCTGTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGTCGCTTGCCGTAACGATg  >  1:339490/1‑236 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                   ATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCA‑AATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGG‑TTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAA‑GGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATG  >  NC_000913/6578‑7033

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