Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,946,790	T→G	G180G (GGT→GGG)	yebB →	DUF830 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,946,790	0	T	G	94.6%	321.1 / 3.8	92	G180G (GGT→GGG)	yebB	DUF830 family protein
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base G (38/49);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCC  >  NC_000913/1946695‑1946886                                                                                                |                                                                                                 aaTGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTgg                                                                                              <  2:719171/100‑1 (MQ=255) aaTGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATGCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTg                                                                                               <  1:1581203/99‑1 (MQ=255)     cAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGtt                                                                                          >  2:651865/1‑100 (MQ=255)     cAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGtt                                                                                          >  1:161766/1‑100 (MQ=255)     cAAAACTAACTGTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGtt                                                                                          >  1:1020160/1‑100 (MQ=255)       aaaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGa                                                                                        <  2:1707473/100‑1 (MQ=255)        aaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGAt                                                                                       <  2:1567756/100‑1 (MQ=255)        aaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGAt                                                                                       <  2:1041971/100‑1 (MQ=255)        aaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGAt                                                                                       <  2:66756/100‑1 (MQ=255)          cTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTc                                                                                     <  1:434366/100‑1 (MQ=255)            cACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCAt                                                                                   >  1:1825066/1‑100 (MQ=255)              cTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGc                                                                                 <  1:1295864/100‑1 (MQ=255)              cTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGc                                                                                 <  2:147267/100‑1 (MQ=255)                  cTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggag                                                                             <  1:1539697/100‑1 (MQ=255)                  cTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggag                                                                             <  1:914054/100‑1 (MQ=255)                  cTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggag                                                                             >  1:1676088/1‑100 (MQ=255)                   tGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggagg                                                                            >  1:995697/1‑100 (MQ=255)                   tGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggagg                                                                            >  2:105996/1‑100 (MQ=255)                   tGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggagg                                                                            >  1:1824651/1‑100 (MQ=255)                    ggAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggaggg                                                                           <  1:936571/100‑1 (MQ=255)                      aaaTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggagggag                                                                         <  2:1691592/100‑1 (MQ=255)                       aaTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGt                                                                        <  2:293536/100‑1 (MQ=255)                        aTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCCCGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTg                                                                       <  1:1130273/100‑1 (MQ=255)                         ttCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTgg                                                                      >  2:977419/1‑100 (MQ=255)                         ttCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTgg                                                                      <  2:1706262/100‑1 (MQ=255)                         ttCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTgg                                                                      <  1:880826/100‑1 (MQ=255)                          tCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGCGGCATCATCCTGGGTTGGTGTTGATTCATGCGGAGGGAGttg                                                                      <  2:1311150/99‑3 (MQ=255)                           cTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGaa                                                                    >  2:1735016/1‑100 (MQ=255)                           cTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGaa                                                                    <  2:486236/100‑1 (MQ=255)                                tCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGcc                                                               <  1:651864/100‑1 (MQ=255)                                  cTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTc                                                             >  1:834202/1‑100 (MQ=255)                                   tGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCa                                                            <  1:1232062/100‑1 (MQ=255)                                     ggTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGc                                                          >  1:1549103/1‑100 (MQ=255)                                      gTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGcc                                                         >  1:1358952/1‑100 (MQ=255)                                      gTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGGTTCATGGGGAGGGAGTGGAAACGCCTCAGcc                                                         >  2:552216/1‑100 (MQ=255)                                         cTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGa                                                      >  1:1925485/1‑100 (MQ=255)                                          tATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGaa                                                     <  1:110339/100‑1 (MQ=255)                                          tATTCCGTGGGAGCGTAAAAACGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGaa                                                     >  1:1190392/1‑100 (MQ=255)                                           aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                    <  2:493624/100‑1 (MQ=255)                                           aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                    >  1:927758/1‑100 (MQ=255)                                            ttCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACt                                                   <  2:1337455/100‑1 (MQ=255)                                            ttCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACt                                                   <  2:1740660/100‑1 (MQ=255)                                             tCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTg                                                  >  1:172789/1‑100 (MQ=255)                                               cGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGAc                                                >  1:1751642/1‑100 (MQ=255)                                               cGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGAc                                                <  2:707759/100‑1 (MQ=255)                                                   ggAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAg                                            >  1:699887/1‑100 (MQ=255)                                                      gCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCg                                         <  1:742897/100‑1 (MQ=255)                                                       cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                        <  2:67636/100‑1 (MQ=255)                                                       cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                        <  2:428939/100‑1 (MQ=255)                                                       cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                        >  1:1719016/1‑100 (MQ=255)                                                        gTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGGGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                        <  1:289771/99‑1 (MQ=255)                                                        gTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                        >  1:583727/1‑99 (MQ=255)                                                        gTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGt                                       >  1:584400/1‑100 (MQ=255)                                                          aaaaCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGtat                                     <  2:995698/100‑1 (MQ=255)                                                          aaaaCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGtat                                     >  2:1071612/1‑100 (MQ=255)                                                           aaaCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGtata                                    >  2:1763081/1‑100 (MQ=255)                                                           aaaCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGtata                                    >  1:272922/1‑100 (MQ=255)                                                           aaaCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGtata                                    <  1:1883589/100‑1 (MQ=255)                                                           aaaCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGtata                                    >  1:674619/1‑100 (MQ=255)                                                             aCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAAc                                  <  2:1826638/100‑1 (MQ=255)                                                              ccGTCACGCCAGCCTGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGTAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACt                                 <  2:890073/100‑1 (MQ=255)                                                                gTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTa                               >  2:1441298/1‑100 (MQ=255)                                                                  cACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAAc                             <  2:526387/100‑1 (MQ=255)                                                                   atgccagccaGTTTGTGGTATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACg                            <  1:1810333/98‑1 (MQ=255)                                                                   aCGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACg                            <  2:272924/100‑1 (MQ=255)                                                                     gccagccaGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCa                          <  2:1358960/100‑1 (MQ=255)                                                                      ccagccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATCACTTAACGCAg                         >  2:1705831/1‑100 (MQ=255)                                                                      ccagccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAg                         >  1:456371/1‑100 (MQ=255)                                                                       cagccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGt                        >  2:1441064/1‑100 (MQ=255)                                                                       cagccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGt                        >  2:382324/1‑100 (MQ=255)                                                                        agccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTc                       <  2:1925498/100‑1 (MQ=255)                                                                          ccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGc                     <  2:529426/100‑1 (MQ=255)                                                                           cagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGcc                    <  2:1296271/100‑1 (MQ=255)                                                                           cagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGcc                    <  2:1901017/100‑1 (MQ=255)                                                                               ttGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCctct                <  2:456373/100‑1 (MQ=255)                                                                                tgtgGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCtctc               <  1:1838137/100‑1 (MQ=255)                                                                                tgtgGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCtctc               >  1:919205/1‑100 (MQ=255)                                                                                 gtgGCATCATCCGGTGTTGGTGTTGATTCATGCGGAGAGAGTGGAAACGCCTCAGTCGGAACTTTCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCg              <  2:1225867/100‑1 (MQ=255)                                                                                 gtgGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCg              >  1:821691/1‑100 (MQ=255)                                                                                 gtgGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCg              >  2:1352876/1‑100 (MQ=255)                                                                                    gCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCa           >  2:283526/1‑100 (MQ=255)                                                                                    gCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCa           >  1:483940/1‑100 (MQ=255)                                                                                    gCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCa           <  1:1390937/100‑1 (MQ=255)                                                                                    gCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCa           >  1:856192/1‑100 (MQ=255)                                                                                     catcatCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAg          >  2:1442566/1‑100 (MQ=255)                                                                                      atcatcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAgg         <  2:1549113/100‑1 (MQ=255)                                                                                        catcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGtt       <  1:977418/100‑1 (MQ=255)                                                                                         atcCGGGTTTGGTGTTGATTAATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc      >  1:1187617/1‑100 (MQ=255)                                                                                         atcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc      <  1:1667107/100‑1 (MQ=255)                                                                                           ccGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAg    >  1:746604/1‑100 (MQ=255)                                                                                           ccGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAg    >  2:1020168/1‑100 (MQ=255)                                                                                            cGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGc   <  2:1824662/100‑1 (MQ=255)                                                                                             ggggTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGcc  >  1:1718231/1‑100 (MQ=255)                                                                                                |                                                                                                 AATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCC  >  NC_000913/1946695‑1946886

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