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

CGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA  >  NC_000913/3556061‑3556355                                                                                                                                                      |                                                                                                                                                     cGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTcccc                                                                                                                                                  <  1:713987/151‑1 (MQ=255)  gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCg                                                                                                                                                 <  1:578003/151‑1 (MQ=255)  gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCg                                                                                                                                                 >  2:894188/1‑151 (MQ=255)  gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCg                                                                                                                                                 >  2:311073/1‑151 (MQ=255)    aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGcc                                                                                                                                               <  2:488171/151‑1 (MQ=255)      ccaccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGcc                                                                                                                                               >  2:1191803/1‑149 (MQ=255)      ccaccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAg                                                                                                                                             <  1:365252/151‑1 (MQ=255)         ccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtgt                                                                                                                                          >  2:938598/1‑151 (MQ=255)               cACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGAt                                                                                                                                    >  2:368343/1‑151 (MQ=255)                 cTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGAttt                                                                                                                                  >  1:779354/1‑151 (MQ=255)                    aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCAc                                                                                                                               >  2:144906/1‑151 (MQ=255)                      gcgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGc                                                                                                                             >  2:1022300/1‑151 (MQ=255)                       cgcTTTCCACCGACAAAAAAAAAGCGTTGGGATATATTTTCACTTGCGCCTCCAAGCGAGGCGATATTCCCCGG‑CCCCT‑GGTCGCGCGGGAGGTTATAAATATTCTGTTCATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCa                                                                                                                                <  2:843122/147‑1 (MQ=255)                        gcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                           >  2:566725/1‑151 (MQ=255)                        gcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                           >  1:1043799/1‑151 (MQ=255)                          ttttCCCCGACAAAAAAAAAGCTTTGGGTAATTTTTTCACTTCCGCCTTAAAGCGAGCGGTTATTCCCCGG‑CCCCT‑GGTCGCGCGGCGGTTTATGAATTTTCTTTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                                          <  2:533515/150‑1 (MQ=255)                           tttCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAAGGTATTACCCGG‑CGCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCACGCCGGAGTAGCGATTTCACGCTCAGc                                                                                                                        >  1:819312/1‑151 (MQ=255)                              ccccGACAAAAAAAAAGCTTTGGGTGATTTTTTCTCTTTCGACTACAAGCGAGACGTTATTCCCCGG‑CCCGG‑GGTCGCGCGGGAGTTTATAAATTTTCTGTTCAGGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                                          <  2:116846/144‑1 (MQ=255)                               aCCGACGACAAAAAAGCGTTCGGCGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCaa                                                                                                                      >  1:2014/1‑149 (MQ=255)                                  gacgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                                 >  1:937867/1‑151 (MQ=255)                                   acaaaaaaaaaCGGTTGGGTCATTTTTTCCCTTCCGACTAAAAGGGAGGCGGTATTCCCGGC‑CCCCG‑GGTGGGGGGGCGGTTTATAATTTTGCGGTACAGGGGGGGAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                                   <  1:574053/143‑1 (MQ=255)                                      acAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                             >  2:956160/1‑151 (MQ=255)                                      acAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                             >  2:241041/1‑151 (MQ=255)                                        aaaaaaGCGTTGGGTGATTTTTTCCCTTTCGACTTAAAGCGAGACGGTATTCCCCGG‑CCCCT‑GGTCGCGCGGCAGTTTATGAATTTTCTTTTCATGCGGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                              <  2:1181588/148‑1 (MQ=255)                                          aaagGCGTTGGGTGATATTTCCCCTGCCGACTACAAGCGAG‑GGGTTATCCCCGG‑CCCCTGGGTCGCGCGGCAGGTTATGATTATTTTGTTCATGAAGAAAAACA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                                         <  1:1193144/151‑1 (MQ=255)                                           aaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                                         >  1:774747/1‑150 (MQ=255)                                            aaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                                       >  1:334120/1‑151 (MQ=255)                                                    ggTTATTTTTTCTTTTTCGCCTTAAAGGGAGACGTTATCCCCCGC‑CCCCG‑GGCCGCGGGGCAGGTTATGAATATTCTTTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTa                                                                                                 <  2:987663/149‑1 (MQ=255)                                                      tattatTTTCCCTTCCGCCTAAAAGCGAGCTGGTATTCCCCGG‑CCCGT‑GGTGGCGCGGCGGGTTATAAATTTTCTTTGCATGGAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                              <  1:872243/147‑1 (MQ=255)                                                         ttttttCCCTTTCGACTCAAAGCGGGGCGTTATCCCCCGC‑CCCTT‑GTCCCCGCGCCAGGTTATAAATATTCTGTACATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                             <  2:1051934/146‑1 (MQ=255)                                                            tttCCCTTTCTACTTCAGGCGAGACGGTATTCCCCGG‑CCCCG‑GGTCGCGCGGCAGGTTATGAATATTCTTTCCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                        <  2:947349/150‑1 (MQ=255)                                                            tttCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                                       >  1:533924/1‑151 (MQ=255)                                                            tttCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                                       >  1:1015320/1‑151 (MQ=255)                                                              tCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtctc                                                                                     >  1:313470/1‑151 (MQ=255)                                                               cctttccgcctaaaaGGGAGCCGTTTTTCCCCGC‑CCCCG‑GGTCCCGGGGCAGTTTATAAATTTTTTGTTCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                                      <  2:234319/136‑1 (MQ=255)                                                                 cTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                                  >  1:506752/1‑151 (MQ=255)                                                                  tttCAACTAAAAGCGAGGGGGTTTTCCCCGC‑CCCCT‑GGGCGGGGGGTAGTTTATAAATTTTCTTTTCA‑GGAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                                    <  2:1106733/148‑1 (MQ=255)                                                                   ttCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                                >  1:359554/1‑151 (MQ=255)                                                                   tccGACTACATGCGAGACGGTATTACCCGG‑CCCGT‑GGGCGCGCGGCAGGTTTTGAATTTTCTGTCCATCCGGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                                    <  2:609857/145‑1 (MQ=255)                                                                    tCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                                >  2:640698/1‑150 (MQ=255)                                                                     cgccTTCAAGCGAGACGGTATTCCCCGC‑CCCCT‑GGCGCCGCGGCAGTTTATGAATTTTCTTTTAAGGCGGGAAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                                  <  2:8232/144‑1 (MQ=255)                                                                      gACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATc                                                                             >  2:678895/1‑151 (MQ=255)                                                                        cTTAAAGGGAGACGTTATCCCCCGC‑CCCCT‑GGTCGCGCGGCATTTTATAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                            <  2:562787/150‑1 (MQ=255)                                                                         ttCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                          >  2:1071378/1‑151 (MQ=255)                                                                            aaGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                        >  2:102333/1‑150 (MQ=255)                                                                                  gaCGTTATTCCCCGC‑CCCCG‑GCTCGCGGGGCGGTTTATAAATTTTTTTTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                           <  1:76149/141‑1 (MQ=255)                                                                                    cGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                               >  1:1224164/1‑151 (MQ=255)                                                                                      gtttcccccGC‑CCCGT‑GGTGGCGGGGCAGGTTATAAATATCGTTTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTcc                                                                <  1:678644/142‑1 (MQ=255)                                                                                       tatcccccGG‑CCCCT‑GGCCGCGCGGCAGGTTATGAATATTCTTTCCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                              <  1:1053810/144‑1 (MQ=255)                                                                                       tATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgcg                                                            >  1:788023/1‑151 (MQ=255)                                                                                          tccccGGCCCCCG‑GGTCGCGCGGTGGTTTATAAATATTTTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                                  <  2:795119/141‑1 (MQ=255)                                                                                           ccccGC‑CCCCT‑GGGCGCGCGGCCGTTTATAAATATTCTGTTAATGAAGGAAAAAT‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                              <  2:333057/144‑1 (MQ=255)                                                                                           ccccGC‑CCCCG‑GGCCGCGCGGAAGGTTATGAATATTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg                                                        <  1:1197569/150‑1 (MQ=255)                                                                                            cccGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc                                                       >  1:823231/1‑151 (MQ=255)                                                                                            cccGC‑CCCCT‑GGCCGCGCGGCAGGTTAGAAATGTTTTTTAAATGGGGGAAAAAA‑CCCCCCGCCAGTGTAGAGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc                                                       <  2:958472/151‑1 (MQ=255)                                                                                                  ccccG‑GGCCGCGGGGAAGTTTATAAATATTCTTTTCATCCAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGt                                                  <  1:938513/151‑1 (MQ=255)                                                                                                   cccT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                                 >  1:52669/1‑151 (MQ=255)                                                                                                    ccT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTg                                                <  2:367062/151‑1 (MQ=255)                                                                                                        ggTCGCGGGGCAGTTTATAAATTTTTTTTTC‑TGCAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                             <  1:936742/151‑1 (MQ=255)                                                                                                         gTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                              >  2:1043505/1‑149 (MQ=255)                                                                                                           cgcgcgctggtttATAAATTTTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                             <  2:637059/137‑1 (MQ=255)                                                                                                                ggCAGTTTATGAATATTCTGTTCAGGGAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                       <  2:608910/149‑1 (MQ=255)                                                                                                                  aatgTTATGAATTTTCTTTCCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                   <  1:1238848/148‑1 (MQ=255)                                                                                                                           aaTTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAc                           <  2:1183415/150‑1 (MQ=255)                                                                                                                           aaTTTTCTGTCCACGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                               <  2:672266/146‑1 (MQ=255)                                                                                                                           aaTATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                          >  1:546251/1‑151 (MQ=255)                                                                                                                               ttCTATGCAGGCAGGGAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc                       <  2:869553/150‑1 (MQ=255)                                                                                                                                     tCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACtt                >  2:1215162/1‑151 (MQ=255)                                                                                                                                     tCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACtt                >  1:1240053/1‑151 (MQ=255)                                                                                                                                       aTGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc              >  1:316167/1‑151 (MQ=255)                                                                                                                                        aggggggAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc               <  2:217559/146‑1 (MQ=255)                                                                                                                                           aGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt            >  1:804641/1‑149 (MQ=255)                                                                                                                                            gaaaaaac‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc              <  1:1022019/143‑1 (MQ=255)                                                                                                                                                    ccTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACcaca  <  1:1151986/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA  >  NC_000913/3556061‑3556355

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