Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,133,542	A→C	E323D (GAA→GAC)	flgE →	flagellar hook protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,133,542	0	A	C	96.8%	214.1 / ‑0.2	62	E323D (GAA→GAC)	flgE	flagellar hook protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (28/32);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

ATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGT  >  NC_000913/1133394‑1133688                                                                                                                                                     |                                                                                                                                                   atatTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACgg                                                                                                                                                  >  1:1194113/1‑151 (MQ=255)      gtgGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTgg                                                                                                                                             <  1:781257/151‑1 (MQ=255)      gtgGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTg                                                                                                                                              <  2:890039/150‑1 (MQ=255)      gtgGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTg                                                                                                                                              <  1:1241003/150‑1 (MQ=255)       tgGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGc                                                                                                                                            >  1:936881/1‑151 (MQ=255)           aaCCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATcc                                                                                                                                        <  2:1135817/151‑1 (MQ=255)            accaccCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCg                                                                                                                                       <  1:989035/151‑1 (MQ=255)              caccCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGa                                                                                                                                      >  2:925630/1‑150 (MQ=255)               accCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAg                                                                                                                                    <  2:528408/151‑1 (MQ=255)                   aGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCg                                                                                                                                 >  1:777466/1‑150 (MQ=255)                   aGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCg                                                                                                                                 >  2:298209/1‑150 (MQ=255)                        ggCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACg                                                                                                                           <  1:308529/151‑1 (MQ=255)                          cTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACgtc                                                                                                                         >  1:741243/1‑151 (MQ=255)                              aaaCCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACgtctggt                                                                                                                     >  1:830017/1‑151 (MQ=255)                                aCCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACgtctggtct                                                                                                                   <  2:402219/151‑1 (MQ=255)                                aCCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGGCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACgtctggtc                                                                                                                    >  1:533764/1‑150 (MQ=255)                                       gATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGc                                                                                                            >  2:241381/1‑151 (MQ=255)                                          cTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAAt                                                                                                         >  1:262183/1‑151 (MQ=255)                                                 gTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCtgg                                                                                                  <  1:786827/151‑1 (MQ=255)                                                      caaatcaaTGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCtggcgtgg                                                                                             >  2:1188714/1‑151 (MQ=255)                                                      caaatcaaTGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCtggcgtgg                                                                                             <  2:870816/151‑1 (MQ=255)                                                        aatcaaTGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCtggcgtggc                                                                                            >  1:1107982/1‑150 (MQ=255)                                                             atgatgACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGtt                                                                                      >  2:146642/1‑151 (MQ=255)                                                                 tgaCGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTgggg                                                                                  >  1:599361/1‑151 (MQ=255)                                                                   acggtacggtTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGAc                                                                                <  2:185106/151‑1 (MQ=255)                                                                    cggtacggtTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGAc                                                                                >  2:555408/1‑150 (MQ=255)                                                                    cggtacggtTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACa                                                                               <  1:628067/151‑1 (MQ=255)                                                                        acggtTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCcg                                                                           <  2:1194186/151‑1 (MQ=255)                                                                                cGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCcgggacggg                                                                    <  2:119625/150‑1 (MQ=255)                                                                                cGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCcgggacggg                                                                    <  1:937844/150‑1 (MQ=255)                                                                                     aCTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACtt                                                              <  2:720646/151‑1 (MQ=255)                                                                                        aTTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTgg                                                           <  2:220922/151‑1 (MQ=255)                                                                                        aTTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTg                                                            <  1:942501/150‑1 (MQ=255)                                                                                             aaCGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCAccc                                                      <  1:802244/151‑1 (MQ=255)                                                                                              aCGAACAAACCCACCTGCTGGGGCAGATTGTCCTGGCGAACTTTGCTAACAACGACGGTCTGGCATCCGAAGGCGACAACGCCTGGTCTGCGACGAAATCTTCTGGCGTGGCGCTGTTGGTGACAGCCGGGACGGGAAACTTTGGCAccc                                                      <  1:1749/150‑1 (MQ=255)                                                                                              aCGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCt                                                     >  1:266541/1‑151 (MQ=255)                                                                                                 aaCAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGAc                                                  >  1:371501/1‑151 (MQ=255)                                                                                                  aCAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGAcc                                                 >  1:1103752/1‑151 (MQ=255)                                                                                                         cAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTg                                          >  1:724321/1‑151 (MQ=255)                                                                                                         cAACTGCTGGGGCAGAGTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGTCCACCg                                             <  1:865219/148‑1 (MQ=255)                                                                                                             tgctgGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCG‑‑‑ACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTgg                                    >  2:322052/1‑150 (MQ=255)                                                                                                              gctgGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCt                                      <  1:742777/150‑1 (MQ=255)                                                                                                               ctgGGGCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTgg                                    >  1:576701/1‑151 (MQ=255)                                                                                                                 ggggCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGaa                                  <  2:936938/151‑1 (MQ=255)                                                                                                                 ggggCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGa                                   <  1:922279/150‑1 (MQ=255)                                                                                                                  gggCAGATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAg                                 <  2:624515/151‑1 (MQ=255)                                                                                                                       gATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTcc                            <  1:1081923/151‑1 (MQ=255)                                                                                                                       gATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTcc                            <  1:228868/151‑1 (MQ=255)                                                                                                                       gATTGTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTcc                            <  2:1283998/151‑1 (MQ=255)                                                                                                                           gTACTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAAc                         <  2:517500/150‑1 (MQ=255)                                                                                                                           gTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACg                        >  2:57037/1‑151 (MQ=255)                                                                                                                             aCTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTc                      >  2:447391/1‑151 (MQ=255)                                                                                                                             aCTGGCGAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTc                      >  2:768035/1‑151 (MQ=255)                                                                                                                                   gAACTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGAtctc                >  2:803779/1‑151 (MQ=255)                                                                                                                                      cTTTGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGt             >  1:873378/1‑151 (MQ=255)                                                                                                                                        ttGCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTa            >  1:1035356/1‑150 (MQ=255)                                                                                                                                         tGCCAACAACGACGGTCTGGCATCCGAAGGCTACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTa            <  2:619435/149‑1 (MQ=255)                                                                                                                                          gCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTaaa          <  2:1213464/150‑1 (MQ=255)                                                                                                                                          gCCAACAACGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTaa           >  2:936312/1‑149 (MQ=255)                                                                                                                                             aacaacGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAAc      >  1:303243/1‑151 (MQ=255)                                                                                                                                              acaacGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACt     <  1:1110164/151‑1 (MQ=255)                                                                                                                                               caacGACGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTg    <  1:456847/151‑1 (MQ=255)                                                                                                                                                 acgacgGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGt  <  2:386643/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGT  >  NC_000913/1133394‑1133688

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