Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	171,343	G→A	A257T (GCC→ACC)	fhuD →	iron(3+)‑hydroxamate import ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	171,343	0	G	A	94.7%	306.8 / 3.1	95	A257T (GCC→ACC)	fhuD	iron(3+)‑hydroxamate import ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (46/44);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGTCTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGGCCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTGCGCGTTCT  >  NC_000913/171247‑171437                                                                                                 |                                                                                               gatcgTCTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCa                                                                                             >  2:1482226/1‑100 (MQ=255) gatcgTCTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCa                                                                                             >  1:422759/1‑100 (MQ=255)   tcgTCTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATg                                                                                           <  1:125551/100‑1 (MQ=255)   tcgTCTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATg                                                                                           <  2:1471105/100‑1 (MQ=255)   tcgTCTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATg                                                                                           <  2:1471062/100‑1 (MQ=255)     gTCTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGAAAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGcc                                                                                         >  2:1298824/1‑100 (MQ=255)      tCTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCg                                                                                        >  2:751691/1‑100 (MQ=255)      tCTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGAAATGGATGCGCTAATGGCAACGCCGCTGTGGCAGGCAATGCCg                                                                                        >  2:2101162/1‑100 (MQ=255)       cTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGt                                                                                       <  2:430275/100‑1 (MQ=255)          gcggcgTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTg                                                                                    <  2:760971/100‑1 (MQ=255)            ggcgTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTc                                                                                  >  1:1984034/1‑100 (MQ=255)                tataAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCcgcg                                                                              >  1:411435/1‑100 (MQ=255)                tataAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCcgc                                                                               >  2:1612198/1‑99 (MQ=255)                 ataAAGACGTTGCTGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCgcgc                                                                             >  1:1539302/1‑100 (MQ=255)                  taAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGcc                                                                            >  1:1320769/1‑100 (MQ=255)                  taAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGcc                                                                            <  2:1286941/100‑1 (MQ=255)                   aaaGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCg                                                                           <  2:1060197/100‑1 (MQ=255)                     aGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGa                                                                         <  2:422757/100‑1 (MQ=255)                       aCGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACg                                                                       <  1:765002/100‑1 (MQ=255)                        cGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGc                                                                      >  2:1205403/1‑100 (MQ=255)                            gATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTc                                                                  <  2:945888/100‑1 (MQ=255)                              tgtgCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGGCCATGCCGTTTGTCCGCGCCGGACGCTTTCAg                                                                <  1:2101169/100‑1 (MQ=255)                              tgtgCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAg                                                                >  2:181620/1‑100 (MQ=255)                                 gCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGGCCATGCCGTTTGTCCGCGCCGGACGCTTTCAgcgc                                                             >  1:30483/1‑100 (MQ=255)                                 gCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAgcgc                                                             >  1:625775/1‑100 (MQ=255)                                   tctGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGt                                                           >  1:1710459/1‑100 (MQ=255)                                    ctGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTa                                                          >  1:1412999/1‑100 (MQ=255)                                      gTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTAcc                                                        <  1:1298827/100‑1 (MQ=255)                                      gTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTAcc                                                        >  1:1901907/1‑100 (MQ=255)                                       ttttGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCt                                                       >  2:865959/1‑100 (MQ=255)                                        tttGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTg                                                      <  2:1450986/100‑1 (MQ=255)                                        tttGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTg                                                      <  1:1843754/100‑1 (MQ=255)                                         ttGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGc                                                     <  2:1901900/100‑1 (MQ=255)                                         ttGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGc                                                     >  1:995008/1‑100 (MQ=255)                                           gATCACGACAACAGCAAAGACATGGATGCGCTAATGTCAACGCCGCTTTGTCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCa                                                    >  2:1951311/1‑99 (MQ=255)                                           gATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGGCCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCa                                                    >  1:430651/1‑99 (MQ=255)                                           gATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAg                                                   >  2:2192455/1‑100 (MQ=255)                                           gATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAg                                                   <  2:880625/100‑1 (MQ=255)                                           gATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAg                                                   >  2:1961030/1‑100 (MQ=255)                                            aTCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGt                                                  >  1:599671/1‑100 (MQ=255)                                            aTCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGt                                                  >  1:599943/1‑100 (MQ=255)                                             tCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTc                                                 <  2:1707637/100‑1 (MQ=255)                                             tCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTc                                                 <  2:646687/100‑1 (MQ=255)                                               aCGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTg                                               >  1:2187595/1‑100 (MQ=255)                                                cGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTgg                                              >  2:2085741/1‑100 (MQ=255)                                                   caacaGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGttt                                           <  1:1967349/100‑1 (MQ=255)                                                   caacaGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGttt                                           <  2:2270669/100‑1 (MQ=255)                                                   caacaGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGttt                                           <  2:622027/100‑1 (MQ=255)                                                   caacaGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGttt                                           <  2:571804/100‑1 (MQ=255)                                                     acaGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTa                                         <  1:1797073/100‑1 (MQ=255)                                                      caGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTAt                                        <  2:1984027/100‑1 (MQ=255)                                                      caGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTAt                                        >  2:1336091/1‑100 (MQ=255)                                                      caGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTAt                                        <  2:2187588/100‑1 (MQ=255)                                                      caGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTAt                                        >  2:35002/1‑100 (MQ=255)                                                       aGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATg                                       <  1:1482230/100‑1 (MQ=255)                                                            aGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCg                                  <  2:781228/100‑1 (MQ=255)                                                             gACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGa                                 >  2:85319/1‑100 (MQ=255)                                                             gACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGa                                 >  2:13594/1‑100 (MQ=255)                                                              aCATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGAc                                <  2:726980/100‑1 (MQ=255)                                                               cATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACg                               <  2:971806/100‑1 (MQ=255)                                                               cATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACg                               <  2:625774/100‑1 (MQ=255)                                                               cATGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACg                               >  2:376375/1‑100 (MQ=255)                                                                aTGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGc                              <  1:181622/100‑1 (MQ=255)                                                                 tGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGct                             >  1:1097627/1‑100 (MQ=255)                                                                 tGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGct                             >  1:760331/1‑100 (MQ=255)                                                                 tGGATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGc                              <  1:1951318/99‑1 (MQ=255)                                                                  ggATGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGctc                            >  2:675184/1‑100 (MQ=255)                                                                    aTGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCtctc                          <  2:1566839/100‑1 (MQ=255)                                                                     tGCGCTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCg                         >  1:965041/1‑100 (MQ=255)                                                                      gcgcTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCgg                        <  1:865961/100‑1 (MQ=255)                                                                       cgcTAATGGCTACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGTAGTCTGGTTTTATGGTGCGACGCTCTCGGc                       <  2:1539298/100‑1 (MQ=255)                                                                       cgcTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGc                       >  2:2105570/1‑100 (MQ=255)                                                                       cgcTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGc                       >  1:835749/1‑100 (MQ=255)                                                                         cTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCaa                     >  1:1835965/1‑100 (MQ=255)                                                                         cTAATGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCaa                     >  1:1835953/1‑100 (MQ=255)                                                                            aTGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGc                  <  1:35003/100‑1 (MQ=255)                                                                            aTGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGc                  <  2:445078/100‑1 (MQ=255)                                                                             tGGCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCa                 >  1:2060688/1‑100 (MQ=255)                                                                              ggCAACGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCAc                >  2:2142969/1‑100 (MQ=255)                                                                                cAACGCCGCTGTGGCAGGCCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACtt              <  2:30482/100‑1 (MQ=255)                                                                                  aCGCCGCTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTtg            <  2:835747/100‑1 (MQ=255)                                                                                   cgccgcTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTtgt           >  1:636897/1‑100 (MQ=255)                                                                                    tccgatctgGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTtgtg          <  1:2031251/93‑1 (MQ=255)                                                                                     ccgcTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTgc         <  2:1097625/100‑1 (MQ=255)                                                                                      cgcTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTgcg        <  1:1336094/100‑1 (MQ=255)                                                                                      cgcTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTgcg        >  2:178827/1‑100 (MQ=255)                                                                                      cgcTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTgcg        <  1:2192463/100‑1 (MQ=255)                                                                                        cTGTGGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTAATGGTGCGACGCTCTCGGCAATGCACTTTGTGcgcg      >  1:1644292/1‑100 (MQ=255)                                                                                         tgtgGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTGCGCGt     <  2:1835960/100‑1 (MQ=255)                                                                                         tgtgGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTGCGCGt     <  2:1835948/100‑1 (MQ=255)                                                                                          gtgGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTGCGCGtt    >  1:1470107/1‑100 (MQ=255)                                                                                          gtgGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTGGTGCGCGtt    >  1:605868/1‑100 (MQ=255)                                                                                           tgGCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTGCGCGTTc   <  1:178829/100‑1 (MQ=255)                                                                                            ggCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTtgtgagatcgga  >  2:2031244/1‑92 (MQ=255)                                                                                            ggCAGACCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTGCGCGTTCt  >  1:534429/1‑100 (MQ=255)                                                                                                 |                                                                                               GATCGTCTGGCGGCGTATAAAGACGTTGATGTGCTCTGTTTTGATCACGACAACAGCAAAGACATGGATGCGCTAATGGCAACGCCGCTGTGGCAGGCCATGCCGTTTGTCCGCGCCGGACGCTTTCAGCGCGTACCTGCAGTCTGGTTTTATGGTGCGACGCTCTCGGCAATGCACTTTGTGCGCGTTCT  >  NC_000913/171247‑171437

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