Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	216,568	C→T	100%	A146T (GCC→ACC)	yaeF ←	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	216,568	0	C	T	97.8%	290.6 / ‑2.5	91	A146T (GCC→ACC)	yaeF	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (48/41);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACG  >  NC_000913/216473‑216665                                                                                                |                                                                                                  gaACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAc                                                                                               <  2:1731797/100‑1 (MQ=255)   aCGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAcgc                                                                                             >  1:84328/1‑100 (MQ=255)     ggATTTAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCgg                                                                                           >  2:1701387/1‑100 (MQ=255)     gcatgcAGTGCGCACCTCTGACCAGTAACCCGCAAGGGGTTAAATGTGACAATGCCGCGTTAGTTATAACCGCTATCTTTGATTTTATTGGTACACGCgg                                                                                           <  2:1426723/95‑1 (MQ=255)        ttCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTg                                                                                         <  2:1548345/99‑1 (MQ=255)          cAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGata                                                                                      >  2:107739/1‑100 (MQ=255)          cAGTGAGCACATCTGCCGAGTCACCATAAAGGGAATAACTTATACAATGACGCGATAGTAATAACCGATATCTTTGATGTCATTGGTTAACGCggcg                                                                                         >  2:1105466/1‑95 (MQ=255)           aGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGatat                                                                                     >  1:2226227/1‑100 (MQ=255)           aGTGAGCACATCTGGCGAATCATCATAAAGAGAATAAATTCAACAAAACCGCGATAGTTATAACCACTATCTTTGATTTTATTGGTAAACGCgatgat                                                                                       >  2:1389769/1‑98 (MQ=255)            gTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCTCTATCTTTGATTTTATTGGTAAACGCGGTGATATc                                                                                    <  2:650380/100‑1 (MQ=255)            gTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATc                                                                                    >  2:406231/1‑100 (MQ=255)            gTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATc                                                                                    <  1:1529176/100‑1 (MQ=255)             tGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCt                                                                                   <  2:1522238/100‑1 (MQ=255)              gAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtg                                                                                  <  2:1304036/100‑1 (MQ=255)              gAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGAGGTGATATCtg                                                                                  >  2:972805/1‑100 (MQ=255)                gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                                <  2:396315/100‑1 (MQ=255)                   caTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCt                                                                             >  2:28094/1‑100 (MQ=255)                      cTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgt                                                                          <  2:1931320/100‑1 (MQ=255)                        ggCGATTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttg                                                                           >  2:1895486/1‑97 (MQ=255)                        ggCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgttg                                                                        >  1:2239371/1‑100 (MQ=255)                        ggCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGACTTTATTGGTAAACGCGGTGATATCTGTGGCttgttg                                                                        >  1:1118648/1‑100 (MQ=255)                          cGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCg                                                                      >  1:1386897/1‑100 (MQ=255)                             gttaCCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCgggg                                                                   <  1:117402/97‑1 (MQ=255)                             gTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCgggg                                                                   >  1:2302185/1‑100 (MQ=255)                              tCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGt                                                                  >  1:891693/1‑100 (MQ=255)                               cACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTa                                                                 >  1:974512/1‑100 (MQ=255)                               cACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTa                                                                 >  2:916171/1‑100 (MQ=255)                                aCCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTaa                                                                <  2:1866729/100‑1 (MQ=255)                                aCCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCCGTGATATCTGTGGCTTGTTGCGTGGTa                                                                 >  2:1966550/1‑99 (MQ=255)                                  cATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGa                                                              >  2:330779/1‑100 (MQ=255)                                   aTAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGAt                                                             >  1:793695/1‑100 (MQ=255)                                     aaaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATcc                                                           <  1:491620/100‑1 (MQ=255)                                     aaaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATcc                                                           >  1:1134926/1‑100 (MQ=255)                                      aaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCg                                                          >  2:1359374/1‑100 (MQ=255)                                      aaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCg                                                          >  1:1029845/1‑100 (MQ=255)                                       aGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCgg                                                         <  1:1409472/100‑1 (MQ=255)                                       aGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGAATTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCgg                                                         >  2:481227/1‑100 (MQ=255)                                        gggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCggg                                                        <  1:2121036/100‑1 (MQ=255)                                        gggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCggg                                                        >  1:351016/1‑100 (MQ=255)                                        gggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCgg                                                         <  1:1896741/99‑1 (MQ=255)                                         ggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGa                                                       >  2:268858/1‑100 (MQ=255)                                         ggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGATAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGa                                                       >  2:672983/1‑100 (MQ=255)                                               aaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTa                                                 >  2:1918420/1‑100 (MQ=255)                                               aaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTa                                                 >  2:129848/1‑100 (MQ=255)                                               aaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTa                                                 <  1:1353386/100‑1 (MQ=255)                                                aTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTaa                                                >  1:950734/1‑100 (MQ=255)                                                 ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTaa                                                >  1:1952383/1‑99 (MQ=255)                                                 ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                               <  1:1936397/100‑1 (MQ=255)                                                 ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                               <  2:2041199/100‑1 (MQ=255)                                                 ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                               >  2:2105609/1‑100 (MQ=255)                                                 ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                               <  2:2154255/100‑1 (MQ=255)                                                 ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                               <  1:605736/100‑1 (MQ=255)                                                    gACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACg                                            <  1:1407047/100‑1 (MQ=255)                                                    gACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACg                                            >  1:1379697/1‑100 (MQ=255)                                                     aCAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACg                                            <  2:2264167/99‑1 (MQ=255)                                                       aaTGCCGCGATAGTTATAACCGCTAACATTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGATGCGGGGTAAGATCCGGGACTCGTAAGACGaaa                                         >  2:563580/1‑100 (MQ=255)                                                        aTGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaaaa                                        >  2:1335710/1‑100 (MQ=255)                                                        aTGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaaaa                                        <  1:2034062/100‑1 (MQ=255)                                                        aTGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGGGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaaaa                                        <  1:507762/100‑1 (MQ=255)                                                           ccGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGc                                      <  1:341749/99‑1 (MQ=255)                                                           ccGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCt                                     >  1:973679/1‑100 (MQ=255)                                                           ccGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCt                                     >  1:869681/1‑100 (MQ=255)                                                            cgcgATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCt                                     <  1:1882959/99‑1 (MQ=255)                                                              cgATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTa                                   <  1:1071482/99‑1 (MQ=255)                                                               gATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATATGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTAt                                  >  2:1293737/1‑99 (MQ=255)                                                                  aGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGTGGTAATATCCGGGACTCGTAAGACGAAAAGCTTATCACt                              <  2:1893897/100‑1 (MQ=255)                                                                  aGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTTGTAAGACGAAAAGCTTATCACt                              >  2:1023020/1‑100 (MQ=255)                                                                  aGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACt                              >  1:896179/1‑100 (MQ=255)                                                                     tataACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCTGGACTCGTAAGACGAAAAGCTTATCACTATg                           >  2:885247/1‑100 (MQ=255)                                                                     tataACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATg                           <  1:1312402/100‑1 (MQ=255)                                                                       taACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCt                         >  2:561113/1‑100 (MQ=255)                                                                         aCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTc                       >  1:703288/1‑100 (MQ=255)                                                                          ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa                      <  2:675485/100‑1 (MQ=255)                                                                          ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa                      >  2:231212/1‑100 (MQ=255)                                                                          ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa                      <  1:1580966/100‑1 (MQ=255)                                                                           cGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCAt                     >  2:1867224/1‑100 (MQ=255)                                                                           cGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCAt                     <  1:579317/100‑1 (MQ=255)                                                                              tATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGc                  <  1:1969504/100‑1 (MQ=255)                                                                                  tttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttt              <  1:2302599/100‑1 (MQ=255)                                                                                  tttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttt              <  2:756031/100‑1 (MQ=255)                                                                                  nttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGTtttt              >  2:2381696/2‑100 (MQ=255)                                                                                    tGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttttt            <  2:664647/100‑1 (MQ=255)                                                                                    tGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttttt            <  1:2381294/100‑1 (MQ=255)                                                                                      aTTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTaa          >  1:1186930/1‑100 (MQ=255)                                                                                       ttttATTGTTAAACGCTGTGATATCTGTGGTTTGTTGAGGGGTAAGATCCGGGACTCTTAAGACGAAAAGCTGATTTTTATGTTTGATCGCTTTTTTAAg         <  2:723904/100‑1 (MQ=255)                                                                                        tttATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTCATCACTATGCTTCATCGCTTTTTTAAgg        >  2:2098019/1‑100 (MQ=255)                                                                                        tttAGTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGCCACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAgg        <  2:129061/100‑1 (MQ=255)                                                                                        tctATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAgg        <  2:1447592/98‑1 (MQ=255)                                                                                          tATTGGTAAACGCGGTGATATCTGTTGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGa      <  1:561116/100‑1 (MQ=255)                                                                                             tGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAAc   >  1:1469471/1‑100 (MQ=255)                                                                                              ggCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACg  >  1:588114/1‑100 (MQ=255)                                                                                                |                                                                                                  GAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAACG  >  NC_000913/216473‑216665

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: ‑