Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	43.1%	85.1 / 19.3	58	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (6/27);  new base G (25/0);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 4.74e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925755‑926050                                                                                                                                                    |                                                                                                                                                        tCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCg                                                                                                                                                     <  2:495537/150‑1 (MQ=255)   aCTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTTG‑GGGGGGGGttt                                                                                                                                                      >  2:1229338/1‑146 (MQ=255)   aCTGCCGTACAGGCAGCTTAGAAATTCACAGGGAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGGGAAAAGACGGAAAAGAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGGGA‑GGGGGGGAGTCg                                                                                                                                                     >  1:494385/1‑148 (MQ=37)     tGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTAA‑GGGGGGGGTTTGAAccc                                                                                                                                                >  2:232473/1‑151 (MQ=255)        cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAAacccc                                                                                                                                              >  1:228072/1‑150 (MQ=255)        cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCg                                                                                                                                             <  1:676766/151‑1 (MQ=255)         gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAAaccccc                                                                                                                                             >  1:1012637/1‑149 (MQ=255)               gCAGCTTAGAAATTCACAGGTAACATACTCCACCCACCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGtt                                                                                                                                      <  1:302593/151‑1 (MQ=255)                cAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTg                                                                                                                                     <  1:942435/151‑1 (MQ=255)                          aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGG‑GGGGGGGGGTTGAAACCCCCCTAAGTTGCCGgat                                                                                                                               >  1:6838/1‑147 (MQ=25)                          aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATacac                                                                                                                           <  2:253080/151‑1 (MQ=255)                           ttCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTTGAACCCCCCCTACGTTGTCGGATacac                                                                                                                           >  2:235202/1‑150 (MQ=25)                              acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTCGAACCCCCCCTACGTTGCCGGATacaca‑c                                                                                                                        >  1:679713/1‑149 (MQ=255)                              acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑Ctt                                                                                                                      <  1:642224/151‑1 (MQ=255)                                 ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTcc                                                                                                                   <  2:132649/151‑1 (MQ=255)                                    aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAg                                                                                                                 <  1:445499/150‑1 (MQ=255)                                    aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGCGA‑GGGGGGGGGTTGAACCCCCCCTACCGTGCCGTATACACACCTTTccgg                                                                                                                 >  1:823593/1‑149 (MQ=21)                                        tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGGGA‑GGGGGGGGGTTGAACCCCCCACACGGTTCCGcatttaacccttcccgg                                                                                                                 >  1:376744/1‑133 (MQ=21)                                               cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTc                                                                                                     <  1:1042095/151‑1 (MQ=255)                                                  gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTCGAACCCCCCACACGTTGTCGTATACACA‑CCTTTCCGGCGGGCTcccttcgc                                                                                                  >  1:489369/1‑145 (MQ=18)                                                     caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCAc                                                                                               <  2:771552/151‑1 (MQ=255)                                                     caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCCCGCACTTTCGTAAGAACTCCTCTTTAAATTTGGCGGGTGG‑GGGGGGGATTCGAACCCCCCAACGGTTGCCGTGTAAACC‑CTTTCCAGGGGGGCTCCTTCTGccc                                                                                                >  2:894189/1‑150 (MQ=255)                                                       ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGGAACCCCC‑CCACGTTGTCGTATAAAAC‑CTTTTCAGGCGGGCTCCCTCTTCCAcc                                                                                              >  2:1191542/1‑149 (MQ=18)                                                       ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTc                                                                                             <  1:687561/151‑1 (MQ=255)                                                        caTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAAACCCCCAACCGTTTTCGTATACACA‑CCTTCCCGGCGGGCTCCTTCAGCCAccc                                                                                             >  2:1150611/1‑148 (MQ=18)                                                                  gCCGTACAGACTGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGAGA‑GGGGGGGGTTCGGAACCCCCAAAAGTTGCCGTATACACA‑CCTTCCAGGCGGGCTCCTTCTGcccccccgccacccccc                                                                                  >  1:1127829/1‑136 (MQ=17)                                                                  gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGG‑GGGGGGGGGTTGAAACCCCCCAACCTTGTCCTATACACA‑CCTTTCAAGCGGGGTCCCTCACCCCCTccggcaaccccc                                                                                  >  2:1148706/1‑140 (MQ=12)                                                                  gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCa                                                                                   <  1:141505/150‑1 (MQ=255)                                                                   ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTTGAAACCCCCCAAAGTTTCCGTATACAAA‑CCTTTCACGCGGGCCCCTTCAGCCACCCCGGCACCCCAcc                                                                                 >  1:944223/1‑151 (MQ=17)                                                                    cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCa                                                                                <  2:1184194/151‑1 (MQ=255)                                                                     gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCaa                                                                               <  1:1008544/151‑1 (MQ=255)                                                                         agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGGCACCTCACCAAATTg                                                                           <  2:563435/151‑1 (MQ=255)                                                                           acagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGtt                                                                         <  1:599279/151‑1 (MQ=255)                                                                            cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGAACCCCCCAAACGTTGCCGTATAAACA‑CTTTCCAGGCGTGCTCCTTCAGCCACCCCGGCACCCCCCCAAAATtttt                                                                        >  1:750427/1‑151 (MQ=14)                                                                              gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTg                                                                      <  2:134976/151‑1 (MQ=255)                                                                                       cGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAACCCCCCCTACGTTGCCGTTTACACA‑CCTTCCTGGCGGGGTCCTTCTGCCACTCGGGCACCCCCCCAAAATGTTTTTCttc                                                                  >  1:228380/1‑144 (MQ=11)                                                                                          aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTCGAACCCCCCATACGTTGCCGTATAAACA‑CTTTCCCGGCGGGCTCCTTCAGCCACCCCGGCACCCCACCCAAATGTTTTTCTTCCcaaac                                                            >  1:1026771/1‑147 (MQ=11)                                                                                          aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGAGTTGAAACCCCCCGACATTGCCGTATACAAA‑CCTTCCAGGCGTGCTCCTTCAGCCACTCGGGCACCCCCCCAAAATTTTTTTCTACCcaa                                                              >  1:849348/1‑147 (MQ=11)                                                                                           aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCa                                                         <  1:144462/151‑1 (MQ=255)                                                                                            aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAACCCCCCATACGTTTTCGTATAAACA‑CTTTTCAGGGGGGGTCCTTCCCCCCCTCGGGCCCCTCCCCCAAATTTTTTTCTACCaaa                                                              >  2:256141/1‑145 (MQ=9)                                                                                            aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGAAACCCCCCTACGTTGTCGGATACACA‑CTTTTCAGGCGGGCTCCTTCCCCCCCCCCGACAACCCCCCAAAATTTTTTTCTACCAAAACTc                                                          >  1:555018/1‑149 (MQ=9)                                                                                              aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGG‑GGGGGGGGATTGAACCCCCCCAACGTTTTCGTAAACACACCTTTCCCCGGGGGC‑CCCTCAGCCACCCCGGGAACTCCCCCAATTGTTTTTCTACCAAAACTc                                                          >  1:582714/1‑147 (MQ=9)                                                                                                    cGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTgg                                                  <  1:1111638/149‑1 (MQ=255)                                                                                                    cGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGc                                                 <  1:880414/150‑1 (MQ=255)                                                                                                         tttCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAACCCCCCATACGTTGCCGTATACAAA‑CCTTTCCGGCGGGGTCCCTCAGCCCCCCCGACACCCCCCCCAAATGTTTTTCTACCAAAACTCCTTGGGGGCAACggg                                           >  1:137071/1‑151 (MQ=9)                                                                                                         tttCGTACGAGCTCTTCTTTAAATATGGC‑GGGGG‑GGGGGGGGGTTGAAACCCCCCGACGTTGTCGTATAAACA‑CTTTTCCTGGGTGCTCCCTTCGCCCCCCCGACACCCCACCCAAATTTTTTTCTACCCAACCCCATGGGGGGGAACggg                                           >  2:1251541/1‑151 (MQ=9)                                                                                                                  aGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctacta                                  <  2:1149067/151‑1 (MQ=255)                                                                                                                   gCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAAACCCCGCTACGTTGCCGTATACACA‑CTTTTCAGGCGGGCTCCTTCAGCCACTCCGGCACCTCACCAAAATGTTTTTCTACCAAACCTCATGGGGGGCAACGGGGGGCTACtat                                 >  1:1002508/1‑151 (MQ=17)                                                                                                                      cttcttTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACtata                                <  2:944207/149‑1 (MQ=255)                                                                                                                         cttTAAATATGGC‑GGTGA‑GGGGGGGGGTTGAAACCCCCATACGTTGCCGGATACACA‑CTTTCCAGGCGGGCTCCTTCAGCCACCCGGACACCCCACCAAAATGTTTTTCTACCAAACCTCATGGGGGGCAACGGGGGGCTACTTTAggggg                           >  1:171340/1‑149 (MQ=17)                                                                                                                            tAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                         <  2:736060/150‑1 (MQ=255)                                                                                                                                   ggC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa                  >  2:217541/1‑150 (MQ=255)                                                                                                                                     c‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               <  1:1045423/151‑1 (MQ=255)                                                                                                                                       ggTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg              >  1:1117088/1‑151 (MQ=255)                                                                                                                                         tGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt            >  1:721015/1‑151 (MQ=255)                                                                                                                                                ggggATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa      <  2:157585/151‑1 (MQ=255)                                                                                                                                                ggggATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa      <  2:374872/151‑1 (MQ=255)                                                                                                                                                ggggATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa      <  2:376740/151‑1 (MQ=255)                                                                                                                                                    aTTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  >  1:722217/1‑151 (MQ=255)                                                                                                                                                    aTTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  1:466733/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        TCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925755‑926050

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