Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,581	0	A	G	42.4%	97.0 / 18.7	61	noncoding (72/88 nt)	serX	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (4/30);  new base G (25/0);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 7.84e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAAAAA  >  NC_000913/1097445‑1097731                                                                                                                                            |                                                                                                                                                           cACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGAt                                                                                                                                              <  1:56875/150‑1 (MQ=255)  aCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTTGAAacccc                                                                                                                                                 >  1:265532/1‑147 (MQ=9)   cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGGACCCCCGATa                                                                                                                                             >  1:45686/1‑150 (MQ=255)    ggTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGGACCCCCCATAc                                                                                                                                            >  2:615573/1‑150 (MQ=255)       cTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGGTTGAAACCCCCGA‑ACGtt                                                                                                                                         >  1:399727/1‑150 (MQ=1)         gTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGGACCCCCCATACGtt                                                                                                                                         >  2:546296/1‑148 (MQ=255)         gTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTg                                                                                                                                        <  2:668073/149‑1 (MQ=255)              aCTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGAAACCCCCATAAGTTGCCGgat                                                                                                                                  >  2:808662/1‑150 (MQ=11)                       cGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑ttt                                                                                                                      <  2:518431/150‑1 (MQ=255)                        gCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTcc                                                                                                                    <  2:360822/151‑1 (MQ=255)                         ccAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCa                                                                                                                   <  1:219300/151‑1 (MQ=255)                            gTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGc                                                                                                                <  2:612037/151‑1 (MQ=255)                              tGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCg                                                                                                               <  2:195510/150‑1 (MQ=255)                               gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TCTCCAGGCGTg                                                                                                             <  2:87010/151‑1 (MQ=255)                                 cGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTCGAACCCCCCATACGTTTCCCCATTCAC‑A‑C‑CTTCCAGGCGTg                                                                                                             >  2:507959/1‑149 (MQ=14)                                    aCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTCGAACCCCCCCTACGTTTGCGTATAAAC‑A‑C‑TTTTCCAGCGTGCTcc                                                                                                         >  2:580960/1‑150 (MQ=14)                                       ccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTCGGAACCCCCATACGTTTCCGCATACAC‑A‑C‑CTTTCACGCGCGCTCCtt                                                                                                       >  2:681674/1‑149 (MQ=14)                                            cGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTCGAACCCCCCCTACGTTGCCGTATACAC‑C‑C‑CTTCCCGGGGGGCTcccttctcccc                                                                                                >  2:174318/1‑142 (MQ=17)                                             ggggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCa                                                                                                 <  1:649060/149‑1 (MQ=255)                                              gggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTTGAAACCCCCCTACGTTTCCGGATACAC‑ACC‑TTTCCA‑GCGGGCTCCTTCCGcccc                                                                                                >  2:69995/1‑147 (MQ=14)                                               ggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTCGAACCCCCCCTACGTTGCCGTATACAC‑A‑C‑CTTCCCGGCGGGCTCCTTCAGcccccc                                                                                              >  2:530538/1‑146 (MQ=18)                                               ggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCg                                                                                             <  2:981632/151‑1 (MQ=255)                                                gCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGAAACCCCCCTAAGTTGCCGTATACAC‑A‑C‑TTTCCCCGGGGGCTCCTTTCGCCAccc                                                                                              >  2:885128/1‑147 (MQ=14)                                                  ttCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTTGAAACCCCCCTACGTTGCCGTATACAC‑A‑C‑TTTTCCGGCGTGCTCCTTCACCCCCTCgg                                                                                            >  2:517959/1‑149 (MQ=17)                                                     tcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGG‑GG‑GGGGGGGTTTCGAACCCCCCA‑ACAGTTTCGTATACAC‑A‑C‑ATTTCCGGGGCGCTCCCTTCGCCACCCgggga                                                                                         >  1:581322/1‑146 (MQ=2)                                                     tcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGAAACCCCCCTACCTTTCCGGATACAC‑A‑C‑TTTTCACGCGGGCTCCTTCAGCCACCCCGACAcc                                                                                       >  2:255276/1‑151 (MQ=17)                                                       aTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GG‑GGGGGGGGTTCGAACCCCCCACACGTTTTCGTATACAC‑A‑C‑TTTCCACGCGTGCTCCCTCAGCCCCCCGGGCAcccc                                                                                     >  2:255524/1‑149 (MQ=17)                                                          cccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTTGGACCCCCCCCACGTTTCCGGATAAAC‑C‑C‑CTTTCCGGCGTGGTCCCTTAGCCACCCCGACACCTcccc                                                                                  >  1:929499/1‑151 (MQ=14)                                                            cccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCaa                                                                                <  2:585787/151‑1 (MQ=255)                                                                    gtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTTGAACCCCCCCTACGTTTCCGCATACAA‑A‑C‑TTTTCCGGGGTGGGCCTTTTCTCACCCCGGCCCCCCACCAAATTTTCGt                                                                        >  1:1043840/1‑151 (MQ=9)                                                                      gCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTc                                                                      <  1:234277/151‑1 (MQ=255)                                                                      gCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTc                                                                      <  1:449467/151‑1 (MQ=255)                                                                       cAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑‑GGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCt                                                                    <  2:355906/151‑1 (MQ=255)                                                                           ataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGt                                                                  <  2:226594/150‑1 (MQ=255)                                                                             aCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GG‑GGGGGGGGGTCGAACCCCCCCTAACTTTCCCGATTAAC‑A‑CTTTTCCAGG‑GTGCTCCTTCTACCACCCTGCCCACCCACCAAACTATCGTTtctttctt                                                               >  2:304835/1‑151 (MQ=2)                                                                                   aaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGAACCCCCCATAAGTTGGCGTATAAAC‑A‑C‑CTTTCCGGCGGGCTCCTTCAGCCACTCGGACACCCCCCCCAAATTTCGTTCCTGTCtt                                                               >  2:420192/1‑145 (MQ=14)                                                                                              aCTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTa                                              <  1:490947/150‑1 (MQ=255)                                                                                                tttCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTTGAACCCCCCCTACGTTTCCGGATAAAC‑A‑C‑TTTTCAGGCGTGGTCCCTCCGCCACCCGGACACCCCCCCCAAATTTTGTTCCTGTTTTTCT‑GGGACGGGGGCCAAtt                                           >  2:191410/1‑151 (MQ=2)                                                                                                 ttCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTCGAAACCCCCATACGTTTTCGTATACAC‑A‑C‑TTTTCCGGGGGGCTCCCTCTCCCACCCCGACACCCCACCCAATTTTCGTTCCTGTCTTGTT‑GGGAAGGGCGCTAAttt                                          >  1:993132/1‑151 (MQ=9)                                                                                                  tCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAAttt                                          <  1:293467/150‑1 (MQ=255)                                                                                                    gTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAgg                                       <  1:622767/151‑1 (MQ=255)                                                                                                       cGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGTTTCAAACCCCCCA‑ACGTTGTCCCATACAC‑A‑C‑CTTTTTGGGGGGCTTCTTTACCCACCCGGGCAACCCACCAAAATTTTGTTTCTGTCTTTTTGGGAACggggccactttttgg                                       >  1:1012634/1‑137 (MQ=2)                                                                                                         aGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTTGAACCCCCCCTACGCTTCCGTATACAC‑A‑C‑CTTTCCGGCGTGCTCCCTCAGCCCCCCGGACACCCCACCAAAATGTTGTTCCTGTCTTGCT‑GGGACGGGCGCCAATTTTGGGaaa                                   >  2:777880/1‑150 (MQ=17)                                                                                                         aGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAAt                                  <  1:934376/151‑1 (MQ=255)                                                                                                          gCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAAt                                  <  2:789251/150‑1 (MQ=255)                                                                                                           ctctTCTTTAAATATGGC‑GGTGGG‑GGGGGGGGTTTGAACCCCC‑CTACGGTTCCGTTTAAACAA‑C‑TTTCC‑GGGGTGGTCCCTCAGCCCCCCGGCCACCCCACCCAAATTTCGTTTCTTTCTTTTT‑GGAACGGGCGCTAATTTTGGGgaa                                   >  1:664398/1‑148 (MQ=9)                                                                                                               tcttTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCAt                               <  1:638269/148‑1 (MQ=255)                                                                                                                  ttAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGAcc                           <  1:78715/149‑1 (MQ=255)                                                                                                                    aaaTATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACCCCTCTCCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAg                       <  1:927589/151‑1 (MQ=255)                                                                                                                     aaTATGGC‑GGT‑GA‑GGGGGGGGTTCGAACCCCCGCTAAGTTTCCGTATACAC‑A‑C‑TTTTCCGGCGGGCTCCTTCAGCCACTCGGGCACCTCACCAAATTTTCGTTCCTGTCTTTCT‑GGAACGGGCGCTAATTTTGGGAAATCATTACCTGAgg                      >  1:285163/1‑151 (MQ=18)                                                                                                                           gC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTcaaca                <  2:330856/151‑1 (MQ=255)                                                                                                                            c‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTcaacaa               <  1:419489/151‑1 (MQ=255)                                                                                                                             tggT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTcaaca                >  1:520874/2‑150 (MQ=255)                                                                                                                              ggT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACaaa              <  2:641321/151‑1 (MQ=255)                                                                                                                              ggT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACaaa              <  1:181435/151‑1 (MQ=255)                                                                                                                               gT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAAc             <  1:338523/151‑1 (MQ=255)                                                                                                                                  gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACtt           >  1:564457/1‑151 (MQ=255)                                                                                                                                   a‑GGGGGTGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGCCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTATTTTAGGGAAATCATGACCTGAGGTCAACAAACttt          >  1:511443/1‑151 (MQ=255)                                                                                                                                      ggggggATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACttttt        <  2:1043618/151‑1 (MQ=255)                                                                                                                                      ggggggATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACttttt        >  1:191397/1‑151 (MQ=255)                                                                                                                                       gggggATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTg       <  2:668481/151‑1 (MQ=255)                                                                                                                                            aTTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGaaaaa  <  1:387604/151‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                           CACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACAC‑A‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCT‑GGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAAAAA  >  NC_000913/1097445‑1097731

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