Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,018,251	T→G	16.2%	S17R (AGC→CGC) ‡	ycbZ ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,018,251	0	T	G	16.2%	206.1 / 38.4	91	S17R (AGC→CGC)	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (37/39);  new base G (9/6);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

AACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGC  >  NC_000913/1018154‑1018350                                                                                                  |                                                                                                    aaCTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTAt                                                                                                   >  1:1448741/1‑100 (MQ=255)    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                                                                                            >  2:1485476/1‑100 (MQ=255)        ccAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATc                                                                                            <  1:187579/100‑1 (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: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:1081897/100‑1 (MQ=255)                  ttGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAgg                                                                                  <  2:1082216/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: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:711921/1‑100 (MQ=255)                                       aCTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTa                                                             >  2:165517/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:2273504/100‑1 (MQ=255)                                                           cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                         >  1:1881230/1‑100 (MQ=255)                                                           cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                         >  2:2046726/1‑100 (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)                                                                                                  |                                                                                                    AACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGC  >  NC_000913/1018154‑1018350

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