Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	53.6%	74.8 / 97.2	70	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (7/25);  new base G (37/0);  total (45/25)
Fisher's exact test for biased strand distribution p-value = 8.11e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-12

CGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925762‑926050                                                                                                                                           |                                                                                                                                                           cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GAAccccc                                                                                                                                                >  1:554932/1‑150 (MQ=255)  gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGa                                                                                                                                              <  2:1107223/151‑1 (MQ=255)   tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTC‑GAAcaccccc                                                                                                                                              >  2:993684/1‑148 (MQ=255)   tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTT‑GAAacccc                                                                                                                                                >  2:242007/1‑148 (MQ=255)        gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTT‑GAAccccccaaa                                                                                                                                            >  1:725557/1‑143 (MQ=255)         cAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGATT‑GAAACCCCGATACgct                                                                                                                                        >  1:588419/1‑148 (MQ=37)         cAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GAACCCCCCATACGtt                                                                                                                                        >  1:967332/1‑150 (MQ=255)          aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTg                                                                                                                                       <  1:1010489/150‑1 (MQ=255)           gCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGcc                                                                                                                                     <  1:538591/151‑1 (MQ=255)             ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTT‑GGAACCCCGAAACGGTGcc                                                                                                                                     >  1:517092/1‑149 (MQ=37)                          ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GAAACCCCCATACGTTTCCGTATaca                                                                                                                              >  2:463507/1‑143 (MQ=255)                           gTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GAACCCCCCATACGTTGCCGGATACAC‑A‑Cttt                                                                                                                      >  1:1063341/1‑148 (MQ=255)                           gTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTC‑GAACCCCCCATACGTTGCCGGATAAAC‑A‑Cttt                                                                                                                      >  1:1082644/1‑148 (MQ=255)                              aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGc                                                                                                                <  2:258868/151‑1 (MQ=255)                                aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCTGAACCCCC‑ATACGTTGCCggttaacaccttt                                                                                                                        >  2:526016/1‑131 (MQ=21)                                  aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GAAACCCCCCTACGTTGCCGTATAAAC‑A‑CCTTTCAGGCGTg                                                                                                             >  2:1029093/1‑150 (MQ=21)                                    tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGGTT‑TAAACCCCCAAACGGTGGCGGATACAA‑ACCTTTccgg                                                                                                                  >  2:297529/1‑142 (MQ=21)                                       aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GAAACCCCCGTACGTTTCCCTATACAC‑A‑CCTTTCAGGCGTGGTcc                                                                                                         >  2:819158/1‑149 (MQ=21)                                        cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTc                                                                                                      <  1:973608/151‑1 (MQ=255)                                          cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTT‑GAACCCCCCATACGTTTCCGGATACAC‑A‑CCTTTCAGGCGGGGTCCtt                                                                                                       >  2:997312/1‑148 (MQ=18)                                            cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGcc                                                                                                  <  1:466086/151‑1 (MQ=255)                                                      tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                                         <  2:632910/150‑1 (MQ=255)                                                           gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCa                                                                                    <  1:981280/150‑1 (MQ=255)                                                             cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAAGGGGGGGGTT‑GTACCCCCCCTACGTTGTCGGAAACAA‑C‑CATTCCCGGCGGGCTCCTTTAGTCAcacacacaacacccc                                                                                  >  1:885529/1‑136 (MQ=12)                                                                acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCaaa                                                                              <  1:124054/150‑1 (MQ=255)                                                                 cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTT‑GGAACCCCCATACGTAGCCGTATAACC‑A‑CTTTCCCTGCGGGCACCCTCACCCACACCGCCAACCCA‑CCaaaat                                                                            >  1:433757/1‑149 (MQ=9)                                                                     cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGGGGGGGGGGTTT‑GAACCCCCCAAACGTTGCCCGATACACAA‑CTTTCC‑GGCGGGCTCCTTCAGCCCCCCCGACACCTCA‑CCCAAATtttt                                                                        >  2:1080771/1‑151 (MQ=14)                                                                     cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTT‑GAACCCCCCCTACGTTGCCGTATACAC‑A‑CTTTTCCGGCGTGCTCCTTCAGCCACTCCGGCACCTCA‑CCCAAATtttt                                                                        >  2:958157/1‑151 (MQ=14)                                                                      agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGtttt                                                                       <  2:1003108/151‑1 (MQ=255)                                                                      agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGAGTC‑TAACCCCCCCTACGGTGCCGGATAAAC‑A‑CTTTCCCGGCGTGCTCCTTCTGCCACCCCGGCAACCCC‑CCAAAATGtttt                                                                       >  1:37393/1‑151 (MQ=14)                                                                           aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTT‑GAACCCCCGCTACGTTGCCGTATACAC‑A‑ATTTCCAGGCGTGCTCCCTCAGCCCCCCGGACACCCCA‑CCAAAATGTTTTTCt                                                                    >  1:438931/1‑149 (MQ=14)                                                                             aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTT‑GAAACCCCCATACGTTGCCGTATAAAC‑A‑CTTTTCCGGCGTGCTCCTTCAGCCACCCGGACACCCCA‑CCAAAATTTTTTTTTAcc                                                                 >  1:487053/1‑150 (MQ=11)                                                                             aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCa                                                                <  1:56749/151‑1 (MQ=255)                                                                               gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTT‑GAACCCCCCATACGTTTCCGGATACAC‑A‑CTTTTCAGGCGTGGTCCTTCAGCCACTCCGGCACCCCA‑CCAAAATGttttttttccca                                                               >  1:29146/1‑142 (MQ=11)                                                                               gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCaaa                                                              <  1:1033999/151‑1 (MQ=255)                                                                                   aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GAAACCCCCATACGTTGCCGTATAAAC‑A‑CTTTCCAGGCGGGCTCCTTCCGCCACTCCGACACCCCC‑CCCAAATGTTTTTCTTACaaaac                                                            >  1:1134019/1‑147 (MQ=11)                                                                                       aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGAGGGGGGGGTTC‑GAAACCCCCATACGCTTCCGTATACAC‑A‑CTTTCCCGGCGGGCTCCTTCAACCACTCCGACACCCCA‑CCAAAATTTTTTTTTACCCAACCCCATgg                                                      >  1:278188/1‑151 (MQ=14)                                                                                        aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GAACCCCCCATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACCCGGACACCCCA‑CCAAATTGTTTTTCTACCAAACCTCATggg                                                     >  1:397699/1‑151 (MQ=255)                                                                                             cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGTGGGGGGGGTTT‑CAACCCCCCCTAAGTTGCCGTATAAAC‑A‑CTTTTCAGGCGTGGTTCTTCCCCCACTCCGGCACCCCA‑CCCAAATTTTTTTCTATCAAAACTCATtggtggg                                                 >  2:218434/1‑149 (MQ=9)                                                                                             cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAATCCTCATGGGTGGCa                                                <  2:1105204/151‑1 (MQ=255)                                                                                              gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCaa                                               <  2:526700/151‑1 (MQ=255)                                                                                               cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTAGGGGGGGCTTT‑GAACCCCCCGTACCTTGCCGTATACCC‑C‑ATTTCCCCGGGGGGTCCTTCACCCACACGG‑CACCCCACCCCCAATGCTTTGCTTCCACACCCCATGGGGGGCaa                                               >  1:822697/1‑150 (MQ=2)                                                                                               cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCTAAGGGGGGGTTC‑TAACCCCCCCTTAGTTGCCGTATAAAA‑A‑CTTTCCAGGCGTGCTCCCTCAGCCCCCCGGACACCCCA‑CCCAAATTTTTTTCTACCAAAACCCCTGGGGGGcga                                               >  2:647035/1‑148 (MQ=2)                                                                                                aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑CGATCCCCCACACGGTGTCGTACAAAC‑C‑CCTTCCAAGCGGGGTCCCTCAGCCACTCCCGCCCCTCA‑ACAAATTTTTTTTCTTCCAAACCCCCTTGGGGGCaa                                               >  1:378474/1‑149 (MQ=2)                                                                                                aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACg                                             <  2:279757/151‑1 (MQ=255)                                                                                                  tttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTT‑GAACCCCCCCTAAGTTGCCCGATAAAC‑A‑CCTTTCCGGCGGGCCCCTTCCGCCCCTCGGGCAACCCC‑CCCAAATGTTTTTCTTCCAAACCTCATGGGGGGGAACggg                                           >  1:173255/1‑151 (MQ=9)                                                                                                  tttCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTT‑TACGCCCCCATCCGTTGGCTTAGACTC‑A‑CATTCCCGGCTTGGATCCTCAAAAACGCCGACACCGCA‑ACCAAATGTTTTTCTACCCACCCTcctaggg                                                    >  2:698822/1‑136 (MQ=2)                                                                                                   ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTT‑GGAACCCCCATACGTTTCCGGATACAC‑A‑CTTTTCAGGCGGGGTCCTTCAGCCCCCCCGGCACCCCACCCAAATTTTTTTG‑TACCAAACCCCATGGGGGGCAACcggg                                          >  1:5126/1‑151 (MQ=9)                                                                                                   ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACggg                                           <  2:198282/150‑1 (MQ=255)                                                                                                      gTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GCACCCCCCATAAGGTGTCGGATACAC‑A‑CTTTCCCGGCGTGCGCCCTCAGCCACTCGGACACCCCA‑CCCAAATGTTTTTTTTCCCAAACTCATGGGGGGGAAAGGGgcg                                        >  1:437747/1‑150 (MQ=2)                                                                                                      gTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GAACCCCCCATACGTTTCCCTATACAC‑A‑CTTTCCAGGCGGGCCCCTTCACCCCCCCGGACCCCCCC‑CCAAATTTTTTTGTTACCCAACCCCATGGGGGGGAACgggg                                          >  1:395683/1‑148 (MQ=14)                                                                                                       tACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTC‑GAACCCCCCATAAGTTTTCGTATACAC‑A‑CCTTCCAGGCGGGCTCCTTCAGCCCCACGGACACCCCC‑CCAAAATGTTTTGTTATCAAACCCCCTGGGTGGGAACGGGgcg                                        >  2:603794/1‑149 (MQ=14)                                                                                                       tACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTC‑GAAACCCCCATACGTTGCCGTATAAAC‑A‑CTTTTCAGGCGTGCTCCCTCAGCCACTCCGACACCCCA‑CCAAAATTTTTTTCTACCAAAACTCATGGGGGGCAAAgggggg                                        >  1:625652/1‑147 (MQ=14)                                                                                                             ctctTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GACCCCCCGATACGTTGCCATATACTC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCTCTCGGACACCTCG‑CCAATGTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGGGCTTCtat                                 <  2:166148/150‑1 (MQ=17)                                                                                                               cttcttTAAATATGGCGGTGAGGGGGGGGTTT‑GAACCCCCGATACGTTGCCGGATACAC‑A‑CTTTTCAGGCGGGCTCCTTCAGCCACCCGGACACCTCA‑CCCAAATGTTTTTCTACCAAACCTCATGGGGGGCAACGGGGGGCTACTATAgg                              >  1:591656/1‑151 (MQ=18)                                                                                                                ttcttTAAATATGGCGGTGAGGGGGGGGGTC‑GAACCCCCCATAAGTTTACGTATAAAC‑A‑ATTTTCCGGGGGGGTCCTTTAGCCACCCGGCCACCCCA‑CCAAATTTTTTTTTTACCCAACCTCATGGGGGGCGACAGGGGGCCACTATAggg                             >  2:393617/1‑151 (MQ=9)                                                                                                                    ttAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACTCCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAg                           <  1:952138/149‑1 (MQ=255)                                                                                                                      aaaTATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTg                        <  1:94074/150‑1 (MQ=255)                                                                                                                       aaTATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGa                      <  1:972066/151‑1 (MQ=255)                                                                                                                            ggCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa                  <  1:242051/150‑1 (MQ=255)                                                                                                                              cGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               >  2:502722/1‑151 (MQ=255)                                                                                                                              cGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               >  1:1135580/1‑151 (MQ=255)                                                                                                                               ggTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg              >  1:465517/1‑151 (MQ=255)                                                                                                                               ggTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg              <  2:693207/151‑1 (MQ=255)                                                                                                                                  gAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           <  2:1055904/151‑1 (MQ=255)                                                                                                                                  gAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           >  2:1003974/1‑151 (MQ=255)                                                                                                                                   aGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGCTGGAGTAAAACGGTc           >  1:1028517/1‑150 (MQ=255)                                                                                                                                     ggggggATTC‑GAACCCCCGATACGTTGCCGTATACAC‑G‑CTTTCCAGACGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTAGCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaag        <  2:762858/151‑1 (MQ=255)                                                                                                                                      gggggATTC‑GAACCCCCGATTCGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga       <  1:196039/151‑1 (MQ=255)                                                                                                                                         ggATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga       >  1:591746/1‑148 (MQ=255)                                                                                                                                           aTTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  1:420797/151‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                           CGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTC‑GAACCCCCGATACGTTGCCGTATACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925762‑926050

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