Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	38.6%	91.6 / 13.7	57	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (31/4);  new base C (0/22);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.28e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACC  >  NC_000913/3556059‑3556357                                                                                                                                                        |                                                                                                                                                       gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTcc                                                                                                                                                      <  1:717638/151‑1 (MQ=255)     cAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGc                                                                                                                                                  <  2:691102/151‑1 (MQ=255)      aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCGACTACAAGCGAGACGGTATTACCCGGCCCCTGGG‑CGCGCGGAAGGGTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCg                                                                                                                                                   >  1:104164/1‑149 (MQ=255)          accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtg                                                                                                                                             >  1:977692/1‑151 (MQ=255)           ccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtgt                                                                                                                                            >  2:298799/1‑151 (MQ=255)               ggCACTGCGGCGCTTTCCACCGACAAAAAAAAGGGGTTCGG‑TGATTTTTGCACTGTCGACTACAAGGGAGACGTTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTAGAAATATTCTTTTCATGCAGGGAAAAC‑CTCCCCGCCAGTGTAg                                                                                                                                          <  2:747341/149‑1 (MQ=255)                 cACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGAt                                                                                                                                      >  2:216805/1‑151 (MQ=255)                          gcTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                             >  1:791935/1‑151 (MQ=255)                            ttttCCCCGCCGAAAAAAAAGCGTTGGG‑TTATTTTTTCACTTCCGACTTCAAGCGAGCCGGTATTCCCCGCCCCCGGGT‑CGCGCGGCAGTTTATAAATTTTCTGTTC‑TTCAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                              <  2:1035289/148‑1 (MQ=255)                            ttttCACCGACGACAAAAACGAGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCATGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                             >  2:36940/1‑149 (MQ=255)                            tttcCACCGACGACAAAAAAGCGTTGGG‑TGATATTTTCCCTTCCGACTTCAAGCGAGAGGTTATTCCCCGCCCCCTGGT‑CGCGCGGCAGGTTATAAATATTCTTTTAATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACg                                                                                                                                <  1:80184/146‑1 (MQ=255)                            tttcCACCGACGAAAAAAAAGCGTTGGG‑TGATTTTTTCCCTTCCGCCTTAAAGCGAGACGTTATTCCCCGCCCCCTGGT‑CGCGCGGCAGGTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                             <  1:367331/149‑1 (MQ=255)                                     acgacAAAAAAGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                                  >  1:451003/1‑151 (MQ=255)                                        aaaaaaagGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGCCCCCGGGC‑CGCGCGGCAGTTTATAAATATTCTTTTCATGCAGGGAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                               <  1:470506/149‑1 (MQ=255)                                        aaaaaaaaGCGTTGGG‑TAATATTTTCACTTCCGACTCAAAGGGAGACGGTTTTCCCCGCCCCCGGGT‑CGCGCGGCAGGTTAGAAATTTTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                                    <  1:658459/144‑1 (MQ=255)                                           aaaaaGC‑TTCGGTTAATATTTTCCTTTTCAACTACAGGTGAAACGGTATT‑CCCGGCCCCCGGTGCGGGCGACATTTTATAAATATTTTGTTCTCGCGGGAAAAAT‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                            <  1:91795/151‑1 (MQ=255)                                            aaaaGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGCCCCCTGGC‑CGCGCGGCGGGTTATGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                                           <  1:802924/151‑1 (MQ=255)                                                   ttCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGG‑AAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                                    >  1:331108/1‑150 (MQ=255)                                                       gttGATTTTTCCTCTTTCGTCTACAAGCGAGACGGTTTTCCCCGGCCCCGGGG‑CGCGCGGCAGTTTTTGAATATTCGTTGCATGCGGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTa                                                                                                   <  2:638648/147‑1 (MQ=255)                                                              ttttCCCTTTCGCCTCCAAGGGAGCCGGTATTCCCCGGCCCCGGGT‑CGCGCGGCGGTTTATGAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                            <  1:930731/149‑1 (MQ=255)                                                                ttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                                        >  1:720457/1‑151 (MQ=255)                                                                 tCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtctc                                                                                       >  1:844476/1‑151 (MQ=255)                                                                 tCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑ATCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                                        >  2:52172/1‑150 (MQ=255)                                                                  cccTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCGGGT‑CGCGCGGCGGTTTATGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                                        <  1:871838/147‑1 (MQ=255)                                                                  cACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                                      >  1:134850/1‑151 (MQ=255)                                                                   aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                                      >  1:1021763/1‑150 (MQ=255)                                                                     tttAGACTTCAAGCCAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGCTATGAATATTCTGATCATGCACGGAACAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCCAATAGAACTTCTCCACg                                                                                   >  2:574234/1‑151 (MQ=255)                                                                      gtCGACTACAAGGGAGGCGGTTTCACCCGGCCCCGGGG‑CGCGCGGCAGGTTATGAATATTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                                   <  2:656918/149‑1 (MQ=255)                                                                         gACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATc                                                                               >  1:217380/1‑151 (MQ=255)                                                                            ttCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                            >  2:118913/1‑151 (MQ=255)                                                                              cAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                          >  1:169766/1‑151 (MQ=255)                                                                               aaGCGGGCCGGTTTTCCCCGCCCCCTGGC‑CGCGCGGCGGTTTATGAATATTCTGTCCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                            <  1:300937/148‑1 (MQ=255)                                                                               aaGCGGGAGGTTATCCCCCGGCCCTGGGG‑CGCGCGGCAGGTTATAAATATTCTGTCCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                              <  1:463315/146‑1 (MQ=255)                                                                               aaGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                          >  2:680826/1‑150 (MQ=255)                                                                                    agaCGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATGTTCTGATCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                                     >  2:589484/1‑150 (MQ=255)                                                                                         gTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                               >  1:989404/1‑151 (MQ=255)                                                                                                      ccTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTg                                                  >  1:121214/1‑151 (MQ=255)                                                                                                                ggggCAGTTTATAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtg                                           <  1:962821/147‑1 (MQ=255)                                                                                                                gcggccggtagagAATTTTCGTTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACGTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                                              <  2:137210/134‑1 (MQ=255)                                                                                                                gcgGCAGGTTATAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                         <  1:161402/151‑1 (MQ=255)                                                                                                                  ggCAGTTTATGAATTTTCTTTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                         <  1:966785/149‑1 (MQ=255)                                                                                                                           tGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                               >  1:103086/1‑150 (MQ=255)                                                                                                                           tGAATATTCTGTTCATCCAGGGAAAAA‑CTCCCCGCCAGAGTAGCGATTTCACGCTCAGCAACATGGCGCGCCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAGGTTCCCGCGCTCGCCAAGTGGCAAGGTGTTAAACCATGcc                               >  2:307000/1‑150 (MQ=255)                                                                                                                           tGAATATGCGGTTAATGCAGAAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                                   <  2:441400/146‑1 (MQ=255)                                                                                                                            gAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAc                             >  1:992228/1‑151 (MQ=255)                                                                                                                                 ttCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg                          <  1:59396/149‑1 (MQ=255)                                                                                                                                 tccTGTTCA‑GGAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                              <  2:120957/143‑1 (MQ=255)                                                                                                                                  tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg                        >  1:747598/1‑150 (MQ=255)                                                                                                                                  tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg                        >  1:31988/1‑150 (MQ=255)                                                                                                                                  tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg                        >  1:348747/1‑150 (MQ=255)                                                                                                                                     gTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc                    >  1:504319/1‑151 (MQ=255)                                                                                                                                         caggcgggAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATAGCGCGGCACACCAGCTAAGAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTCCAAGGTGTTAAACGATGCCACGGGCGAGACt                   <  2:102222/144‑1 (MQ=255)                                                                                                                                           gCAGAGAAAAA‑CTCCCCGCCAGTTTAGCGCTTTCACCCTCATCAACATAGCGCGGCACACCAGCTAATAGCACTTCCCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGACTCTCCGGGCGCGACTTCCGt              >  2:1051264/1‑151 (MQ=255)                                                                                                                                              ggaaaaac‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGc           <  1:794027/148‑1 (MQ=255)                                                                                                                                                 aaaac‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGcc          <  1:522389/149‑1 (MQ=255)                                                                                                                                                    aa‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCcaccac     >  1:790691/1‑151 (MQ=255)                                                                                                                                                       cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCacac   >  2:426779/1‑151 (MQ=255)                                                                                                                                                        tCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACAcc  >  1:559842/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACC  >  NC_000913/3556059‑3556357

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