Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,018,254	C→A	15.8%	D16Y (GAT→TAT)	ycbZ ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,018,254	0	C	A	15.8%	266.0 / 31.7	95	D16Y (GAT→TAT)	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base C (37/43);  new base A (9/6);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

TGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACT  >  NC_000913/1018157‑1018353                                                                                                  |                                                                                                    tGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCgg                                                                                                   <  1:1838247/100‑1 (MQ=255) tGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCgg                                                                                                   >  1:1300737/1‑100 (MQ=255)  gCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGt                                                                                                  >  1:10653/1‑100 (MQ=255)  gCTCCAGCGCAAAGTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGt                                                                                                  >  1:1749175/1‑100 (MQ=255)   cTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTa                                                                                                 >  1:702936/1‑100 (MQ=255)    tCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTAt                                                                                                >  1:2562325/1‑100 (MQ=255)     ccAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATc                                                                                               <  1:187579/100‑1 (MQ=255)     ccAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATc                                                                                               >  2:1485476/1‑100 (MQ=255)       agtgcAAATTGCAGCCGCGGTTGAGTATCCCTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAg                                                                                             <  2:1302344/97‑1 (MQ=255)       aGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAg                                                                                             >  2:1334888/1‑100 (MQ=255)         cgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGa                                                                                           >  1:890631/1‑100 (MQ=255)         cgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGa                                                                                           <  1:987320/100‑1 (MQ=255)          gcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaa                                                                                          <  2:1744074/100‑1 (MQ=38)           cAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaa                                                                                         <  2:2137296/100‑2 (MQ=38)           cAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaa                                                                                         <  2:1325591/100‑2 (MQ=38)            aaaTTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAAcc                                                                                        <  1:1626982/100‑1 (MQ=255)            aaaTTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAAcc                                                                                        <  2:1949651/100‑1 (MQ=255)              aTTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAg                                                                                      >  1:2272947/1‑100 (MQ=255)               ttGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAgg                                                                                     <  2:1082216/100‑1 (MQ=255)               ttGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAgg                                                                                     <  2:1081897/100‑1 (MQ=255)               gtGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGTAAATATTTCCTGATAGCTATCGGTATCAGGAACCAgg                                                                                     <  2:1393395/99‑1 (MQ=255)                tggaGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGt                                                                                    <  1:1674391/97‑1 (MQ=255)                 gCAGCCGCGGTTGAATATCACGGAAAAAAAGATCGTTTGCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaatatttg                                                                                   <  2:1916875/100‑8 (MQ=21)                  cAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaatatttgc                                                                                  <  2:2537783/100‑9 (MQ=25)                   aGCCGCGCTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCAc                                                                                 >  1:1371875/1‑100 (MQ=255)                    gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCAc                                                                                 >  2:86081/1‑99 (MQ=255)                    gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACg                                                                                <  2:2487983/100‑1 (MQ=255)                    gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACg                                                                                >  1:877794/1‑100 (MQ=255)                    gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACg                                                                                <  1:825821/100‑1 (MQ=255)                        cgGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCAt                                                                            <  1:1435563/100‑1 (MQ=255)                         ggTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATg                                                                           <  2:820781/100‑1 (MQ=255)                             gAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAg                                                                       <  2:130044/100‑1 (MQ=255)                             gAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAg                                                                       <  2:167638/100‑1 (MQ=255)                              aatATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGt                                                                      >  1:1754245/3‑100 (MQ=255)                              aGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGt                                                                      <  2:778686/100‑1 (MQ=255)                                tATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGttt                                                                    <  1:726085/100‑1 (MQ=255)                                tATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGttt                                                                    >  1:1219259/1‑100 (MQ=255)                                  tCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCg                                                                  <  1:1358373/100‑1 (MQ=255)                                   cACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGGCACGCCATGCAAGTTTCGt                                                                 <  2:404917/100‑1 (MQ=255)                                    aCTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTa                                                                >  2:165517/1‑100 (MQ=255)                                    aCTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTa                                                                >  2:711921/1‑100 (MQ=255)                                       gAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATg                                                             >  1:1423414/1‑100 (MQ=255)                                        aaTAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATgg                                                            <  2:1748119/100‑1 (MQ=255)                                        aaTAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTGCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATgg                                                            <  2:2600459/100‑1 (MQ=255)                                        aaTAAAGGATCGTTTTCGTCAAACAAATGTGGCTGAGCAAATATCTCCTGATAGCTAGCGGGAACAGGAACCAGGTCACGCCATGCAAGTGTCGTAATgg                                                            <  2:1178204/100‑1 (MQ=255)                                          tAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTc                                                          >  2:2273491/1‑100 (MQ=255)                                            aaGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCaa                                                        <  2:2364753/100‑1 (MQ=255)                                                 tCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTg                                                   >  1:2438231/1‑100 (MQ=255)                                                  cGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGa                                                  <  2:1016590/100‑1 (MQ=255)                                                   gTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGAt                                                 >  1:888490/1‑100 (MQ=255)                                                    ttttCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATg                                                <  2:2182742/100‑1 (MQ=255)                                                       tCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGttt                                             >  1:2491541/1‑100 (MQ=255)                                                        cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                            >  1:1881230/1‑100 (MQ=255)                                                        cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                            >  2:2046726/1‑100 (MQ=255)                                                        cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                            <  1:2273504/100‑1 (MQ=255)                                                        cGTAAATCAAATGTGGCTGAGCAAAGATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                            <  2:633925/100‑1 (MQ=255)                                                         gTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGttttt                                           <  1:822345/100‑1 (MQ=255)                                                          tcaatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttttt                                          <  2:890626/100‑1 (MQ=255)                                                           caatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTa                                         >  1:1208661/1‑100 (MQ=255)                                                           caatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCAAGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTa                                         >  2:2245736/1‑100 (MQ=255)                                                           caatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTAACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTa                                         >  2:697835/1‑100 (MQ=255)                                                            aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGTCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                        <  2:917680/100‑1 (MQ=255)                                                            aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                        <  2:894247/100‑1 (MQ=255)                                                            aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                        <  1:1257913/100‑1 (MQ=255)                                                            aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                        <  2:888485/100‑1 (MQ=255)                                                            aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                        <  1:1670483/100‑1 (MQ=255)                                                                        atGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCaaa                            >  2:764055/2‑100 (MQ=255)                                                                          gAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAg                           <  1:2451360/99‑1 (MQ=255)                                                                           aGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGcc                         >  1:2401225/1‑100 (MQ=255)                                                                             gatctggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgc                       <  2:2528495/93‑1 (MQ=37)                                                                             cAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgc                       >  2:501172/1‑100 (MQ=255)                                                                              acctggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcg                      >  2:896480/7‑100 (MQ=37)                                                                              aaaTATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcg                      <  2:2052252/100‑1 (MQ=255)                                                                                ctggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcg                      >  1:1819239/5‑98 (MQ=37)                                                                                 tggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGAtt                   >  1:1909717/4‑100 (MQ=38)                                                                                 tctttCCTGATACCTATCGGTATCAGGAACCCGGTCACGACATGCAAGTTTCGTAATGGTCAAGGTGGATGTTTTCTAGTCTGTTTTCAAagcagcaa                     >  2:1749228/3‑96 (MQ=255)                                                                                 tatTGCCTGATAGCTATCGGTATCCGGAAACGGGTCACGCCATGGAAGGTTCGTAATGGTCCAAGTACATGTTTTTTGGACTGGTGTCTAAGCCGCGac                    <  2:699831/99‑2 (MQ=255)                                                                                  ggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcagatc                  >  1:2528508/3‑95 (MQ=37)                                                                                   tttCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat                 >  2:1598653/1‑100 (MQ=255)                                                                                   gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat                 >  1:1385636/2‑100 (MQ=38)                                                                                   gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat                 >  1:19450/2‑100 (MQ=38)                                                                                   gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat                 >  2:1816321/2‑100 (MQ=38)                                                                                    ttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtata                >  1:576623/1‑100 (MQ=38)                                                                                     tCCTGATACCGATAGCTATCAGGAACAAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtata                >  2:485489/1‑99 (MQ=38)                                                                                      ccTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAcc              >  2:682821/1‑100 (MQ=255)                                                                                      ccTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAcc              >  1:1982218/1‑100 (MQ=255)                                                                                       cTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCg             <  2:1939596/100‑1 (MQ=255)                                                                                          aTAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTaa          >  1:1633114/1‑100 (MQ=255)                                                                                           tAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAAc         >  1:176387/1‑100 (MQ=255)                                                                                             gCTATCGGTATCAGGAACC‑GGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCgg      >  2:297990/1‑100 (MQ=255)                                                                                               tATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGc     <  1:2288704/100‑1 (MQ=255)                                                                                               tATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGc     <  1:2457789/100‑1 (MQ=255)                                                                                               tATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGc     <  2:1448733/100‑1 (MQ=255)                                                                                               tATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTAGTAATGGTCAAAGTTGATGTTTTTCAGTCTGTTGTCAAAGCCGCGATTATCCCGTAACCg       >  2:1533594/1‑98 (MQ=255)                                                                                                 tCGGTATCAGGAACCAGGTCACGTCATGCAAGTTTCGTAATGGGCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCAc   <  2:2499909/100‑1 (MQ=255)                                                                                                  cGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACt  >  1:210615/1‑100 (MQ=255)                                                                                                  |                                                                                                    TGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACT  >  NC_000913/1018157‑1018353

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