Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,018,251	T→G	15.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	15.2%	225.0 / 36.0	92	S17R (AGC→CGC)	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (41/37);  new base G (7/7);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

CAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGC  >  NC_000913/1018153‑1018350                                                                                                   |                                                                                                    caACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTa                                                                                                    >  2:2008812/1‑100 (MQ=255)   aCTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATc                                                                                                  <  1:2502960/100‑1 (MQ=255)   aCTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATc                                                                                                  <  1:2249191/100‑1 (MQ=255)   aCTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATc                                                                                                  <  2:1801218/100‑1 (MQ=255)     tGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCgg                                                                                                <  1:2517977/100‑1 (MQ=255)     tGCTCCAGCGCAAATGGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCgg                                                                                                <  1:496172/100‑1 (MQ=255)      gCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGt                                                                                               >  2:1302561/1‑100 (MQ=255)      gCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCATCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGt                                                                                               >  1:699233/1‑100 (MQ=255)       cTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGCTAGCTATCGGTa                                                                                              >  1:2153624/1‑100 (MQ=255)        tCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTAt                                                                                             >  1:2585205/1‑100 (MQ=255)        tCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTAt                                                                                             >  1:2007481/1‑100 (MQ=255)         ccAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATc                                                                                            <  1:1427923/100‑1 (MQ=255)          cAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCa                                                                                           <  2:1306325/100‑1 (MQ=255)           aGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAg                                                                                          <  2:2585224/100‑1 (MQ=255)           aGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAg                                                                                          >  2:196463/1‑100 (MQ=255)            gcgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAgg                                                                                         <  1:725338/100‑1 (MQ=255)             cgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGa                                                                                        >  1:324392/1‑100 (MQ=255)             cgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGa                                                                                        <  1:1145464/100‑1 (MQ=255)             cgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGa                                                                                        <  1:488168/100‑1 (MQ=255)              gcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaa                                                                                       <  1:1303945/100‑1 (MQ=38)              gcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaa                                                                                       <  2:1003529/100‑1 (MQ=38)                   ttGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATCTTTCCTGATAGCTATCGGTATGAGGAACCa                                                                                    >  1:1234332/1‑98 (MQ=255)                    tGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGt                                                                                 >  2:2586244/1‑100 (MQ=255)                    tGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGt                                                                                 >  1:1296128/1‑100 (MQ=255)                     gCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGTGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTc                                                                                <  1:1412992/100‑1 (MQ=255)                     gCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTc                                                                                >  1:1262476/1‑100 (MQ=255)                      cAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCa                                                                               >  2:984419/1‑100 (MQ=255)                      cAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaatatttgc                                                                               <  2:2567949/100‑9 (MQ=25)                       aGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaatatt                                                                                  <  2:954033/96‑6 (MQ=37)                        gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACg                                                                             <  1:984416/100‑1 (MQ=255)                         ccGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGc                                                                            <  2:2291044/100‑1 (MQ=255)                         ccGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGc                                                                            >  2:862031/1‑100 (MQ=255)                          cgcgGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGcc                                                                           >  2:1370070/1‑100 (MQ=255)                             ggatGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATCTGCCTGACCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATg                                                                        <  2:117366/97‑1 (MQ=255)                             ggTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATg                                                                        >  2:529546/1‑100 (MQ=255)                             ggTGGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATg                                                                        <  1:2359196/100‑1 (MQ=255)                                   gTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGtt                                                                  <  1:352859/100‑1 (MQ=255)                                     aTCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTc                                                                >  1:2138317/1‑100 (MQ=255)                                     aTCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTc                                                                <  2:1163062/100‑1 (MQ=255)                                      tCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCg                                                               >  2:1894842/1‑100 (MQ=255)                                      tCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCg                                                               <  2:1289186/100‑1 (MQ=255)                                      tCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCg                                                               <  1:1143892/100‑1 (MQ=255)                                        aCTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTa                                                             >  2:307122/1‑100 (MQ=255)                                         cTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTaa                                                            >  2:2347355/1‑100 (MQ=255)                                           gAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATg                                                          >  1:1844837/1‑100 (MQ=255)                                             aTAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGt                                                        >  2:2511688/1‑100 (MQ=255)                                                  ggATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAg                                                   >  2:1483455/1‑100 (MQ=255)                                                     tCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTg                                                >  2:958638/1‑100 (MQ=255)                                                     tCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTg                                                >  2:1891298/1‑100 (MQ=255)                                                      cGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGa                                               >  1:912187/1‑100 (MQ=255)                                                       gTTTTTGTGAATCAAATGTGGCTGGGCAAAGATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGAt                                              <  1:1766898/100‑1 (MQ=255)                                                       gTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGAt                                              <  1:1767290/100‑1 (MQ=255)                                                        ttttCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACAAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATg                                             >  2:1324712/1‑100 (MQ=255)                                                        ttttCGTCAATCAAAGGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGAtn                                             <  1:2089328/100‑2 (MQ=255)                                                            cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttt                                         >  2:2004039/1‑100 (MQ=255)                                                              tcaatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttttt                                       >  2:2312976/1‑100 (MQ=255)                                                              tcaatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttttt                                       <  2:721434/100‑1 (MQ=255)                                                               caatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTa                                      >  2:1229254/1‑100 (MQ=255)                                                               caatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTa                                      >  2:2279612/1‑100 (MQ=255)                                                                aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                     >  2:2061878/1‑100 (MQ=255)                                                                  tcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTc                                   >  2:1760707/1‑100 (MQ=255)                                                                    aaaTGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCtg                                 >  2:949664/1‑100 (MQ=255)                                                                     aaTGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCtgt                                <  1:1575567/100‑1 (MQ=255)                                                                        gtgGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCtgttgt                             >  2:162510/1‑100 (MQ=255)                                                                         tgGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAATTTGATGTTTTTTAGTCTGTTGTc                            <  2:442275/100‑1 (MQ=255)                                                                          ggCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCa                           <  2:2579171/100‑1 (MQ=255)                                                                          ggCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCa                           <  2:105030/100‑1 (MQ=255)                                                                              gAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGc                       <  2:1534015/100‑1 (MQ=255)                                                                              gAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGc                       <  2:2145233/100‑1 (MQ=255)                                                                              gAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGc                       <  1:2483720/100‑1 (MQ=255)                                                                              gAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAAGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGc                       >  2:2480009/1‑100 (MQ=255)                                                                              gAGCAAATATTTCCTGATAGCTATCGGTAGCAGGAACCAGGTCACGCTATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGc                       <  1:1930912/100‑1 (MQ=255)                                                                                 cAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgc                    >  2:301393/1‑100 (MQ=255)                                                                                  acctggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcg                   >  1:231938/7‑100 (MQ=37)                                                                                   aaTATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcg                   <  1:724366/99‑1 (MQ=255)                                                                                   aaTATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGa                  <  1:856925/100‑1 (MQ=255)                                                                                    ctggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGAt                 >  2:479611/5‑100 (MQ=37)                                                                                      ggttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATta               >  1:98437/3‑100 (MQ=38)                                                                                       gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTTGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat              >  1:872584/2‑100 (MQ=38)                                                                                       gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat              >  1:1827900/2‑100 (MQ=38)                                                                                       gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat              >  1:932638/2‑100 (MQ=38)                                                                                        ttCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtata             <  2:329941/100‑1 (MQ=255)                                                                                         tCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAc            <  1:2342557/100‑1 (MQ=255)                                                                                         tCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCCAGTTTCGTACTGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAc            >  1:583111/1‑100 (MQ=255)                                                                                          ccTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAcc           >  1:1734368/1‑100 (MQ=255)                                                                                           cTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCg          >  1:1978414/1‑100 (MQ=255)                                                                                           cTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCg          <  2:699236/100‑1 (MQ=255)                                                                                           cTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCg          <  1:1041925/100‑1 (MQ=255)                                                                                              aTAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTaa       <  2:10758/100‑1 (MQ=255)                                                                                                 gCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCg    <  1:1788114/100‑1 (MQ=255)                                                                                                 gCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCg    >  2:1601318/1‑100 (MQ=255)                                                                                                 gCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCg    >  1:1473251/1‑100 (MQ=255)                                                                                                   tATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGc  >  2:1230154/1‑100 (MQ=255)                                                                                                   tATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGc  <  2:1003487/100‑1 (MQ=255)                                                                                                   |                                                                                                    CAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGC  >  NC_000913/1018153‑1018350

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