Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,492,834	0	T	G	41.9%	118.6 / 21.3	76	G89G (GGT→GGG)	tsgA	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (6/37);  new base G (31/0);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 6.98e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAC  >  NC_000913/3492686‑3492977                                                                                                                                                     |                                                                                                                                                   tGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGgtt                                                                                                                                                  >  2:383344/1‑151 (MQ=255)  gCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGttt                                                                                                                                                  >  2:235030/1‑150 (MQ=255)      gTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATg                                                                                                                                              <  2:930781/150‑1 (MQ=255)       tttCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGGTGCCGGgttgttgtt                                                                                                                                            >  1:83174/1‑151 (MQ=255)        ttCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGGCGGgtt                                                                                                                                                  >  1:823288/1‑144 (MQ=255)            aGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGcc                                                                                                                                       <  2:727734/151‑1 (MQ=255)             gTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCca                                                                                                                                      <  2:63605/151‑1 (MQ=255)              tATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGGTGGGTTGATGttttggca                                                                                                                                      >  2:1387921/1‑144 (MQ=255)                  aGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGcc                                                                                                                                 <  1:677339/151‑1 (MQ=255)                   gTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCt                                                                                                                                <  1:782893/151‑1 (MQ=255)                       tACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTCCGGGTTTATGTTTAGCCCCAGCCTgg                                                                                                                              >  2:1316386/1‑149 (MQ=255)                        atcttcacattaGTAAACGACGGAATTTTAATCTCTATCTTGCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAATGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGGCGGGTTGATGTTTAgagaaagactgggg                                                                                                                            >  1:291107/3‑137 (MQ=255)                         ccttcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCt                                                                                                                          <  2:842575/151‑1 (MQ=255)                            tcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTTAGCCACAACCTGGGGCTGtt                                                                                                                       >  2:1246104/1‑151 (MQ=255)                             caccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACAGCCTGGCGCTGtt                                                                                                                       >  2:220680/1‑150 (MQ=255)                               ccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTtctc                                                                                                                    <  1:956082/151‑1 (MQ=255)                                cttcCTCACCGCCGGCATTTTAGTCACTAACTACCTCAACGCCCGGCTGATGGAAATCGCCCCGGTGAGAACCCAGCTACGTTTTGGCTTTCCCCTGATGGTGCTGGCGGTTGCCAGCTTGATGTTCAGCCACGGCCTGGCTCTTTtct                                                                                                                     <  1:416916/149‑1 (MQ=255)                                  tcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTcgg                                                                                                                  <  2:610192/150‑1 (MQ=255)                                  tcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGAGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTcg                                                                                                                   <  2:70225/149‑1 (MQ=255)                                   ccTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTcggc                                                                                                                 <  2:233572/150‑1 (MQ=255)                                       aaCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTGTGTTTAGCCACCGCCTGGCGCTGGTCTcggcgg                                                                                                               >  1:1016659/1‑148 (MQ=255)                                       aaCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGATGCCGGGTTTTTGTTCAGCCACACCCCCGCGCGGTTCCCGgcggcg                                                                                                             >  1:889899/1‑150 (MQ=255)                                        aCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGAt                                                                                                           <  2:189145/151‑1 (MQ=255)                                          gCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGATGCGGGGTTGATGTTCAGCCACAGCCTGTCGCTGTTCTcggcgg                                                                                                               >  1:822083/1‑145 (MQ=255)                                               cATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCTCTGTTCTCGGCTGCGATGTTCa                                                                                                      <  1:306525/149‑1 (MQ=255)                                                aTTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATtc                                                                                                   <  2:983449/151‑1 (MQ=255)                                                aTTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGGTGTTCAGCCCCAGCCCGGCGCTGTTCTCGGCGGCGGTGTTTAtt                                                                                                    >  2:372219/1‑150 (MQ=255)                                                aTTTTAATCTCTATCTTACTCAAAGAATGGCTGATGGAAATCGTCACGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGAGGTTGCCGGGTTGATGTTCAGGCAAAGCCTGGGGCTGTTTTTGGTGGGGGTGtttttt                                                                                                    >  1:705020/1‑146 (MQ=255)                                                    tAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAGCGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCg                                                                                               <  1:118617/150‑1 (MQ=255)                                                     aaTCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCgg                                                                                              <  1:250400/150‑1 (MQ=255)                                                       tctcTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATTTTATGACCAAGCCTGGGGGTGCTCTTGGCGGGGGTGTTCATTT‑TTggg                                                                                             >  1:900214/1‑149 (MQ=255)                                                           tATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGGTTTCGGCGGGGATGTTCATTT‑TCGgggggg                                                                                         >  2:815729/1‑149 (MQ=255)                                                           tATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGGTGCCGGGTTGATGTTCAGCCACAGCCCGGCGCTGGTTTCCGCGGGGATGGTTATTA‑TCGGggagg                                                                                         >  1:562631/1‑149 (MQ=255)                                                            aTCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGggtgg                                                                                         <  1:1126231/148‑1 (MQ=255)                                                            aTCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACCGCCTGGCGCTGTTCTCGGCGGGGGGGTTCTTTC‑TCGgggggg                                                                                         >  1:739346/1‑148 (MQ=255)                                                            aTCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTCAGCCACAACCTGGCGCTGTTCTCGGCGGGGGTGTTCATTC‑TCGgggggg                                                                                         >  2:459669/1‑148 (MQ=255)                                                             tCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGgggggg                                                                                         >  2:179023/1‑147 (MQ=255)                                                             tCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGACGGGTTTATGTTTTGCCACCGCCCGGCGCTGTTCTCCGCGGGGGGGGTTCTTC‑TCGgggggg                                                                                         >  1:805261/1‑147 (MQ=255)                                                              cTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTCCGGGTTTATGTTCAGCCCCCGCCTGGCGCTGTTTTCGGCGGCGGGGGTTATTC‑TTGgggggg                                                                                         >  1:1321245/1‑146 (MQ=255)                                                              cTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTTTGTTCAGCCACCACCTGGCGCTGTTCTCGGCGGCGGAGGTTATTT‑TCGGGGGGGGCa                                                                                      >  1:108276/1‑149 (MQ=255)                                                               ttCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGTTGTTCAGCCACCGCCTGGCGCTGTTCTCGGCGGCGAGGTTTTTTC‑TCGGGGGGGTGAg                                                                                     >  1:108213/1‑149 (MQ=255)                                                               ttCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGGTCAGCCACCGCCTGGCGATGTTTTCGGCGGGGGTGTTTATTC‑CCGGGGGGGGCAGCg                                                                                   >  1:1252358/1‑151 (MQ=255)                                                                tCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCg                                                                                   <  2:760698/150‑1 (MQ=255)                                                                    cAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCAt                                                                               <  2:473496/150‑1 (MQ=255)                                                                    cAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTCAACCACAGCCTTGCGCTGTTTTCGGGGGGGGCGTTCGTTC‑TCGGGGTGGGCAg                                                                                     >  2:213660/1‑144 (MQ=255)                                                                      aCGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGACGGGTTGATGTTTAGCCCCAACCCCGCGCTGTTCTCGGGGGCGGGGTTCTTTC‑TCCgggggggccggcg                                                                                  >  1:861118/1‑138 (MQ=255)                                                                        gCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCCCAGCCTGGCGCTGTTTTCGGCGGGGATGTTTATTC‑TCGGGGGGGGCAGCGGCCTcccc                                                                           >  2:834362/1‑150 (MQ=255)                                                                         ccTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCAt                                                                         <  1:647321/151‑1 (MQ=255)                                                                            ggCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATg                                                                        <  2:69864/149‑1 (MQ=255)                                                                               tGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTTAGCCACAGCCTGGCGCTGTTCTCGGCGGGGATGTTCCTTT‑TCGGGGGGGGCAGCGGCATCACCATGTCg                                                                     >  2:884646/1‑149 (MQ=255)                                                                                gATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTTAGCGGCATCACCATGTCGAtt                                                                  <  2:1017410/151‑1 (MQ=255)                                                                                 aTGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTg                                                                 <  1:138853/151‑1 (MQ=255)                                                                                   ggAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTGTGTTTTGCCCCACCCTGGCGCTGGTCTCGGCGGGGGTGTTCCTTCTTCGGGGGGGAC‑GCGGCAACACCCTGTCGCTTg                                                                 >  1:511431/1‑149 (MQ=255)                                                                                       aTCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑t                                                          <  2:1066223/151‑1 (MQ=255)                                                                                         cGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGGGCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑tt                                                         <  1:1242079/150‑1 (MQ=255)                                                                                            cccGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTg                                                     <  2:496749/151‑1 (MQ=255)                                                                                               gTTGAAAACGCATTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTa                                                  <  1:274717/151‑1 (MQ=255)                                                                                                  gAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAaca                                               <  2:1017526/151‑1 (MQ=255)                                                                                                    aaaCGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTCAGCCCCATCCTGGGGCTGTTCTCGGGGGGGATGTTTATTT‑TCGGGGGGGTCAGCGGCATCACCCAGTCGCGTGGTGCA‑TTCCTGGGTacac                                              >  1:220526/1‑150 (MQ=255)                                                                                                    aaaCGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGCTGGCGGATTGAAGGTCAAACAAAACCCCGCGGTGTTCTCCGCGGCCCCATTCCCTC‑TCGGGGGGGGTAGCGGCGTCCCCCTGGCGAATGGAAAA‑TTTCCGGTGaca                                               >  2:824203/1‑149 (MQ=255)                                                                                                       cGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAAt                                          <  2:698656/151‑1 (MQ=255)                                                                                                          aGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTa                                       <  2:1231088/151‑1 (MQ=255)                                                                                                          aGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTCTGTTCAGCCCCAGCCTGGCGCCGTTCTCGGCGGCGACGGTCATTC‑TCGGGGGGGGCAGCGGCCTCACCCTGTCCATTGGTACA‑TTTCCGGTAACACAAATGt                                        >  1:1379109/1‑150 (MQ=255)                                                                                                            ttACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTATGTTCAGCCCCAGCCCGGCGCTGTTCTCGGCGGCGGTGTTCATTT‑TCGGGGGGGGCAGCGGCATCACCATGTCGGTTGGTACA‑TTTCTGGGAACACAAATGTA‑t                                     >  2:314090/1‑150 (MQ=255)                                                                                                              aCGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTCCGGGTTGATGTTCAGCCACAACCTGGCGCTGTTCTCGGCGGCGATGGTCATTC‑TCGGGGGGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGGAACACAAATGTA‑Tg                                    >  2:918234/1‑149 (MQ=255)                                                                                                               cGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAg                                 <  2:1352940/151‑1 (MQ=255)                                                                                                                      gCTTTCTCCTGATGGTGCTGGCGGTTTCCGGGTTGATGTTCCGCCACCCCCCGGCGCTGTTTTCGGCGGCGCTGTTCATTC‑TCGGGGGGGTCAGCGGGCTTCCCCTTTCGGTTGGTACA‑TTCCCGGGTAAACAAAAGGA‑TTAAggg                               >  1:807344/1‑146 (MQ=255)                                                                                                                      gCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCGTCa                          <  1:1217813/151‑1 (MQ=255)                                                                                                                      gCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTCAGCCCCAGCCTGGCGCTGTTTTCGGCGGCGGTGTTCATTC‑TCGGGGGGGTCAGCGGCATCACCCTGTCGGTTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCg                             >  2:334945/1‑148 (MQ=255)                                                                                                                        tttCTCCTGATGGTGCTGGCGGGTGCCCTGTTGATTTTTCCCACCCGCCTGGCGCTGTTCCCGGCGGTGATGGTTTTTT‑TCGGGGGGGTGCGCGGCCACACCATTTCCGTTGGGTCATTTCCTGGGAA‑ACAAAAGTATTGAAggcgcccgc                         >  1:511394/1‑144 (MQ=255)                                                                                                                          tctcCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGTGGTGGTCAGCGGCATCTCCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑GGAAGGGCTCCAgc                        <  1:854984/149‑1 (MQ=255)                                                                                                                             ccTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCGTCAGcgcg                     <  2:78477/149‑1 (MQ=255)                                                                                                                                 aTGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCGTCAGCGCGGTTccc               <  1:1316556/151‑1 (MQ=255)                                                                                                                                  tGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACAGCCTGGCGCTGTTCTCGGCGGCGGTGTTTCTTT‑TCGGGGGGGGCAGCGGCATCACCCTGTCGGTTGGGACA‑TTCCTGGGAAAACCAATGGA‑TGAAGGGCGTCCGCGCgggtgcc               >  1:17259/1‑145 (MQ=255)                                                                                                                                        tGGCGGTTGCCGGTTTGATGGTCAGCCCCCGCCTGGCGCTGTTCTCGGCGGGGGTGTTCATTC‑TCGGGGTGGTCAGCGGAATCACCATGTCGATTGGTACA‑TTTCTGGTAACACAAATGTA‑TGCAGGGCGTCAGCGCGGTTCCCGCCTGtt        >  1:882629/1‑151 (MQ=255)                                                                                                                                           cGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCGTCAGCGCGGTTCCCGCCTGTTAtt     <  2:152367/151‑1 (MQ=255)                                                                                                                                              ttGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAc  >  2:73266/1‑151 (MQ=255)                                                                                                                                              ttGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAc  >  2:356379/1‑151 (MQ=255)                                                                                                                                               tGCCGGTTTGATGTTCAGCCCCAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCGTCAGCGCGTTTCCCGCCTGTTATTTAc  <  1:833704/150‑1 (MQ=255)                                                                                                                                               tGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAc  <  2:1204622/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTC‑TCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTA‑TGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAC  >  NC_000913/3492686‑3492977

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