Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,138,152	C→T	T261I (ACC→ATC)	flgJ →	putative peptidoglycan hydrolase FlgJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,138,152	0	C	T	95.1%	113.2 / NA	43	T261I (ACC→ATC)	flgJ	putative peptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (16/23);  minor base G (1/1);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAAC  >  NC_000913/1137905‑1138398                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       gcACCGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAAGGCGTAACCCGCGCTACGCCGCCGTGACGAtcgc                                                                                                                                                                                                                                                     >  2:165660/1‑251 (MQ=255)     tGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCcgcg                                                                                                                                                                                                                                                 >  1:310895/1‑251 (MQ=255)        aattgggTTGGGGGCACCGGCAAACCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTGGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGt                                                                                                                                                                                                                                              <  2:425692/246‑1 (MQ=255)              gTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGaa                                                                                                                                                                                                                                        <  1:31810/251‑1 (MQ=255)                      aaaCGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGTTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGgc                                                                                                                                                                                                                                <  1:165660/250‑1 (MQ=255)                          ggaaaaTCCGCCGCGAAAACGGCGAGCCGAGCTATAACCGGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAATTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAg                                                                                                                                                                                                                            <  1:438589/248‑1 (MQ=255)                           acAAACCCGCCGCGAAAACGGCGAGCCGATCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCCCCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAgg                                                                                                                                                                                                                           <  2:441234/250‑1 (MQ=255)                               aTCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGAGCGCCGCGAGTGCGGAACAGGGGGAGCAGGcccc                                                                                                                                                                                                                       >  2:409498/1‑250 (MQ=255)                                        gAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTGCGCCGCCGCGACGATCGCCGCGCGTCCGCAACAGGGAGCGCAgcacccaccggaacc                                                                                                                                                                                                              >  2:105839/1‑237 (MQ=255)                                               gCGAGCCGGGCTATAACCTGTTGGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGAGCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTAt                                                                                                                                                                                                       <  2:310895/251‑1 (MQ=255)                                                   gccgagcTATAACCTGTTTGGTGTCAAAGCCTCTGGAAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTCCGCGTCTACAGCCC‑GTATCTGGAATCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCa                                                                                                                                                                                                   <  2:289371/251‑1 (MQ=255)                                                       atgatctcccTGTTTGGGGTCAAAGCCTCGGGCAACTGGAAAGGGCCAGTTACCGAAAGCACCACGACTGAATAGAAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTGGCGTCTACAGCCC‑GTACCTGGAAGCCTGGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCAcc                                                                                                                                                                                                 <  2:145953/241‑1 (MQ=255)                                                              aCCTTTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGACAGTTACCGAAATCACCACGACTGAATATGAAAACGGCGAAGCGGAGAAAGTAACAGCGCAGATTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCCACGCCGCCGTGACGACCGCCGCGGGGGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTAAGCCAccaggtcgga                                                                                                                                                                                         >  1:145953/1‑242 (MQ=255)                                                               ccTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACt                                                                                                                                                                                       >  1:124205/1‑251 (MQ=255)                                                                                     gtcacctgtaaaGGGCCAGTTACTGAAATCACCACGATTGAATATGAATACGGCGAATCGATGAAAGTAAAAGCGAAGTTGCGCGTCTACAGCTC‑GTACCTGGAAGCCTTGCCGCATAACTTTGGGTTGTTAACGCGTAACCCGCGCAACGCCGCCGTGACGATCGCCGCTAGTGCGGAACAGGGGGCGCAGTCCCTACAGGACGCGGTCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACAt                                                                                                                                                                 <  1:88854/242‑1 (MQ=255)                                                                                            ggAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTca                                                                                                                                                           <  1:178435/250‑1 (MQ=255)                                                                                                           cTGAAATCACCACGCCTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGCAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCAACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGAt                                                                                                                                            >  2:65019/1‑250 (MQ=255)                                                                                                           cTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGaa                                                                                                                                                                                                                                        <  1:181856/158‑1 (MQ=255)                                                                                                           cTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGaa                                                                                                                                                                                                                                        >  2:181856/1‑158 (MQ=255)                                                                                                                  caccacGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAACCTCACCAACATGATTCAGCAGATGAAATCGATAAGCGAc                                                                                                                                     >  2:279588/1‑250 (MQ=255)                                                                                                                           tGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGc                                                                                                                            >  1:166437/1‑250 (MQ=255)                                                                                                                                tggAAAACGGGGAAGGGATGATAGTAAAAGCGAATTTTCGCGTTTACAGCTC‑GTTTCTGGAATCCGTGTCGGATTCCGTGGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCCCTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAAc                                                                                                                       <  2:381017/248‑1 (MQ=255)                                                                                                                                atGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAAcc                                                                                                                      <  1:64198/251‑1 (MQ=255)                                                                                                                                   aaaaCGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAACCTCACCACCATGATTCAGCAGATGAAATCGATAAGCGCCAAGGTGAGCAAAACCTAc                                                                                                                   >  2:215788/1‑251 (MQ=255)                                                                                                                                    aaaCGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACa                                                                                                                  <  2:166437/251‑1 (MQ=255)                                                                                                                                        ggCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGcc                                                                                                                                                                                                    >  1:256924/1‑165 (MQ=255)                                                                                                                                        ggCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTGGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGcc                                                                                                                                                                                                    <  2:256924/165‑1 (MQ=255)                                                                                                                                            aagcgaagAAAGTAAAAGGGGAGTTTCGTGTCTACATCTC‑GTATCGTGAAGACTTGTCGGATGACGTGGGGCTGTTAACGCGGATCCCGCGGTACTCCGCTGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTAGTCCCGCAAACTCACCAACATGATTCATCAGATGAAATCGATAAGCGACAAGGTGATCAAAACCTACAGTATGAAc                                                                                                          <  2:339891/251‑1 (MQ=255)                                                                                                                                            aacagaagaaATGAAAAGCGAAGTTTCTCTTCTTCAG‑TCTGTATTTTGAAGCCTTGGCGGAGTACTTTGGGTTGTTAACGCGTAACACCCGCTTCGCCTCCGGGAATATCGCCTCTAGTGCGTACCAGGTGGCGCAGGCCCTACAGGACGAGGTCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCTACAATGTGATCAAAACCTACAGTTTGAAc                                                                                                          <  2:339879/247‑1 (MQ=255)                                                                                                                                               cgaagAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACAtt                                                                                                       <  1:105839/251‑1 (MQ=255)                                                                                                                                                 aagaaAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTAGCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGa                                                                                                     <  1:52299/251‑1 (MQ=255)                                                                                                                                                  agaaAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGAt                                                                                                    >  1:189742/1‑251 (MQ=255)                                                                                                                                                       gTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGCACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGCTTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCt                                                                                               >  2:455621/1‑251 (MQ=255)                                                                                                                                                        aaaaaGGGAAGGTACGCGTCTACTGCTC‑GTATCTGGAAGCCTTGCCGGATGACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTg                                                                                              <  2:14087/250‑1 (MQ=255)                                                                                                                                                           aaGCGAAGTTTCGCGTCTACAGCTC‑GTTCCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTACCCCGCGCTACGCCGCCGTGACGATCGCCGCTAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAACCCTACAGTATGAACATTGATAATCTGTTc                                                                                           <  2:124205/251‑1 (MQ=255)                                                                                                                                                                          tCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTc                                                                            >  1:86975/1‑251 (MQ=255)                                                                                                                                                                            tACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTa                                                                                                                    <  2:408337/209‑1 (MQ=255)                                                                                                                                                                            tACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTa                                                                                                                    >  1:408337/1‑209 (MQ=255)                                                                                                                                                                                           ggAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCt                                                                <  1:221663/247‑1 (MQ=255)                                                                                                                                                                                                   tGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATAc                                                     <  1:256752/250‑1 (MQ=255)                                                                                                                                                                                                      cGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGGATGAACATTGATAATCTGTTCTGAATAACTCAAGACCGGCGGGTCGCTGCCggacatactatg                                                 >  2:147807/1‑240 (MQ=255)                                                                                                                                                                                                                                aaCCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTc                        <  1:129639/250‑1 (MQ=255)                                                                                                                                                                                                                                                ccgTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGAtt        >  1:281998/1‑250 (MQ=255)                                                                                                                                                                                                                                                  gTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAACTCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGACGGCTTATAAGGAACCTCCATGTCCAGCTTGATtaat     >  2:151438/1‑251 (MQ=255)                                                                                                                                                                                                                                                     acgatCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAAc  <  1:3771/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       GCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTC‑GTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAAC  >  NC_000913/1137905‑1138398

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