Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,018,256	G→C	16.0%	T15S (ACC→AGC)	ycbZ ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,018,256	0	G	C	16.0%	260.7 / 42.4	94	T15S (ACC→AGC)	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base G (37/42);  new base C (9/6);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

TGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTAC  >  NC_000913/1018157‑1018355                                                                                                    |                                                                                                    tGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCgg                                                                                                     >  1:1300737/1‑100 (MQ=255) tGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCgg                                                                                                     <  1:1838247/100‑1 (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:987320/100‑1 (MQ=255)         cgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGa                                                                                             >  1:890631/1‑100 (MQ=255)          gcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaa                                                                                            <  2:1744074/100‑1 (MQ=38)           cAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaa                                                                                           <  2:1325591/100‑2 (MQ=38)           cAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaa                                                                                           <  2:2137296/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                                                                                  >  1:877794/1‑100 (MQ=255)                    gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACg                                                                                  <  1:825821/100‑1 (MQ=255)                    gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACg                                                                                  <  2:2487983/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:1219259/1‑100 (MQ=255)                                tATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGttt                                                                      <  1:726085/100‑1 (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                                              >  2:2046726/1‑100 (MQ=255)                                                        cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                              >  1:1881230/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:19450/2‑100 (MQ=38)                                                                                   gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat                   >  1:1385636/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                >  1:1982218/1‑100 (MQ=255)                                                                                      ccTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAcc                >  2:682821/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)                                                                                                    gTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTAc  >  1:2539306/1‑100 (MQ=255)                                                                                                    |                                                                                                    TGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTAC  >  NC_000913/1018157‑1018355

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