Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,130,803	C→T	A21T (GCG→ACG)	flgA ←	flagellar basal body P‑ring formation protein FlgA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,130,803	0	C	T	95.0%	101.9 / NA	42	A21T (GCG→ACG)	flgA	flagellar basal body P‑ring formation protein FlgA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (14/24);  minor base G (1/1);  total (16/26)
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

ATCGCGGCAACCACATAATTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGCCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGT  >  NC_000913/1130559‑1131047                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      aTCGCGGCAACCACATAATTCCCTGTGGCCTGTACATTAACCGGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGATATTGCTCGTCGTACt                                                                                                                                                                                                                                                 <  1:331532/251‑1 (MQ=255)     cgGCAACCACATAATGTCCTGGGGCCGGTACATTAACCTGTAAATATCGTGTGTCGGTACCGCAGCGTGCCAACACATTAACATTGCCCCACAGGCGGG‑AAGTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCAGCACGTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATGGCGATGTGAGATTGCTCGTCGTGCTTAAc                                                                                                                                                                                                                                             <  2:276400/251‑1 (MQ=255)       gCAACCACATAATTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGGAATT‑TTGCTCATCGGAAGAAATGGCTGCTCGCATGGCGGGAGTAGATTGGGCGCTGTAAGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTTTGCAATTGCGATGTGAGAATGCTCGTCGTACTTaaagg                                                                                                                                                                                                                                           >  1:46981/1‑248 (MQ=255)                                       aaCCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTAc                                                                                                                                                                                                           <  2:399697/251‑1 (MQ=255)                                               aaTATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTTCGTTttatt                                                                                                                                                                                                   <  2:37129/251‑1 (MQ=255)                                                atatCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTAGCGCGAGGATCGCCACGCTACGTTttatta                                                                                                                                                                                                  <  1:169011/251‑1 (MQ=255)                                                     gTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCCCTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACCCTACGTTTTATTAT‑CAGc                                                                                                                                                                                            >  1:269002/1‑251 (MQ=255)                                                           cGTTACCGCATCGTGCCAACACATTCACATTGACCCACAGGCGGT‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCTTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTc                                                                                                                                                                                      <  2:62109/251‑1 (MQ=255)                                                              tACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCCTCGAAAAAAAAGGCTCCTCGCAAGGCGGAAGGGGATTGGGCGCTGTACGAATAGAAAAACCAACCTCCTCACTTTCCCCCGCGCGGTGGGCGCTTAAAAAGATGTTCAAATGCGCTGTGAGATTGCGCGGCGCGCTTTAAGGGCTGAACAGAATCGCGATGATCGGCACGCTACGCTTT‑TTATTCAGCAttttatcc                                                                                                                                                                                   >  1:303767/1‑247 (MQ=255)                                                                ccGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGcccc                                                                                                                                                                                 <  2:66858/251‑1 (MQ=255)                                                                 cGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGATGGGCAATTGCGATGTGAGATTGCTCGTCGTACTTTACGGACTGAACAGTATCGCGATGATCGCCACGCCACGGTTTATTTT‑AAGCATTTTTGccccc                                                                                                                                                                                >  1:435606/1‑251 (MQ=255)                                                                           cAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCt                                                                                                                                                                      <  1:392398/251‑1 (MQ=255)                                                                                    cacattGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGa                                                                                                                                                             <  1:126495/251‑1 (MQ=255)                                                                                    cacattGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGa                                                                                                                                                             <  1:212415/251‑1 (MQ=255)                                                                                                 aGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTcc                                                                                                                                                >  1:312579/1‑251 (MQ=255)                                                                                                      gg‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAAt                                                                                                                                           >  1:101069/1‑251 (MQ=255)                                                                                                         aaTTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATg                                                                                                                                          >  1:69836/1‑250 (MQ=255)                                                                                                             attGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAAt                                                                                                                                                           >  1:583433/1‑229 (MQ=255)                                                                                                             attGCTCATCGAAAGCAATGGCGGCTCGCATGGCGGTATGAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGAGTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCTCGCTACGTTTAATTAT‑CAGCATTTTCGCCCCCAGCCATGTCTACAACGTGAAt                                                                                                                                                           <  2:583433/229‑1 (MQ=255)                                                                                                                      cGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCAGCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGCGAATTGTACCTTTCCGCAATGACCATCAACGGCAt                                                                                                                            <  2:559066/251‑1 (MQ=255)                                                                                                                       gAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCgat                                                                                                                                                                                                                         <  2:598579/158‑1 (MQ=255)                                                                                                                       gAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCgat                                                                                                                                                                                                                         >  1:598579/1‑158 (MQ=255)                                                                                                                               tGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCg                                                                                                                    >  1:151979/1‑250 (MQ=255)                                                                                                                                         cATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCACAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGc                                                                                                         <  1:518195/251‑1 (MQ=255)                                                                                                                                         cATGGCGGTAGTAGATGGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTATCTGTCCGCAATGACCAGCAACGGCATAAATAGCGACCCATTTTGc                                                                                                         <  2:45029/251‑1 (MQ=255)                                                                                                                                            ggCGGTAGTAGATTTGGCGCTGTACGGATTGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGTCGCTAAAAAAGTTGTGCAATTGCTAGGTGAGATTGCTCGTCGTTCTTAACGGACTGAACAGTAGCGCGATGATCGCCACGCTACGTTTTATTCT‑CAGCATTTTCGCCCCCAGCCAGTGCTTCAACGTGCATTGTACCTGTCCTCAATGACCATCAACTGCATAAATATCGACTCATTTTGCGt                                                                                                       <  2:414556/250‑1 (MQ=255)                                                                                                                                             gCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGttt                                                                                                     <  1:453172/251‑1 (MQ=255)                                                                                                                                                        ttGGGCGCTGTTCGAATCGCAACACGACCCTCATCACTTACCCCCCTGAGTTGGGCGCTAAAAAAGTTGTGCACTTGAGATCTGAGCCTGCTCGTCGCACTTCACGGACTGAACAGTAACGCGAAGATCGCCACGCTACGTGTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTCCAACGTGAGTTGTCCCTGTCCGCAATGACCATCAACGGCATCAATAACGACCCATTTTGCGTTTATTCCGCCGAt                                                                                          >  2:556889/1‑251 (MQ=255)                                                                                                                                                                      aaaaaaaaCACGAACCTCATCACTTGCCCCCGCGGGGTTGGCGATAAAAAAGTTGGGCAATGGCGAGGTGAGATCGCTCGTCGTACTTAACGGACTGAACATTATCGCGATGATCGTCACGCTACGTTTTATTAT‑CATCATTTTCGCCCCCAGCCATTTCTACAACGTGAATGGTACCGGTCCGCAATGACCATCAACGGCATAAATAGCGACCCAGTTTGTGTTAATTCCGCCGATAACGCGCTCGTaaa                                                                            <  2:303767/246‑1 (MQ=255)                                                                                                                                                                            aaCACGAACCTCATCTCTTACCCCCGCGAGTTTGGCGCTAAAAAAGTTGGGCAATTGCGATGTGAGATTGCGCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACTCTACTTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTTTACCTGTCCGCAATGACCATCAACCGCAAAAATAGCGACCCTTTGTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCAtg                                                                      <  2:381376/251‑2 (MQ=255)                                                                                                                                                                             acacGAACCACATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGGGAGATTGCCCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCAttt                                                                     <  2:101069/251‑1 (MQ=255)                                                                                                                                                                               acGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGAGGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCAACGCTACGTTTTATTAG‑CAGCATTTTCGCCCCCAGCCATTTCTACAACTTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTaa                                                                   <  1:41609/251‑1 (MQ=255)                                                                                                                                                                                         tcaCTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGc                                                          >  1:296506/1‑250 (MQ=255)                                                                                                                                                                                              taccccatagAGTTGCGCGCTCAAAATGTTGTGCAATTGAGCTGTGAGATACCTCGTCGCACTTAACGGATTGAACAGTATCGCGATGATAGCCACGCTACGTTTTATTAT‑AAGCATTTTCACCCGCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAACAAGCGACCCATTTTGCGTTTATTCCGCCGATACGGCGCGCGCAAAGGCTTTTAAGATGATGGCAgcctt                                                    >  2:569571/10‑247 (MQ=255)                                                                                                                                                                                                cccccGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGCCGTACTTAACGGACTGAACAGTATAGCGATGATCGCCACGCTACGTTTTATTAT‑AAGCATTTTCGCCCCCCGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAAACGCGACCCATTTTGAGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAAtttt                                                  >  2:203734/1‑251 (MQ=255)                                                                                                                                                                                                 ccccGCGAGGTGGGCGCAAAAAAAGTTGTGCAATTGCGAGGTGAGATTGCTCGCCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTg                                                 <  2:473367/251‑1 (MQ=255)                                                                                                                                                                                                         gTTGGGCGCTAAAAAAGTTGTGCAATTGCGAGGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGc                                         <  2:93888/251‑1 (MQ=255)                                                                                                                                                                                                           ggggCGCTAAAAAAGTTGTGCAATTGCGATGGGAGATTGCTCGTCGTACTTAACGAACTGACCAGTACCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAAGGACAATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCTCCGATAACGCGCGCGTAAAGGCATTTAAGCTGAGGGCAGAATTTTGATACCTGCgg                                       <  2:468944/250‑1 (MQ=255)                                                                                                                                                                                                                   aaaaaaGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGat                                >  1:620352/1‑250 (MQ=255)                                                                                                                                                                                                                             ggCAATTGCGATGGTAGATTGCTCGTCGTACTTAACGGACTGAACATTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCGCGATaa                     <  2:544049/250‑1 (MQ=255)                                                                                                                                                                                                                                         gtgAGATTGCTCGGCGTACTTAACGGACTGAACAGTATCTCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGAcgccg          <  2:269002/250‑1 (MQ=255)                                                                                                                                                                                                                                                tGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGCCCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAACGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGt  >  2:47956/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      ATCGCGGCAACCACATAATTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGG‑AATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGCCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTAT‑CAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGT  >  NC_000913/1130559‑1131047

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: ‑