Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	207,483	A→C	17.9%	V786V (GTA→GTC)	dnaE →	DNA polymerase III alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	207,483	0	A	C	17.9%	217.1 / 43.0	99	V786V (GTA→GTC)	dnaE	DNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (34/47);  new base C (6/12);  total (40/59)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGG  >  NC_000913/207385‑207581                                                                                                   |                                                                                                   tACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTaa                                                                                                   >  1:966253/1‑100 (MQ=255)  aCGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAAt                                                                                                  >  1:1620188/1‑100 (MQ=255)   cGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATg                                                                                                 <  2:1334948/100‑1 (MQ=255)   cGGATTTAACAAATCGCACTATGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATg                                                                                                 >  1:689214/1‑100 (MQ=255)     gATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAc                                                                                               >  1:52113/1‑100 (MQ=255)           aCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgc                                                                                         <  1:1269532/100‑2 (MQ=255)           aCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgc                                                                                         <  2:273776/100‑2 (MQ=255)           aCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGa                                                                                         <  1:249982/100‑1 (MQ=255)            cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                        <  2:702688/100‑3 (MQ=255)            cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                        <  2:27561/100‑3 (MQ=255)            cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGat                                                                                        <  1:846030/100‑1 (MQ=255)             aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcca                                                                                       <  2:2290482/100‑4 (MQ=255)             aaaTCGCACTCTGCGGCCTATGCGTTGGTGTCATATCAAACGTTATGGCTGAAAACGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAccgccctta                                                                                       >  1:1059327/1‑96 (MQ=255)              aaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTGATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcca                                                                                       <  1:1654686/99‑4 (MQ=255)               aTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccata                                                                                     <  2:271146/100‑6 (MQ=255)               aTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATg                                                                                     <  1:395014/100‑1 (MQ=255)               aTCGCACTCTGCCGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                      <  2:2322146/99‑5 (MQ=255)                tCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATgg                                                                                    <  2:1296084/100‑1 (MQ=255)                tCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATgg                                                                                    <  1:1418544/100‑1 (MQ=255)                tCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATgg                                                                                    <  2:56049/100‑1 (MQ=255)                tCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATg                                                                                     <  1:8736/99‑1 (MQ=255)                 cGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGa                                                                                   <  1:2320042/100‑1 (MQ=255)                 cGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCGGAAAGCTCACTATCCTGCGGAGTTTATGGCACCGGTAATGACCGCCGATATgg                                                                                    >  1:379894/1‑99 (MQ=255)                 cGCACCATGCGGCTTATGCATTGCTTTCACATCAACCGTTATGGGTGAAAGCGCACTATCCTGCGGAGTCTATGGCGGCGGTCATTAccgccgcc                                                                                        <  1:2262290/95‑3 (MQ=255)                  gCACTCTTCGGCCTATGCTTTGGTGTCATATCAAACGTTACGGCTGAAAGCGCACTACCCTGACGAGTTCATGGCGGCGGTCATTGccgccgccataaac                                                                                  <  1:2418780/100‑9 (MQ=255)                  gCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGac                                                                                  <  1:2201117/100‑1 (MQ=255)                  gCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGac                                                                                  >  1:1275745/1‑100 (MQ=255)                  gCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGac                                                                                  >  1:1403258/1‑100 (MQ=255)                  gCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGac                                                                                  <  1:2422452/100‑1 (MQ=255)                   cACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                 <  1:2019128/100‑1 (MQ=255)                    aCTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccataaactc                                                                                <  2:218151/100‑11 (MQ=38)                     gtctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTACGGAGTTTATCGCGGCGGTAACAACCGCCGAGATGGaca                                                                                 >  1:1113831/2‑98 (MQ=255)                     gtctGCGGCCTATAATTTGGTGCCATATCAAACGTTATGGCTGCAAGCGCACTATACTGCGGAGTTTAAGGCGGCAGAAAAGACTGCCCAAGTGGac                                                                                  >  1:1270497/2‑97 (MQ=255)                     ctctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGacaac                                                                               >  2:349282/1‑100 (MQ=255)                       ctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAacac                                                                             >  2:2086800/1‑100 (MQ=255)                       ctGCGGCCTATGCTTTGGTGTCATATCAAACGCTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAacac                                                                             >  2:2461265/1‑100 (MQ=255)                         gCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCg                                                                           >  2:234569/1‑100 (MQ=255)                             ccTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGaa                                                                       <  2:276070/100‑1 (MQ=255)                                 tGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggtg                                                                   <  2:1485333/100‑1 (MQ=255)                                   cTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggtggt                                                                 >  2:31414/1‑100 (MQ=255)                                   cTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTACGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggtg                                                                   >  1:554187/1‑98 (MQ=255)                                    tttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGgtggtg                                                                <  1:1216083/100‑1 (MQ=255)                                     ttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGtggtgg                                                               <  2:1838118/100‑1 (MQ=255)                                      tGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTggg                                                              >  2:1606171/1‑100 (MQ=255)                                       ggTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGt                                                             >  2:1718183/1‑100 (MQ=255)                                          gtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtg                                                          >  1:13807/1‑100 (MQ=255)                                          gtCATATCAAACGATATGGCTGAAAGCGCACTGTCCTGCGGAGCTTATGGCGGCGGTAATGACCGCCGCTATGGACAACACCGAGAAGGCGGTGGGACtg                                                          >  1:1168972/1‑100 (MQ=255)                                            cATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggt                                                        >  2:1001181/1‑100 (MQ=255)                                             atatCAATCGTCATGGCTCAAAGCCCACCACCATGAGGAGTGTAAAGCGTCGGTCATGACCGCCGATATGCACAACACCGAGAAGGTGGTGGGTCtggt                                                        <  1:2211240/99‑1 (MQ=255)                                             atatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtg                                                       <  1:746124/100‑1 (MQ=255)                                             atatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtg                                                       <  2:1385454/100‑1 (MQ=255)                                              tatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                      <  1:920408/100‑1 (MQ=255)                                              tatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                      <  1:2348917/100‑1 (MQ=255)                                               atCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTATTGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGa                                                     <  1:1126593/100‑1 (MQ=255)                                               atCAAACGTTATGGCTGAAAGCGCACTAGCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGa                                                     <  1:2461248/100‑1 (MQ=255)                                                 cAAACGTTATGGCTGAAAGCGCATTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATg                                                   >  2:111206/1‑100 (MQ=255)                                                 cAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATg                                                   <  1:2568022/100‑1 (MQ=255)                                                  aaaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGa                                                  >  2:252094/1‑100 (MQ=255)                                                  aaaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGa                                                  <  1:2052289/100‑1 (MQ=255)                                                  aaaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGa                                                  >  2:2010320/1‑100 (MQ=255)                                                   aaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAg                                                 <  2:1341065/100‑1 (MQ=255)                                                     cGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTg                                               <  2:588059/100‑1 (MQ=255)                                                         aTGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTgg                                           <  1:1346714/100‑1 (MQ=255)                                                           ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCg                                         >  2:2518106/1‑100 (MQ=255)                                                           ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCg                                         <  2:265814/100‑1 (MQ=255)                                                           ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCg                                         <  1:1718174/100‑1 (MQ=255)                                                            gCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCgg                                        >  2:2477118/1‑100 (MQ=255)                                                            gCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCgg                                        >  2:670752/1‑100 (MQ=255)                                                              tGAAAGCGCACCATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGGAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGAt                                      >  1:1359458/1‑100 (MQ=255)                                                                 aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATggg                                   <  2:92520/100‑1 (MQ=255)                                                                   gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGc                                 <  1:1606162/100‑1 (MQ=255)                                                                   gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGc                                 <  2:1162453/100‑1 (MQ=255)                                                                   gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGc                                 <  1:897277/100‑1 (MQ=255)                                                                    cgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGc                                 >  1:651318/1‑99 (MQ=255)                                                                    cgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCt                                >  2:676544/1‑100 (MQ=255)                                                                       aCTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGaa                             >  1:1949937/1‑100 (MQ=255)                                                                          aTCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAAt                          >  2:1528033/1‑100 (MQ=255)                                                                            ccTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATcc                        <  1:1143924/100‑1 (MQ=255)                                                                            ccTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATcc                        <  2:52114/100‑1 (MQ=255)                                                                              tGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTg                      <  2:1910719/100‑1 (MQ=255)                                                                              tGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTg                      <  1:1009904/100‑1 (MQ=255)                                                                              tGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTg                      <  2:2427672/100‑1 (MQ=255)                                                                               gCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGc                     <  2:616318/100‑1 (MQ=255)                                                                               gCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGc                     <  2:345039/100‑1 (MQ=255)                                                                                 ggAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGcc                    <  1:234313/99‑1 (MQ=255)                                                                                   aGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacc                 <  1:844237/100‑1 (MQ=255)                                                                                   aGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacc                 <  1:346207/100‑1 (MQ=255)                                                                                      ttATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAg               >  1:2591195/1‑99 (MQ=255)                                                                                       catatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGat             >  2:2001894/6‑100 (MQ=255)                                                                                       catatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGat             >  2:113102/6‑100 (MQ=255)                                                                                       catatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGCCTGGCGGATGAGTGCAGCCGGACGGGGCTGAAAATCCTGccacca                >  1:1281517/6‑97 (MQ=255)                                                                                        aTGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGata            <  2:379800/100‑1 (MQ=255)                                                                                          ggcggcGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGAtata          <  2:1017982/100‑1 (MQ=255)                                                                                          atcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGAtata          >  2:1460538/3‑100 (MQ=255)                                                                                          atcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGAtata          >  2:2200139/3‑100 (MQ=255)                                                                                           tcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGACTGGTGAATGAGTGCTGGCGCATGGGGCTGAAAATCCTGCCACCAGATATaa         >  1:1858995/2‑100 (MQ=255)                                                                                             ggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAAc       <  1:111205/100‑1 (MQ=255)                                                                                                 gTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCg   >  1:58219/1‑100 (MQ=255)                                                                                                 gTAATGACCGCCGATATGGACAACACCGAGAAGGAGGTGGGTCTGGTGGATGAGAGCTGGCGGATGGGGCTGAAAATCCTGCCCCCAGATATAAACTCCg   >  1:849538/1‑100 (MQ=255)                                                                                                  tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTAGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg  <  2:1545736/100‑1 (MQ=255)                                                                                                   |                                                                                                   TACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGG  >  NC_000913/207385‑207581

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