Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	44.0%	98.9 / 30.2	75	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (33/9);  new base C (0/33);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 2.50e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTG  >  NC_000913/3556060‑3556347                                                                                                                                                        |                                                                                                                                             ccGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTccc                                                                                                                                           >  1:934583/1‑151 (MQ=255)   gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCg                                                                                                                                         <  2:706095/151‑1 (MQ=255)   gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCg                                                                                                                                         <  2:1121177/151‑1 (MQ=255)     aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGcc                                                                                                                                       <  2:831601/151‑1 (MQ=255)       ccaccTCGGCACTGACGCGCTTTTCCCCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCa                                                                                                                                      <  2:607243/150‑1 (MQ=255)        caccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCCCTTTCGA‑CTGCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAg                                                                                                                                     <  1:143858/150‑1 (MQ=255)         accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGC‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAgtg                                                                                                                                   >  2:369592/1‑151 (MQ=255)            tCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTCCCCTTTCGA‑CTACAAGCGAGACGTTATT‑CCCCGCCCCT‑GGTCCGCGCGGCAGTTTATGAATTTTCTGTTCATG‑AAGGAAAAACCCCCCCGCCAgtgt                                                                                                                                  <  1:601476/149‑1 (MQ=255)              ggCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCg                                                                                                                              >  1:408769/1‑151 (MQ=255)                   tgccgcgTTTTCCCCCAACGAAAAAAAAGCGTTCGTTGATATTTTCTCTTCCGA‑CTTAAAGGGAGACGTTATTCCCCGGCCCCG‑GGT‑CGCGGGGCAGTTTATGAATTTTTTGTCCATG‑AAGGAAAAACCCCCCCGCCAGTGTAGCGAttt                                                                                                                          <  2:467547/147‑1 (MQ=255)                       gcgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGc                                                                                                                     >  2:341959/1‑151 (MQ=255)                        cgcTTTCCCCCGACGACAAAAAAGCGTTGGGTGATTTTTTCTCTTTCAA‑CTCAAAGCGAGACGTTATCCCCCGGCCCCG‑GGG‑CGCGCGGCAGGTTATGAATATTCTGTTCAGG‑AAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACg                                                                                                                      <  1:535500/149‑1 (MQ=255)                            tttCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                                >  1:963009/1‑151 (MQ=255)                              tCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCa                                                                                                               >  2:175832/1‑150 (MQ=255)                               cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                             >  1:153168/1‑151 (MQ=255)                                 ccaaggaaaaaaaaGCTTTGGATAATTTTTTCA‑TTCCGATCTAAAACCGAGCCGGTTTTCCCCGGCCCCT‑GGC‑CGCGCGGCGGGTTATGAATTTTCTTTTCAAG‑AAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                           <  2:220876/143‑1 (MQ=255)                                 ccGACGAAAAAAAAGCGTTGGGTGATATTTCCACTTCCGC‑CTTCAAGCGAGCCGGTATTCCCCGCCCCCT‑GGT‑CGCGCGGCAGGTTATAAATTTTCTGTTCATG‑CAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                           <  2:1277152/151‑1 (MQ=255)                                    acgaaAAAAACGCGTTGGGTTATATTTTCACTTCCGC‑CTACAAGCGAGACGTTATTCCCCGCCCCCG‑GGT‑CGCGCGGGAGTTTAGAATTTTTCTGTTCATG‑GAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAg                                                                                                                 <  1:685054/142‑1 (MQ=255)                                     cgacAAAAAGGCGTTCGGTGATATTTTCCCTTTCGA‑CTTCAAGCGAGCCGTTATTCCCCGGCCCCT‑GGC‑CGCGCGGCAGGTTATAAATATTCTGTTAATG‑AAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                         <  1:735081/149‑1 (MQ=255)                                     cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg                                                                                                       >  1:1245582/1‑151 (MQ=255)                                     cgaaAAAAAGGCGTTGGGTGATTTTTTCCCTTGCGC‑CTTCAGGCGAGCCGGTATTCCCGGGCCCCT‑GGT‑CGCGCGGCAGTTTATAAATTTTCTTTTCATG‑CAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                          <  2:864518/148‑1 (MQ=255)                                       acAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                     >  1:1001630/1‑151 (MQ=255)                                       aaaaaaaaGCGTTCGGTGATATTTTCCCTTTAGA‑CTTCAAGCGAGCGGGTATTCCCCGGCCCCT‑GGT‑GGCGCGGCAGGTTATGAATATTTTTTTCATG‑CAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                           <  2:1014921/143‑1 (MQ=255)                                        cAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                    >  1:170456/1‑151 (MQ=255)                                        aaaaaagGCGTTGGGTAATTTTTTCTTTTTCAC‑CTTCAAGCGAGACGGTATCCCCCGGCCCTG‑GGC‑CGCGGGGGGGTTTATAAATTTTTTGTTCATG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                    <  1:922692/150‑1 (MQ=255)                                              aGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                              >  1:347770/1‑151 (MQ=255)                                                 tttCGGTTATTTTTTCGCTGTCGA‑TGTCAAGCAGGTCTTTATTCCCCGCCCCCG‑GGC‑CGCGCGGGAGTTTACAAGTTTTTTTTTCAGT‑AAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                             <  1:1295070/148‑1 (MQ=255)                                                  ttgggTAATTTTTCCCCTTTCGC‑CTTAAAGGGAGGGGTTATCCCCCGCCCCCG‑GGC‑CGCGCGGCGGTTTATAAATTTTCTGTTCAGG‑CGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAg                                                                                            <  2:15948/146‑1 (MQ=255)                                                     ggTTATTTTTTCACTTTCGA‑CTTCAAGCGAGACGGTATCCCCCGGCCCCT‑GGT‑CGCGCGGCAGTTTATGAATTTTCTGTTCAGG‑GAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAg                                                                                            <  1:1028668/146‑1 (MQ=255)                                                     ggTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                        >  1:715523/1‑150 (MQ=255)                                                            ttttCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                >  2:1067444/1‑151 (MQ=255)                                                             tttCTCTTGCGA‑CTCCAAGCGAGACGGTATCCCCCGGCCCCT‑GGT‑CGCGCGGCGGTTTATGAATATTCTGTTCATG‑CAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                <  2:519518/150‑1 (MQ=255)                                                              ttCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                              >  2:1021435/1‑151 (MQ=255)                                                               tCACTTTCGA‑CTGCAAGCGAGACGGTATTCCCCGGCCCCT‑GTG‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                               <  1:407148/149‑1 (MQ=255)                                                                    ttgctacttaaacagagaCGTTATCCCCCGCCCCCG‑GGC‑CCGGCGGCGGGTTATAAATTTTCTGTCAATG‑AAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                           <  1:271417/135‑1 (MQ=255)                                                                          ctaaaaGCGAGACGTTATTCCCCGGCCCCG‑GGT‑CGCGCGGCAGTTTATGAATATTCTTTTCATG‑CAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                    <  2:75999/146‑1 (MQ=255)                                                                          cTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                   >  1:737000/1‑151 (MQ=255)                                                                           ttaaaGCGAGACGTTTTTCCCCGGCCCCT‑GGT‑CGCGCGGCAGTTTATAAATATTCTGTTCAGG‑CGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                    <  2:1224112/146‑1 (MQ=255)                                                                              aggCGAGACGGTATCCCCCGCCCCCG‑GGT‑CGCGCGGCGGTTTATAAATATTCTGTTCATG‑CGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                    <  2:729615/144‑1 (MQ=255)                                                                              aaGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                >  2:616801/1‑150 (MQ=255)                                                                                gggggaggttATTCCCCGGCCCCT‑GGT‑CGCGGGGCAGGTTATAAATTTTTTGTTCATG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaa                                                              <  2:519951/141‑1 (MQ=255)                                                                                  gagaCGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                           >  2:255494/1‑151 (MQ=255)                                                                                   agaCGTTTTTCCCCGGCCCCG‑GGC‑CGCGCGGCAGGTTATGAATATTTTTTTCAGG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATg                                                            <  2:497720/149‑1 (MQ=255)                                                                                    gccgttATTCCCCGCCCCCT‑GGT‑CGCGCGGCAGGTTATAAATATTCTGTTCAGG‑CAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                          <  2:231121/145‑1 (MQ=255)                                                                                    gaCGTTTTTCCCCGCCCCCT‑GGT‑CCCGCGTCAGGTTATAAATATTCTGTTAATG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                          <  2:873671/150‑1 (MQ=255)                                                                                    gaCGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTc                                                         >  2:274232/1‑151 (MQ=255)                                                                                      cGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                       >  1:1221339/1‑151 (MQ=255)                                                                                      cGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                       >  1:1065263/1‑151 (MQ=255)                                                                                       ggTATTCCCCGCCCCCG‑GGC‑CGCGGGGCAGGTTATGAATTTTCTGTTCAGG‑AAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTc                                                         <  2:679377/148‑1 (MQ=255)                                                                                       ggTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                      >  1:1048193/1‑151 (MQ=255)                                                                                       ggTATCCCCCGCCCCCG‑GGT‑CGCGCGGCAGTTTATAAATATTCTGTTCATG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                      <  2:1028415/151‑1 (MQ=255)                                                                                        gTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                     >  2:643488/1‑151 (MQ=255)                                                                                         tATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgcg                                                    >  1:1089408/1‑151 (MQ=255)                                                                                          tttcccccGCCCCTGGGT‑CGCGCGGGAGTTTATAAATATTCTTTTCATG‑CGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATg                                                            <  2:82596/139‑1 (MQ=255)                                                                                           ttccccGGCCCCT‑GGT‑CGCGCGGCAGGTTATAAATTTTCTTTTCATG‑CAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCgcgc                                                   <  2:909693/147‑1 (MQ=255)                                                                                           ttccccGCCCCCG‑GGT‑CGCGGGGCAGGTTATAAATTTTCTGTTCATG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt                                                             <  2:1124133/137‑1 (MQ=255)                                                                                              cccGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc                                               >  1:431388/1‑151 (MQ=255)                                                                                               ccGCCCCGT‑GGT‑CGCGCGGCGGGTTATGAATTTTCTGTTCATG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                              <  1:1176528/151‑1 (MQ=255)                                                                                                  cccccG‑GGT‑CGGGCGGGAGGTTATGAATATTCTGTACAGG‑AAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                           <  2:434579/150‑1 (MQ=255)                                                                                                             cgcgcgGCAGTTTATGAATTTTCTTTTCATG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                     <  1:582702/148‑1 (MQ=255)                                                                                                               cgcgGCAGGTTATGAATATTCTTTTCATG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAgg                                   <  2:483193/148‑1 (MQ=255)                                                                                                               cgcgGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtgt                                >  2:1000569/1‑151 (MQ=255)                                                                                                                     aGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                          >  2:363913/1‑151 (MQ=255)                                                                                                                      ggTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                          >  1:1039536/1‑150 (MQ=255)                                                                                                                        ttATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                       >  2:360235/1‑151 (MQ=255)                                                                                                                          aTGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACCATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                     >  1:34295/1‑151 (MQ=255)                                                                                                                             aaTATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                  >  1:575987/1‑151 (MQ=255)                                                                                                                             aaTATTCTGTCCATG‑CGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACGTGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                  <  1:251724/151‑1 (MQ=255)                                                                                                                              atatTCT‑TTCAAGCCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                    <  1:848525/148‑1 (MQ=255)                                                                                                                                    tgtaaaTG‑CAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc               <  2:686957/142‑1 (MQ=255)                                                                                                                                    tGTTCAGG‑AAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag            <  2:1541/150‑1 (MQ=255)                                                                                                                                        cATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc       >  2:94152/1‑151 (MQ=255)                                                                                                                                             cAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGCGTTAAACGATGCCACGGGCGAGACTTCCGtt   >  2:360415/1‑151 (MQ=255)                                                                                                                                              aGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTg  >  1:746887/1‑151 (MQ=255)                                                                                                                                              aGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACGCCAGCTAATAGAACTTCTCCACGCAGCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCg     <  1:138426/148‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                             CCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGA‑CTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATG‑CAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTG  >  NC_000913/3556060‑3556347

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