Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,019,349	G→A	6.1%	G322S (GGC→AGC)	fliK →	flagellar hook‑length control protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,019,349	0	G	A	6.1%	226.4 / 3.8	66	G322S (GGC→AGC)	fliK	flagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base A (2/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCA  >  NC_000913/2019212‑2019487                                                                                                                                          |                                                                                                                                           aaGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGc                                                                                                                                          <  2:408978/140‑1 (MQ=255) aaGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGc                                                                                                                                          <  2:10685/140‑1 (MQ=255)  aGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCg                                                                                                                                         <  2:248196/140‑1 (MQ=255)    tGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGaa                                                                                                                                       <  1:661456/140‑1 (MQ=255)    tGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGaa                                                                                                                                       <  2:647285/140‑1 (MQ=255)    tGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGaa                                                                                                                                       <  1:813075/140‑1 (MQ=255)    tGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGaa                                                                                                                                       <  2:296026/140‑1 (MQ=255)    tGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGGGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGaa                                                                                                                                       <  2:679241/140‑1 (MQ=255)      gatgatAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAg                                                                                                                                     >  1:598964/1‑140 (MQ=255)      gatgatAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAg                                                                                                                                     >  2:854890/1‑140 (MQ=255)      gatgatAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAg                                                                                                                                     >  2:1083812/1‑140 (MQ=255)                aGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGt                                                                                                                           <  1:138098/140‑1 (MQ=255)                aGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGt                                                                                                                           <  1:181184/140‑1 (MQ=255)                 ggCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTc                                                                                                                          >  1:101708/1‑140 (MQ=255)                 ggCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTc                                                                                                                          >  1:634044/1‑140 (MQ=255)                 ggCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTc                                                                                                                          >  1:1090958/1‑140 (MQ=255)                 ggCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTc                                                                                                                          >  1:1090911/1‑140 (MQ=255)                           cAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGc                                                                                                                <  1:226685/140‑1 (MQ=255)                                     cACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCcag                                                                                                      >  1:727895/1‑140 (MQ=255)                                     cACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCcag                                                                                                      >  1:727853/1‑140 (MQ=255)                                            tCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCTGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCTAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCcagcaacagc                                                                                               >  1:784462/1‑140 (MQ=255)                                            tCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCcagcaacagc                                                                                               >  1:262929/1‑140 (MQ=255)                                            tCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCcagcaacagc                                                                                               >  1:262812/1‑140 (MQ=255)                                            tCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCcagcaacagc                                                                                               <  1:219760/140‑1 (MQ=255)                                                  tGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGcc                                                                                         >  1:435112/1‑140 (MQ=255)                                                     aCGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTAGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAAc                                                                                      >  1:582932/1‑140 (MQ=255)                                                     aCGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTAGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAAc                                                                                      >  1:583816/1‑140 (MQ=255)                                                             cccTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCa                                                                              <  2:232238/140‑1 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGTAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAAc                                                                    >  1:941772/1‑140 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAAc                                                                    >  2:457202/1‑140 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAAc                                                                    >  2:458107/1‑140 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAAc                                                                    >  2:329957/1‑140 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAAc                                                                    >  1:464489/1‑140 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAAc                                                                    >  2:458147/1‑140 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAAc                                                                    >  2:1095904/1‑140 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTGAGGGCAAAGCAACAGCAGTGGCGAAAGCTGTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGacc                                                                    >  2:552756/1‑138 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGGTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCTAAGCCAACGCACAGCAAACCATGAAc                                                                    >  2:594465/1‑140 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATGCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTGTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGacc                                                                    >  2:552893/1‑138 (MQ=255)                                                                       agcGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCAGTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAAc                                                                    >  2:325161/1‑140 (MQ=255)                                                                        gcgcTGCCGGTACTGCGCACGCAGCTGGCCGAAAGAGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGAAATCCATGAAcc                                                                   >  2:648780/1‑140 (MQ=255)                                                                         cgcTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACct                                                                  >  1:1087423/1‑140 (MQ=255)                                                                                   aCTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCggggg                                                        >  1:907534/1‑140 (MQ=255)                                                                                     tGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGaa                                                      <  1:465205/140‑1 (MQ=255)                                                                                      gcgcACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAg                                                     >  2:188669/1‑140 (MQ=255)                                                                                        gcacgcaGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAgac                                                   >  1:1226169/1‑140 (MQ=255)                                                                                         cacgcaGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAgacg                                                  >  2:706702/1‑140 (MQ=255)                                                                                                   gCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCt                                        <  2:1226169/140‑1 (MQ=255)                                                                                                       aaaGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCg                                    <  2:860694/140‑1 (MQ=255)                                                                                                       aaaGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCg                                    <  2:860498/140‑1 (MQ=255)                                                                                                       aaaGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCg                                    <  1:202638/140‑1 (MQ=255)                                                                                                        aaGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCgg                                   <  1:1129522/140‑1 (MQ=255)                                                                                                            ggCATTGAGTTCGGGCAATGCAACATCAGTGGCGTAAGCTTTTGTGGTCTGCAGCAGGCCGCTTCCCAGCATCAGCTAAGCCAAAGCACAGCAAACCATGACCCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTcc                               <  2:472497/140‑1 (MQ=255)                                                                                                              cATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAACAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTccca                             >  1:774436/1‑139 (MQ=255)                                                                                                              cATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAACAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTccca                             >  1:774501/1‑139 (MQ=255)                                                                                                               aTTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTcccat                            <  2:150356/140‑3 (MQ=255)                                                                                                               aTTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTcccat                            <  2:150156/140‑3 (MQ=255)                                                                                                               aTTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTcccat                            <  2:539228/140‑3 (MQ=255)                                                                                                               aTTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTcccat                            <  2:408353/140‑3 (MQ=255)                                                                                                               aTTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTcccat                            <  1:366455/140‑3 (MQ=255)                                                                                                                 tCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCAtct                          <  1:445322/140‑1 (MQ=255)                                                                                                                 tCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCAtct                          <  2:1089339/140‑1 (MQ=255)                                                                                                                       aGGGCAAAGCAACATCAGTAGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTAc                    <  2:582932/140‑1 (MQ=255)                                                                                                                       aGGGCAAAGCAACATCAGTAGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTAc                    <  2:583816/140‑1 (MQ=255)                                                                                                                             aaGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGc              >  2:477249/1‑140 (MQ=255)                                                                                                                                   cATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTaa        <  1:184175/140‑1 (MQ=255)                                                                                                                                         tGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCa  >  1:651059/1‑140 (MQ=255)                                                                                                                                          ggCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTc                                      <  2:41095/103‑1 (MQ=255)                                                                                                                                          ggCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTc                                      >  1:41095/1‑103 (MQ=255)                                                                                                                                          ggCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTc                                      >  1:41172/1‑103 (MQ=255)                                                                                                                                          ggCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTc                                      <  2:41172/103‑1 (MQ=255)                                                                                                                                          |                                                                                                                                           AAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCA  >  NC_000913/2019212‑2019487

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