Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,257,698	0	T	G	41.2%	122.3 / 21.5	80	V195G (GTT→GGT)	lexA	transcriptional repressor of SOS regulon
Reads supporting (aligned to +/- strand):  ref base T (10/37);  new base G (33/0);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 6.08e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCC  >  NC_000913/4257553‑4257846                                                                                                                                                  |                                                                                                                                                        gCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTc                                                                                                                                                    <  2:949997/151‑1 (MQ=255) gCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTa                                                                                                                                                       >  2:230596/1‑148 (MQ=255)        ttGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACgg                                                                                                                                             >  2:823684/1‑151 (MQ=255)         tgatgaCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTTGCAACgg                                                                                                                                             >  2:458095/1‑150 (MQ=255)            tgaCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGc                                                                                                                                            >  2:107337/1‑148 (MQ=255)             gaCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACt                                                                                                                                        <  1:1172303/151‑1 (MQ=255)             gaCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACt                                                                                                                                        <  2:210355/151‑1 (MQ=255)               cGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGGCCTTCGTCAGCAGAGCTTCTCCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTgg                                                                                                                                      <  1:751211/151‑1 (MQ=255)                 aaGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCt                                                                                                                                    <  2:424050/151‑1 (MQ=255)                 aaGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCt                                                                                                                                    <  1:822980/151‑1 (MQ=255)                  aGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTTTTTGCAACGGCGACTgg                                                                                                                                      >  1:888410/1‑148 (MQ=255)                    ttACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTg                                                                                                                                   >  2:927223/1‑149 (MQ=255)                     tACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTa                                                                                                                                 <  2:926653/150‑1 (MQ=255)                      aCCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAAcagaggatggccgtt                                                                                                                                 >  2:1168878/1‑135 (MQ=255)                        cGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTAATCGCCACGGGGACTGGCtttta                                                                                                                                >  1:388362/1‑144 (MQ=255)                            aaGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATAtct                                                                                                                       <  1:854489/151‑1 (MQ=255)                            aaGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAAAGGCGACTGGGCGTa‑aaa                                                                                                                             >  2:409269/1‑144 (MQ=255)                                  cTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCCACGGCGACTGGCTGTA‑AC‑ATATctct                                                                                                                     >  2:696567/1‑147 (MQ=255)                                  cTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATctct                                                                                                                     >  1:142782/1‑147 (MQ=255)                                   tGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGATGGGGGTATTCGCAACGGCGACTGGCTGTA‑AA‑ATAtcctt                                                                                                                     >  1:366646/1‑143 (MQ=255)                                    gAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AA‑ATATCTCTGAGAcc                                                                                                               >  2:1225535/1‑151 (MQ=255)                                      aaaaaCAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGGGACcgc                                                                                                             >  2:784843/1‑151 (MQ=255)                                          aCAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATg                                                                                                         <  2:654028/151‑1 (MQ=255)                                          aCAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATg                                                                                                         <  2:531648/151‑1 (MQ=255)                                          aCAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGAATGGCTGTA‑AC‑ATATCTCTGAGACCGCGATg                                                                                                         >  2:693543/1‑151 (MQ=255)                                           cAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGGTGGGGGTTATCGCCACGGGCGCTGGCCGTAGACGATTTCTCTgagt                                                                                                                 >  1:578360/1‑143 (MQ=255)                                            aGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATgcc                                                                                                       <  2:913698/151‑1 (MQ=255)                                                   tAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGGGGTATTCGCTACGGCGACTGGCCGTA‑AC‑CAATTTCTTAGACcgcgttaccccc                                                                                                    >  2:1155824/1‑139 (MQ=255)                                                     aaGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTTTTTGCAACGGCGGCTTGCTTTT‑AC‑ATATCTCTGGGGGCGCGGTTCCGCCTGGCgg                                                                                              >  1:414452/1‑150 (MQ=255)                                                      aGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTc                                                                                             <  1:666534/151‑1 (MQ=255)                                                      aGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTc                                                                                             <  1:556941/151‑1 (MQ=255)                                                      aGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTc                                                                                             <  2:833884/151‑1 (MQ=255)                                                       gTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTcg                                                                                            <  2:699775/151‑1 (MQ=255)                                                       gTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTcg                                                                                            <  1:1233431/151‑1 (MQ=255)                                                        tCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGGGACCGCGATGCCGCCCGgcggcg                                                                                            >  1:27507/1‑150 (MQ=255)                                                         cGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTcgcg                                                                                          <  2:879111/151‑1 (MQ=255)                                                         cGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTcgcg                                                                                          <  2:730000/151‑1 (MQ=255)                                                          gAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTTTTTGCCACGGCGGCTGGCCGTG‑AC‑ATATCTCTGAGACCGCGGTGCCGCCCGgcggcg                                                                                            >  1:675278/1‑148 (MQ=255)                                                           aaCTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGGGGTTTTCGCAAAGGCGACTGGCTGTA‑AC‑ATATCTCTGGGACCCCGATGCCGCCTGGCGTCGCgg                                                                                         >  2:1080148/1‑150 (MQ=255)                                                           aaCTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCCACGGGGAGTGGCTGTA‑AA‑AAATCTCTGAGATCGCGATGCCGGCTGGGGtggggg                                                                                         >  2:523275/1‑145 (MQ=255)                                                             cTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCCACCGCGACTGGCTGTA‑AA‑ATATCTCTTAGACCCCGATGCCGCCTGGCGGCGCGggtt                                                                                      >  1:84173/1‑151 (MQ=255)                                                                ttGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCCCAACGGCGACCGGCTGTC‑AC‑ATATTTCTGTGACCGCGGTTCCGCCTGGGGTCGCGggttgt                                                                                    >  1:968361/1‑150 (MQ=255)                                                                        aaaTAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTA‑AA‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGGTTGTTTTTCAtc                                                                            >  1:911373/1‑150 (MQ=255)                                                                         aaTAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTTGCAACGGCGACTGGCTGTA‑AC‑ATATTTCTGGGACCGCGATGCCGCCCGGCGGCGCGGGTTGttttt                                                                                >  1:117438/1‑145 (MQ=255)                                                                           tAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTTGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGGCGCCCGGCGCCGCGGGTTGTTTTTCAt                                                                             >  1:1078235/1‑146 (MQ=255)                                                                           tAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACCGCGCCTGGCTGTA‑AA‑ATCTCTCTTAGACCGCGATTCCGCCCGGGGGCGCGGGTTGttttt                                                                                >  1:1229411/1‑143 (MQ=255)                                                                             gCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCAtctct                                                                         <  1:1130964/148‑1 (MQ=255)                                                                                aGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTtcat                                                                   <  2:479632/151‑1 (MQ=255)                                                                                aGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTA‑AA‑ATATCTCTGAGACCGCGATGCCGCCCGGCGTCGCGGGTTGTTTTTCATCTCTCtt                                                                      >  1:492780/1‑148 (MQ=255)                                                                                 gTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGGCGCGGGTTGTTTTTCATCTCTCtt                                                                      >  2:898280/1‑147 (MQ=255)                                                                                   ttAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTTGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCCGGCGTCGCGGGTTGTTTTTTCTCTCTCtt                                                                      >  1:1233093/1‑145 (MQ=255)                                                                                   ttAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGGTGGGGGTATTCGCAACGGCGGCTGGCTGTG‑AC‑ATATCTCTGGGACCGCGATCACGCCCGGGGTCGCGGGTTGGTTTTTAtctttttttat                                                                   >  1:705652/1‑140 (MQ=255)                                                                                    tAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTTGCAAAGGCGGCTGGCTGTT‑AA‑ATATCTCTGAGGCCGCGGTGCCGGCCGGGGTCGCGGGTTGTTTTTCCTCTCTCTTCCTCAgg                                                               >  1:867272/1‑151 (MQ=255)                                                                                     aaaCCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCCGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCCGGCGTCGCGGGTTGTTTTTTATCTCTCTtcatcacg                                                               >  1:98135/1‑148 (MQ=255)                                                                                       aCCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCtt                                                            <  2:1049457/151‑1 (MQ=255)                                                                                       aCCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGGCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCtt                                                            >  2:39111/1‑151 (MQ=255)                                                                                        ccAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTg                                                           <  1:226426/151‑1 (MQ=255)                                                                                          aaTTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCCGGCGGCGCGGTTTGTTTTTTATCTCTCTTCATCAGGCTTGt                                                          >  1:840000/1‑150 (MQ=255)                                                                                              gTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGc                                                      <  1:605239/150‑1 (MQ=255)                                                                                                cGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGGTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCa                                                     >  2:1085053/1‑149 (MQ=255)                                                                                                 gTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTTTTCGCAACCGCGAATGGCTGTA‑AC‑ATATCTCTGGGACCCCGCTACCGCCTCGCGTCGCGGGTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTTCAttg                                                  >  1:938043/1‑149 (MQ=255)                                                                                                 gTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATgg                                                  <  1:676644/151‑1 (MQ=255)                                                                                                   tGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGc                                                 <  2:825112/150‑1 (MQ=255)                                                                                                    gACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTTATCTCTCTTCATCAGGCTTGTCTGCATGGC‑a                                               >  2:462806/1‑150 (MQ=255)                                                                                                     aCCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTTGCAACGGCGGCCGGCTGTA‑AC‑ATATCTCTGAGGCCGGGATGCCGCCCGGCGTCGCGGGTTGTTTTTTATCTCTCTTTCTCAGGCCTTTCTGCATGGC‑Att                                             >  1:589046/1‑151 (MQ=255)                                                                                                          cGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGAC‑ATTCCTc                                         <  2:441113/150‑1 (MQ=255)                                                                                                             cagcagAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACt                                      <  1:1225471/150‑1 (MQ=255)                                                                                                                   agCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACCGGCTGGA‑AA‑ATTTTTCTGAGACCCCGATGCCGCCCGGCGGCGCGGGTTGTTTTTTATTTCTCTTTTTTAGGCTTTTCTGCATGGCAATTCCTccttcct                                   >  1:1060288/1‑142 (MQ=255)                                                                                                                        cACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑TC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATaa                           <  1:246869/150‑1 (MQ=255)                                                                                                                        cACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATaa                           <  2:933188/150‑1 (MQ=255)                                                                                                                         aCCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGGCGCGGGTTGTTTTTCATTTCTCTTCATCAGGCTTGTCTGCATGGC‑CTTTCCCACTTCATCTGATaaa                          >  2:1030232/1‑150 (MQ=255)                                                                                                                           cATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAg                         <  2:423069/149‑1 (MQ=255)                                                                                                                              tGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAGCAct                     <  2:651372/150‑1 (MQ=255)                                                                                                                               gAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAGCActc                    <  1:198191/150‑1 (MQ=255)                                                                                                                                aaGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAGCACTCTg                  <  2:772320/151‑1 (MQ=255)                                                                                                                                       ggCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAGCACTCTGGCAtc             <  1:458084/149‑1 (MQ=255)                                                                                                                                          ggTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGcc        <  1:39108/151‑1 (MQ=255)                                                                                                                                              ggggTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTa     <  2:383172/150‑1 (MQ=255)                                                                                                                                              ggggTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTAc    >  1:1000374/1‑151 (MQ=255)                                                                                                                                                 gTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTAccc  <  1:689599/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                        GCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTA‑AC‑ATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGC‑ATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCC  >  NC_000913/4257553‑4257846

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