Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,133,069	A→G	T166A (ACT→GCT)	flgE →	flagellar hook protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,133,069	0	A	G	99.0%	379.9 / NA	97	T166A (ACT→GCT)	flgE	flagellar hook protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (51/45);  minor base T (0/1);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGGGTTACCCGGCAACCGGTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGAT  >  NC_000913/1132921‑1133215                                                                                                                                                     |                                                                                                                                                   gACGGGTTACCCGGCAACCGGTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCt                                                                                                                                                  >  4:225285/1‑151 (MQ=255) gACGGGTTACCCGGCAACCGGTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCt                                                                                                                                                  >  4:170712/1‑151 (MQ=255)     ggTTACCCGGCAACCGGTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGt                                                                                                                                                >  1:1126324/1‑149 (MQ=255)             ggCAACCGGTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTc                                                                                                                                      <  4:9443/151‑1 (MQ=255)              gCAACCGGTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCa                                                                                                                                     <  1:89135/151‑1 (MQ=255)              gCAACCGGTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCa                                                                                                                                     <  1:50789/151‑1 (MQ=255)               cAACCGGTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTc                                                                                                                                      >  2:352655/1‑149 (MQ=255)                    ggTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGcc                                                                                                                                >  1:779081/1‑150 (MQ=255)                    ggTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCa                                                                                                                               <  3:399602/151‑1 (MQ=255)                      tACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGc                                                                                                                             >  1:300259/1‑151 (MQ=255)                       aCGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGc                                                                                                                             >  1:1070565/1‑150 (MQ=255)                        cgccgcCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCa                                                                                                                            >  4:164916/1‑150 (MQ=255)                         gccgccGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAAt                                                                                                                          >  1:64965/1‑151 (MQ=255)                          ccgccgACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATg                                                                                                                         >  4:127230/1‑151 (MQ=255)                           cgccgACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGc                                                                                                                        <  2:106604/151‑1 (MQ=255)                            gccgACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCg                                                                                                                       >  2:251268/1‑151 (MQ=255)                              cgACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGa                                                                                                                     <  3:191705/151‑1 (MQ=255)                               gACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGAt                                                                                                                    <  3:88262/151‑1 (MQ=255)                                aCTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATa                                                                                                                   <  2:32049/151‑1 (MQ=255)                                 cTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATa                                                                                                                   <  2:1227690/150‑1 (MQ=255)                                   aTTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTAAGCGCCAGCAATGCGGATAGCt                                                                                                                >  2:402529/1‑151 (MQ=255)                                    ttCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCGCCCCGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTTCTGTTTCGCCATTCCGCGCCAGCAATGCTGATa                                                                                                                   <  3:299997/147‑1 (MQ=255)                                     tCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCtat                                                                                                              >  2:193098/1‑151 (MQ=255)                                      cAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCtata                                                                                                             <  1:705139/151‑1 (MQ=255)                                        gCAAGGGGCGAATCCGACCAATATTTCGATCCTGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCACATCAACCTGAATTCCAGTGCTCCGCTTCCTGCTGTTACGCCTTTCAGCGCCAGCAATGCGGATAGCtata                                                                                                             >  2:1200150/1‑149 (MQ=255)                                        gCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAAc                                                                                                           >  4:242401/1‑151 (MQ=255)                                         cAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACa                                                                                                          <  1:358090/151‑1 (MQ=255)                                            ggggCGAGTCCGACCAAGATTACGATCCCGAATACCCTGATGGCAGCGAAAACTACCGCCACGGCATCGATGCAGATCACCCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACaaa                                                                                                        <  2:483587/150‑1 (MQ=255)                                             gggCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACaaaa                                                                                                       >  1:255925/1‑150 (MQ=255)                                                 gAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAg                                                                                                    >  1:370754/1‑149 (MQ=255)                                                 gAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGt                                                                                                  <  1:88790/151‑1 (MQ=255)                                                  aaTCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGtt                                                                                                 <  1:858864/151‑1 (MQ=255)                                                      cGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATTAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGt                                                                                             >  1:202396/1‑151 (MQ=255)                                                        aCCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGa                                                                                           <  1:748588/151‑1 (MQ=255)                                                        aCCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGa                                                                                           <  3:373177/151‑1 (MQ=255)                                                            atatTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGt                                                                                       <  3:117068/151‑1 (MQ=255)                                                              atTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGttt                                                                                     >  2:622284/1‑151 (MQ=255)                                                                ttCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTc                                                                                   >  1:548166/1‑151 (MQ=255)                                                                ttCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTc                                                                                   >  1:135401/1‑151 (MQ=255)                                                                  cGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCg                                                                                  <  1:965601/150‑1 (MQ=255)                                                                  cGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGa                                                                                 <  2:490555/151‑1 (MQ=255)                                                                   gATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGAc                                                                                >  2:36675/1‑151 (MQ=255)                                                                     tCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGAc                                                                                >  2:393887/1‑149 (MQ=255)                                                                        cGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTc                                                                            >  2:941532/1‑150 (MQ=255)                                                                         gAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAg                                                                          <  3:205008/151‑1 (MQ=255)                                                                          aaTACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGtttt                                                                                    <  1:1068276/140‑1 (MQ=255)                                                                          aaTACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGtttt                                                                                    >  2:1068076/1‑140 (MQ=255)                                                                             aCCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGt                                                                       >  1:83915/1‑150 (MQ=255)                                                                              cccTGATGGCAGCGAAAACCACCCCCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTa                                                                      <  1:298895/150‑1 (MQ=255)                                                                                  gATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATg                                                                   <  1:1172618/149‑1 (MQ=255)                                                                                  gATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGc                                                                  <  4:206956/150‑1 (MQ=255)                                                                                  gATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCt                                                                 >  2:1172408/1‑151 (MQ=255)                                                                                    tGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTc                                                                >  2:213949/1‑150 (MQ=255)                                                                                       cAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTcatg                                                             >  2:1039673/1‑150 (MQ=255)                                                                                           gAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTcatgacatg                                                        >  2:261120/1‑151 (MQ=255)                                                                                           gAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTcatgacatg                                                        >  2:840170/1‑151 (MQ=255)                                                                                                  accaccACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTc                                                  <  1:943849/150‑1 (MQ=255)                                                                                                  accaccACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCt                                                 >  1:593286/1‑151 (MQ=255)                                                                                                   ccaccaCGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCt                                                 >  1:1155787/1‑150 (MQ=255)                                                                                                    cacctcGACATCGATGCAGATCACCCTCAATTCCATTGATCCGCTTCCTGCTGATACGCCATACAGCGCCAGCAATGCAGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGAGTAATGCTCATGGCATGAGCTTCTAc                                               <  2:321975/151‑1 (MQ=255)                                                                                                      ccacGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACtt                                             >  4:528/1‑151 (MQ=255)                                                                                                       cacGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACttt                                            <  2:549896/151‑1 (MQ=255)                                                                                                        acGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTtg                                           <  4:38312/151‑1 (MQ=255)                                                                                                             aTCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGa                                        >  2:647303/1‑149 (MQ=255)                                                                                                               cGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGa                                     <  3:85052/150‑1 (MQ=255)                                                                                                               cGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGAc                                    >  2:618305/1‑151 (MQ=255)                                                                                                                gATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGAcc                                   >  1:300217/1‑151 (MQ=255)                                                                                                                gATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGAcc                                   <  1:889120/151‑1 (MQ=255)                                                                                                                   gCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCggg                                <  2:391612/151‑1 (MQ=255)                                                                                                                   gCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCgg                                 <  3:528/150‑1 (MQ=255)                                                                                                                   gCAGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCgg                                 <  4:196635/150‑1 (MQ=255)                                                                                                                     acaTCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCgggg                               >  2:196774/3‑150 (MQ=255)                                                                                                                     aGATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATGACAATAAATGTTCGGTGACTGATTTCGACAGTCCGGGGAATGCTCATGACATGAGTGTCTACTTTGTGAAGACCgggg                               >  2:771711/1‑150 (MQ=255)                                                                                                                      gATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGat                             >  2:701868/1‑151 (MQ=255)                                                                                                                      gATCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGat                             <  4:159541/151‑1 (MQ=255)                                                                                                                       aTCAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGata                            <  2:924147/151‑1 (MQ=255)                                                                                                                         cAACCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGataat                          >  1:1218804/1‑151 (MQ=255)                                                                                                                           aCCTGAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGAtaataa                        <  1:267191/151‑1 (MQ=255)                                                                                                                               gAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTTAGGGTAATGCTCATGACATGGGCGTCTTCTTTGTGAAGACCGGGGATAATAACTgg                    <  1:128109/151‑1 (MQ=255)                                                                                                                               gAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTgg                    <  1:847463/151‑1 (MQ=255)                                                                                                                               gAATTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTg                     <  4:50705/150‑1 (MQ=255)                                                                                                                                 aTTCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCa                  <  2:1242252/151‑1 (MQ=255)                                                                                                                                   tCCAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAgg                >  1:1094755/1‑151 (MQ=255)                                                                                                                                     cAGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTc              >  2:170402/1‑151 (MQ=255)                                                                                                                                      aGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCt             <  1:251305/151‑1 (MQ=255)                                                                                                                                      aGTGATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCt             >  4:444345/1‑151 (MQ=255)                                                                                                                                         gATCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTac           >  2:413478/1‑150 (MQ=255)                                                                                                                                          aTCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTacac         <  2:1026553/151‑1 (MQ=255)                                                                                                                                          aTCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTacac         >  2:171049/1‑151 (MQ=255)                                                                                                                                           tCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACAcc        >  3:280117/1‑151 (MQ=255)                                                                                                                                           tCCGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACAcc        >  2:1247307/1‑151 (MQ=255)                                                                                                                                            ccGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACAccc       <  1:770649/151‑1 (MQ=255)                                                                                                                                            ccGCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACAccc       >  2:975549/1‑151 (MQ=255)                                                                                                                                              gCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACAccc       >  1:803767/1‑149 (MQ=255)                                                                                                                                              gCTTCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAg     <  4:406018/151‑1 (MQ=255)                                                                                                                                                 tCCTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGAt  >  2:673565/1‑151 (MQ=255)                                                                                                                                                  ccTGCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGAt  <  2:1158627/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GACGGGTTACCCGGCAACCGGTACGCCGCCGACTATTCAGCAAGGGGCGAATCCGACCAATATTTCGATCCCGAATACCCTGATGGCAGCGAAAACTACCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGAT  >  NC_000913/1132921‑1133215

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