Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,159	0	G	A	45.3%	69.1 / 35.5	64	intergenic (+139/‑75)	thrC/yaaX	L‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (31/4);  new base A (0/29);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 2.98e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCG  >  NC_000913/5010‑5305                                                                                                                                                         |                                                                                                                                                   atcatcAGTAAAATCTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgaga                                                                                                                                                   <  1:248051/151‑1 (MQ=255)    atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCTGCTTTTTTCTGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGAGGTAACTTAGAGATTAGGATTGCGgaga                                                                                                                                                   >  1:709740/1‑148 (MQ=255)    atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                 >  2:363009/1‑150 (MQ=255)      caGTAAAATCTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                              >  2:80939/1‑151 (MQ=255)        gTAAAATCTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaaca                                                                                                                                             >  1:937803/1‑150 (MQ=255)          aaaaTCTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                           <  1:649072/150‑1 (MQ=255)          aaaaTCTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                          >  2:464915/1‑151 (MQ=255)           aaaTCTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                         >  1:531764/1‑151 (MQ=255)               cTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGt                                                                                                                                     >  2:832064/1‑151 (MQ=255)                     aTTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCAt                                                                                                                               >  2:943675/1‑151 (MQ=255)                       tATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                             >  1:4709/1‑151 (MQ=255)                        aTCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                             >  2:959142/1‑150 (MQ=255)                         tctcAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                             >  2:228383/1‑149 (MQ=255)                         tctaAATCAGGCGGGT‑TTGGCTTTTAGGCACCCCGGTTTTTTTAAGAAAAAATTAGGGAGAAAAAGGACGGGAAAAAAGG‑AGAAA‑TTCTCAATAAGTGGGTTACTTTAGAGATTAGGATGGGGGAAAAAAACACCCGCCGTtctc                                                                                                                                 <  1:851643/145‑1 (MQ=255)                             aaTCAGGCGGGG‑TTTGTTTTTATGCAGCCCGGTTTTTTTATAAAAAAATTATGGAAAAAAATGACAGAAAAAAAGA‑AAAAA‑TTCTCAAAAAATGCGTTAATTTAGAGATTAGGATGGCGGAGAAAAACAACCGCCGTTCTCATCGAGTAAt                                                                                                                       <  1:643154/151‑1 (MQ=255)                             aaTCAGGCGGGG‑TTTGTTTTTATGCAGCCCGGTTTTTTTATAAAAAAATTATGGAAAAAAAGGACAGGAAAAAAGG‑AAAAA‑TTCAAAAAAAAGGCGGTAACTAAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAAt                                                                                                                       <  2:134344/151‑1 (MQ=255)                                   gCCGGGTTTTGTTTTTATGCACCCCGTTTTTTTTATAAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGA‑AAAAA‑TCCTAAATAAATGCGGTAACTAAGAGATTGGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                     <  2:1289/148‑1 (MQ=255)                                      gggttttttttttATGCACCCCGGTTTTTTTAAAAAAAAATTAGGAAAAAAAATGAGGGGAAAAAAGA‑AAAAA‑TTCTAAAAAAATGGGTTAACTTAGATATAAGGTTGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                 <  2:31581/140‑1 (MQ=255)                                      ggg‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                               >  2:742347/1‑150 (MQ=255)                                               ttttATGCAGCCCGGTTTTTTTATAAAAAAATTATGGAAAAAAATGACAGGAAAAAAGG‑AGAAA‑TTCTAAATAAATGCGTTAACTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                        <  2:333529/149‑1 (MQ=255)                                                 ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                     >  2:531336/1‑150 (MQ=255)                                                 ttATGCACCCCGGTTTTTTTATAAAAAAATTATGGAAAAAAATGACAGAAAAAAAGG‑AAAAA‑TTCTAAAAAAATGCGGTAACTTAAAGATTAGGTTGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                     <  2:968666/150‑1 (MQ=255)                                                  tATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                    >  2:215358/1‑150 (MQ=255)                                                   aggccgCCCGGTTTTTTTATAAAGAAATTATGGAAAAAAAGGACGGGGAAAAAGG‑AAAAA‑TTCTAAATAAATGCGTAAATTTAGAAATGAGGAGTGCGGAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                   <  1:924105/145‑1 (MQ=255)                                                      cacccgggtttttttATAAAAAATTTAGGGAGAAAAAGGAC‑GGGAAAAAGGAAAAAA‑TTCAAAATAAGTGCGGTAACTTAGAGATTAGGTTTGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACgg                                                                                                <  1:1119964/141‑1 (MQ=255)                                                       aGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                              >  2:32259/1‑151 (MQ=255)                                                         cccgttttttttAGAAAAAAATTAGGGAAAAAAAGGACAGGAAAAAAGG‑AAAAA‑TTCTAAAAAAATGCGTTAATTTAGAGATTAGGATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                               <  2:871976/142‑1 (MQ=255)                                                         cccgttttttttAAAAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGG‑AAAAA‑TTCTAAATAAATGCGTTAACTTAGAGATTAGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                   <  2:727718/138‑1 (MQ=255)                                                          cgggtttttttAAAAAAAATTTAGGAAAAAAAATGA‑AGGAAAAAAAG‑AAAAATTTCCCAATAAGTGCGGTAATTTAGAAAGAAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                              <  1:75570/143‑1 (MQ=255)                                                            ggtttttttAAAAAAAATGTATGGAAAAAAATGAGGGGAAAAAAGG‑AAAAT‑TTCCCAAAAACTGGGTTAATTTAGATATAAGGATTGCGGAAAAAAACATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                          <  2:691391/147‑1 (MQ=255)                                                                    taaaaaaaaaTTATGGAAAAAAATGACAGGAAAAAAGG‑AAAAA‑TTCTCAAAAAATGCGGTAACTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                  <  2:776480/143‑1 (MQ=255)                                                                        aagaaATTAGGGAAAAAAATGAAGGGAAAAAAGG‑AAAAA‑TGCTAAAAAAGTGCGGTATCTAAGAGTTTGGGATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                             <  2:107657/151‑1 (MQ=255)                                                                        aaaaaaTTATGGAAAAAAATGACAGGGAAAAAGG‑AAAAA‑TTCTAAAAAAATGCGGTAACTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                 <  2:752296/144‑1 (MQ=255)                                                                         agaaATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc                                                                            >  2:536448/1‑151 (MQ=255)                                                                           aaaTTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                          >  2:385865/1‑151 (MQ=255)                                                                           aaaTTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                          >  2:546606/1‑151 (MQ=255)                                                                            aaTTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                         >  2:1031767/1‑151 (MQ=255)                                                                                 tGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaa                                                                     >  1:1035481/1‑150 (MQ=255)                                                                                  ggAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaa                                                                    >  2:759241/1‑150 (MQ=255)                                                                                  gaaaaaaaaGGACGGGAAAAAAGG‑AAAAA‑TTCTAAATAAAGGCGGTAATTTAGAGATGGGGAGTGCGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaa                                                                    <  1:645452/146‑1 (MQ=255)                                                                                    agaAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                   >  2:269377/1‑149 (MQ=255)                                                                                    aaaaaaaTGACAGAAAAAAAGA‑AAAAA‑TTCAAAAAAAGGGCGTTAACTTAGAAATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg                                                                 <  2:1067138/149‑1 (MQ=255)                                                                                    aaaaaaaGGACAGGAAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAATTTAGAGATTAGGATTGGGAAAAACAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg                                                                 <  1:68412/149‑1 (MQ=255)                                                                                     aaaaaagggaagaaaaaaaGG‑AAAAA‑TTCTAAAAAAGTGCGGTATCTTAAAGATAAGTTTGGAGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa                                                                <  1:62940/137‑1 (MQ=255)                                                                                      aaaaaaGAAAGGAAAAAAGG‑AAAAT‑TTCTAAAAAATGGTGTTACCTTAGAGATTAGGATGGAGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                               <  1:1101085/151‑1 (MQ=255)                                                                                       aaaaTGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                              >  1:231717/1‑151 (MQ=255)                                                                                       aaaaTGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                              >  1:975309/1‑151 (MQ=255)                                                                                                ggaaaaaaGG‑AAAAA‑TTCTAAATAAATGGGTTAACTGAGAGATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                       <  1:1099933/146‑1 (MQ=255)                                                                                                  gAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                   >  1:288795/1‑151 (MQ=255)                                                                                                    aaaaGG‑AAAAA‑TTCTCAATAAATGTGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                 <  1:1057262/151‑1 (MQ=255)                                                                                                      aaGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                               >  2:590181/1‑151 (MQ=255)                                                                                                       aGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                              >  1:1142284/1‑151 (MQ=255)                                                                                                           aaaaa‑TTCTAAAAAAATGCGTTAACTTAGAAATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                           <  2:161139/149‑1 (MQ=255)                                                                                                                  tctcAATAAATGCGGTAACTTAGAGATTAGGATTGAGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                     >  1:683193/1‑151 (MQ=255)                                                                                                                      aaaaaaGGCGGTAACTAAGAGATTAGGTTGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                  <  1:514031/147‑1 (MQ=255)                                                                                                                      aaTAAATGCGTTAATTAAGAAATAGGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                      <  1:442888/146‑1 (MQ=255)                                                                                                                      aaTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                  >  1:175098/1‑150 (MQ=255)                                                                                                                      aaTAAAGCCGGTAACTTGGAGATTAGGATTGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                     <  1:463330/147‑1 (MQ=255)                                                                                                                       aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                >  2:1116879/1‑151 (MQ=255)                                                                                                                          aaggCGTTAACTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                             <  2:27806/148‑1 (MQ=255)                                                                                                                               gttAACTAAGAAATTAGGATGGGGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                           <  1:328420/146‑1 (MQ=255)                                                                                                                                   aCTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                    >  2:511794/1‑151 (MQ=255)                                                                                                                                      tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg                 >  2:276568/1‑151 (MQ=255)                                                                                                                                           aTTAGGATTGCGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag             <  1:32249/150‑1 (MQ=255)                                                                                                                                               gagtggaggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc               <  1:625077/137‑1 (MQ=255)                                                                                                                                                      ggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCg  <  2:891042/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                   ATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGG‑TTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGG‑AGAAA‑TTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCG  >  NC_000913/5010‑5305

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