Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,401,459	0	T	C	27.8%	126.3 / 22.1	72	D625G (GAT→GGT)	csrD	regulator of CsrB and CsrC decay
Reads supporting (aligned to +/- strand):  ref base T (38/14);  new base C (0/20);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 4.46e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTG  >  NC_000913/3401321‑3401586                                                                                                                                           |                                                                                                                                gCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATc                                                                                                                                >  2:127347/1‑140 (MQ=255)      cAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCtt                                                                                                                           >  1:552645/1‑140 (MQ=255)        gACAAGCACACTCCCGTGAGCCGCAGCGCGCATTTTTCTACGTGAAAACGGATTAAACGGAGGTTTAAACCGAGTATTTTTGTGAATATTTTTTCACGTTAGTATAAAGTGCCTGTGAGGACGCAAAAAAACCCCCTTGc                                                                                                                         <  2:251788/140‑1 (MQ=255)        gACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAGGTGGCTGTGAGGACGCGAAAAAACCCCCTTGc                                                                                                                         <  1:577017/140‑1 (MQ=255)          cAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGccc                                                                                                                       <  1:1124640/140‑1 (MQ=255)              cacaCTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcc                                                                                                                               >  2:776570/1‑128 (MQ=255)              cacaCTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcc                                                                                                                               <  1:687644/128‑1 (MQ=255)              cacaCTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcc                                                                                                                               >  2:687644/1‑128 (MQ=255)              cacaCTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcc                                                                                                                               <  1:776570/128‑1 (MQ=255)              cacaCTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCt                                                                                                                   >  1:688556/1‑140 (MQ=255)               ccccTCCCAGGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAACCGGCAGGTTAAACGGAGTATCTTTGTAAATATTTTTTCCCGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTg                                                                                                                  <  2:1196340/137‑1 (MQ=255)                   tCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGc                                                                                                                         <  1:211511/129‑1 (MQ=255)                   tCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGc                                                                                                                         >  2:211511/1‑129 (MQ=255)                    cccATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCCCGTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAATCCCCTTGCCCGCCTGTAaca                                                                                                             <  2:787094/140‑1 (MQ=255)                       aTGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcg                                                                                                          >  1:424418/1‑140 (MQ=255)                       aTGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAATCCCCTTGccc                                                                                                                       >  1:423368/1‑127 (MQ=255)                       aTGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAATCCCCTTGccc                                                                                                                       <  2:423368/127‑1 (MQ=255)                        tGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcgc                                                                                                         <  2:688556/140‑1 (MQ=255)                         gagaCGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACcgcg                                                                                                        <  2:419361/140‑1 (MQ=255)                         gagaCGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAATCCCCTTGCCCGCCTGTAACACcgcg                                                                                                        <  2:1150796/140‑1 (MQ=255)                             cGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGa                                                                                                    >  2:555074/1‑140 (MQ=255)                               cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAAt                                                                                                  >  2:556405/1‑140 (MQ=255)                                aGCGCGCATTATTTTACGTAAAAACGGATAAACCGGCAGTTAAACCCGATTATTTTTGGAAATATTTTTCCACGTTAGTATCAAGGGGCTGGGAGGACGCAAAAAATTCCCCTTCCCCGCCTGTAACACCGCGCTGAATc                                                                                                 <  2:861997/140‑1 (MQ=255)                                  agacagttattCTACGTAAAAACGGATTAACCGGCGGATTAAACCGATTATCTTTGTGAATATTTTTTCCCTTTAGTATCAGGTGGCTGTGAGGACGCAAAAAAACCCCCTGGCCCGCCTGTAACACCGCGCTGAATCAg                                                                                               <  2:254653/134‑1 (MQ=255)                                    cgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAggg                                                                                             >  1:518659/1‑140 (MQ=255)                                    cgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCGGTGAGGACGCGAAAAGATCCCCTTTCCCGCCTGTAACACCGCGCTGAATCAggg                                                                                             >  1:99906/1‑140 (MQ=255)                                       attattCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCt                                                                                          <  1:1012150/140‑1 (MQ=255)                                             cTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACt                                                                                    <  2:424418/140‑1 (MQ=255)                                                 gTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCt                                                                                >  2:577953/1‑140 (MQ=255)                                                 gTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCt                                                                                >  1:265504/1‑140 (MQ=255)                                                 gTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCt                                                                                >  1:165143/1‑140 (MQ=255)                                                   gAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTctgtct                                                                              >  1:254653/1‑135 (MQ=255)                                                   gAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTc                                                                              >  2:171123/1‑140 (MQ=255)                                                     aaaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTcgc                                                                            >  2:267406/1‑140 (MQ=255)                                                      aaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTcgcg                                                                           <  1:127347/140‑1 (MQ=255)                                                           aTTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAacgc                                                                      >  2:692726/1‑140 (MQ=255)                                                            ttAACCGGCAGGTTAACCCGAGTATCTTTGTGAATTTTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAacgca                                                                     <  2:165143/140‑1 (MQ=255)                                                                         aaacccGATTATTTTTGTGAATTTTTTTTCACTTTAGTATCAGGTGGCTGTGAGGACGCAAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTa                                                        <  2:979400/136‑1 (MQ=255)                                                                          aaaCCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTaa                                                       >  2:970259/1‑140 (MQ=255)                                                                             cggAGTATTTTTGTAAATATTTTTTCACGTTAGTATCAAGGGGCTGTGGGGACGCGAAAAACCCCCCTCGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAAcc                                                    <  2:552645/138‑1 (MQ=255)                                                                               gAGTATTTTTGTAAATTTTTTTTCACTTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTg                                                  <  2:1211393/140‑1 (MQ=255)                                                                                aGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTgg                                                 >  1:621628/1‑140 (MQ=255)                                                                                 gTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTggg                                                >  2:1130125/1‑140 (MQ=255)                                                                                 gTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTggg                                                >  1:545969/1‑140 (MQ=255)                                                                                  tATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGt                                               >  1:1211393/1‑140 (MQ=255)                                                                                  tATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGt                                               >  1:833451/1‑140 (MQ=255)                                                                                   ttttttGTAAATTTTTTTTCACTTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTg                                              <  2:265504/137‑1 (MQ=255)                                                                                    tCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGc                                             <  1:73158/140‑1 (MQ=255)                                                                                     ttttGTGAATTTTTTTCCACTTTAGTATAAAGTGGCTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCt                                            <  2:621628/139‑1 (MQ=255)                                                                                        tgtgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGt                                         <  1:692726/140‑1 (MQ=255)                                                                                          tgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTcc                                       >  1:348646/1‑140 (MQ=255)                                                                                          tgAATATTTTTTCACGTTAGTATCAAGTGGCAGGGAGGAGGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTcc                                       <  1:555074/140‑1 (MQ=255)                                                                                            aaTATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGATTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCg                                     >  1:1000330/1‑140 (MQ=255)                                                                                               atTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAg                                  <  1:970259/140‑1 (MQ=255)                                                                                               atTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAg                                  <  1:361247/140‑1 (MQ=255)                                                                                                 tttttCCCGTTAGTATCAGGTGGCTGTGAGGACGCAAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCa                                <  2:382487/140‑1 (MQ=255)                                                                                                 tttttCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCt                                            >  2:135436/1‑128 (MQ=255)                                                                                                 tttttCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCt                                            <  1:135436/128‑1 (MQ=255)                                                                                                    ttCCCGTTAGTATCAAGTGGCTGTGAGGCCGCAAAAAAACCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGATTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGc                             <  2:1000330/140‑1 (MQ=255)                                                                                                          ttAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTccac                       >  1:473969/1‑140 (MQ=255)                                                                                                            aGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTccacca                     <  2:473969/140‑1 (MQ=255)                                                                                                                 aaggTGGTTGGGAGGCCGCAAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCtt                <  2:85456/137‑1 (MQ=255)                                                                                                                   aGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCa                                >  2:564027/1‑122 (MQ=255)                                                                                                                   aGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCa                                <  1:564027/122‑1 (MQ=255)                                                                                                                     tGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGaa            >  2:59151/1‑140 (MQ=255)                                                                                                                     tGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGaa            >  1:1105354/1‑140 (MQ=255)                                                                                                                     tGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGaa            >  2:1159810/1‑140 (MQ=255)                                                                                                                      ggCTGTGAGGCCGCAAAAAAACCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAc           <  2:821527/140‑1 (MQ=255)                                                                                                                       gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc          >  2:665811/1‑140 (MQ=255)                                                                                                                        cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCttctgtctc         >  1:2393/1‑133 (MQ=255)                                                                                                                        cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca         >  1:720582/1‑140 (MQ=255)                                                                                                                        cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca         >  1:821527/1‑140 (MQ=255)                                                                                                                        cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca         >  1:85456/1‑140 (MQ=255)                                                                                                                              ggACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCt   <  1:1089924/140‑1 (MQ=255)                                                                                                                               gACGCAAAAAACCCCCATAACCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTg  <  1:267406/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                GCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTG  >  NC_000913/3401321‑3401586

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