Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	42.6%	78.9 / 24.7	61	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (9/26);  new base G (26/0);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 7.21e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGA  >  NC_000913/925753‑926048                                                                                                                                                       |                                                                                                                                                        gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTc                                                                                                                                                     <  1:279286/151‑1 (MQ=255)          cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCG‑T‑GGGGGGGG‑TGGCAAccccc                                                                                                                                             >  2:431985/1‑150 (MQ=38)           gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑TTCGAAccccc                                                                                                                                             >  1:1035093/1‑149 (MQ=255)             aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAG‑G‑GGGGGGGG‑GTTGAAacccccca                                                                                                                                          >  2:941272/1‑147 (MQ=38)              cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATAc                                                                                                                                        <  1:633276/151‑1 (MQ=255)               aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑GTCGAAcccccccta                                                                                                                                         >  2:72519/1‑145 (MQ=255)                 gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGtt                                                                                                                                     <  1:630697/151‑1 (MQ=255)                   aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGc                                                                                                                                   <  1:976439/151‑1 (MQ=255)                        aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑GTTGGACCCCCCATAAGTTGCCg                                                                                                                                 >  2:95667/1‑148 (MQ=37)                          aaaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTCCGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATac                                                                                                                            <  1:579791/151‑1 (MQ=255)                              tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAA‑GGGGGGGG‑TCGGAACCCCCCA‑ACGTTGCCGTATacaca‑c                                                                                                                       >  1:533823/1‑151 (MQ=25)                              tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGG‑A‑GGGGGGGG‑TTTGAAACCCCCATACGTTGCCGtataaa                                                                                                                           >  2:213896/1‑146 (MQ=25)                                acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑GTTCAACCCCCCATAAGTTTCCCTCTaaaca                                                                                                                         >  2:966766/1‑148 (MQ=25)                                  aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑TTCGAACCCCCCATACGTTGCCGTATaaacaactt                                                                                                                     >  1:254074/1‑146 (MQ=255)                                   ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTcc                                                                                                                  <  2:310426/151‑1 (MQ=255)                                   ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAG‑A‑GGGGGGGG‑GTTGAACCCCCGCAACGTTGccccattaacaccttc                                                                                                                    >  1:758401/1‑136 (MQ=25)                                      aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGG‑G‑GGGGGGGG‑TTTGAAACCCCCATACGTTTCCGTTTACAAA‑Cttt                                                                                                                    >  2:922435/1‑147 (MQ=21)                                      aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAg                                                                                                                <  2:872597/150‑1 (MQ=255)                                        cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCg                                                                                                             <  2:261324/151‑1 (MQ=255)                                            cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGc                                                                                                          <  1:119042/150‑1 (MQ=255)                                               cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑GTCGAAACCCCGCCACCTTGCCGTTTAAACA‑CTTTTCAggggggggcc                                                                                                       >  1:495520/1‑142 (MQ=21)                                                 cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑TTTTGACCCCCCATACGTTTCCGCATACAAC‑ATTTTCAGGGGGGCTcc                                                                                                       >  1:488473/1‑148 (MQ=21)                                                   cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAg                                                                                                  <  1:109684/151‑1 (MQ=255)                                                         ccaTGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTc                                                                                            <  1:181381/151‑1 (MQ=255)                                                               tCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑TTCGAAACCCCCCTAAGTTGCCGCATAAACA‑CCTTCCCGGCGTGCTCCTTCAGCCACTccggaa                                                                                       >  2:228100/1‑145 (MQ=17)                                                                  cTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑TTCGGACCCCCCATACGTTTCCGTATACACA‑CTTTCCAGGCGGGCTCCTTCAGCCCCTCCGACA‑CCTc                                                                                  >  2:972566/1‑151 (MQ=17)                                                                    gCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCta                                                                                  <  2:349379/149‑2 (MQ=255)                                                                    gCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCAc                                                                                <  1:417314/151‑1 (MQ=255)                                                                     ccGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGG‑A‑GGGGGGGG‑TTTGAACCCCCCACACGTTTCCGTATAAACACCTTTCCCGTCGGGC‑CCCTCAGCCACCCCGGCA‑ccccccc                                                                               >  2:965916/1‑146 (MQ=12)                                                                      cGTACAGACAGATAAAATGCG‑‑AAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGG‑A‑GGGGGGGG‑GTTGAACCCCCGAAACAGTGCCGGATAAACA‑ATTTTCAGGGGTGCTCCTTCAGGCACCCGGGCA‑ACCCAcc                                                                               >  2:184086/1‑149 (MQ=14)                                                                          cagacagaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAAtt                                                                          <  1:938836/151‑1 (MQ=255)                                                                           agacagaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTg                                                                         <  1:691861/151‑1 (MQ=255)                                                                            gacagaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑GTTGAAACCCCCATAAGGTGCCGGATACACA‑CCTTCCCGGCGGGCTCCTTCAGCCCCCCCGACA‑CCCCCCCaaaattt                                                                        >  2:1008613/1‑147 (MQ=14)                                                                            gacagaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGt                                                                        <  2:649763/151‑1 (MQ=255)                                                                              cagaTAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTT‑G‑GGGGGGGA‑TTCCAACCCCCCTTACGATGCCGTATTCAAC‑CCTTCCAAGGGTGCGCCCTCATCCCCCCGGGGA‑CCTCACCCAAATGttt                                                                      >  2:117575/1‑151 (MQ=9)                                                                                   aaaaTGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGACCCCCCGATCCGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGTACC‑CCTCACCAAATTGTTTTGCt                                                                  <  1:58903/150‑1 (MQ=255)                                                                                   aaaaTGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAG‑AGGGGGGGGG‑TTGTAACCCCCCA‑ACGTTGCCCTATAAACA‑CCTTCCAGGCGGGCTCCTTCACCCCCCCCGGCACCCCCCCCAAATtttttt                                                                     >  1:430080/1‑148 (MQ=9)                                                                                      aTGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAG‑A‑TGGGGGGA‑GTTGCACCCCCCCAACGTTTGCGCATACACA‑CTTTCCCGGCCGGCCCCCTCCGCCCCTCCGACA‑CCCCCCCCAAATTTTTTGCTTcc                                                               >  2:578051/1‑150 (MQ=9)                                                                                          g‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑TTTGAACCCCCCCTACGTTTCCGGATAAAAA‑CCTTCCAGGCGGGGCCCTTCCGCCACCCCGGCACCCCCCCCAAATTTTTTTgttccaaa                                                             >  2:552958/1‑142 (MQ=2)                                                                                            aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑TTCGAACCCCCGATACGTTGCCGTATAAACA‑CTTTTCAGGCGGGCTCCTTCAGCCACCCGGACA‑CCTCACCAAATTTTTTTGCTACCCAAcc                                                          >  2:309430/1‑150 (MQ=14)                                                                                              aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑TTCGAACCCCCGCTACGTTTTCGTATATACA‑CTTTTCCGGCGTGGCCCTTTCGCCACCCGGGAA‑CCCCCCCCAAATGTTTTTCTATCaaaac                                                          >  1:707767/1‑146 (MQ=2)                                                                                                aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAG‑G‑GGGGGGGAGTTTGAA‑CCCCGCAAAGTTTCCCTATACACA‑CTTTCCCGGCGGGCTCATTCCGCCACACCGGCA‑CCCCCCCCAAATGTTTTTTTAccccacc                                                          >  1:485876/1‑141 (MQ=2)                                                                                                  aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGG‑TTCGGACCCCCGAAAACGTGCCGTAAACAAA‑CTTTCCAGGGGTGCCTCCTCCGCCCCTCCGACACCCCCACCAAATT‑TTTTGCTTACAAACCTCATggg                                                   >  1:146624/1‑151 (MQ=9)                                                                                                     cTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑AGGGGGGGGG‑TTCAAACCCCC‑AAACGTTGCCCGATAAACA‑CCTTCCCGGCGTGCTCCCTCTCCCCCTCGGGCA‑CCCCACCCAAATATTTTTCCTCCAAAcccccttggggg                                                >  2:51631/1‑141 (MQ=2)                                                                                                       cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTAATGGGTGGc                                               <  1:724570/150‑1 (MQ=255)                                                                                                             ttCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATCCGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACC‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGTCAACggg                                         <  1:164959/150‑1 (MQ=255)                                                                                                                   cGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctac                                  <  2:218682/151‑1 (MQ=255)                                                                                                                   cGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑GTTGAACCCCCCCTACGCTGGCGTATACAAA‑CCTTCCAGGCGGGCTCCTTCACCCACTCCGGCA‑CCCCACCCAAATGTTTTTCTACCAAAACTCATGGGGGGCAACGGGGGGctac                                  >  2:7214/1‑151 (MQ=17)                                                                                                                        tcttctTTAAATATGGCGGTG‑A‑GGGGGGGG‑TTCGAACCCCCCACACCTTTGCGTATAAACA‑CTTTTCAGGGGGGGTCCTTCACCCCCCCCGACA‑CCCCACCCAATTTTTTTTCTTCCCAACCCCCTGGGTGGGAACGGGGGGCCCCtata                              >  1:50788/1‑150 (MQ=9)                                                                                                                        tcttctTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAAACCCCCCTACGTTTCCGGATACACA‑CTTTCCAGGCGGGCTCCTTCCGCCACCCGGACA‑CCCCCCCAAAATGTTTTTCTACCAAACCCCATGGGGGGCAAAGGGGCGctac                                  >  2:750362/1‑146 (MQ=17)                                                                                                                          ttcttTAAATATGGCGGTG‑AGGGGGGGGG‑TTGAAACCCCCGAT‑TGGTGTCGTATTCACA‑CTTTTCATGCGTGCTCCCCCTCCCCCCCTGACA‑CCCCACCACAATCTTTTTTTACCAACCCCCCTTGGGGGGAACGGGGGGCTCCcatatg                            >  1:174479/1‑148 (MQ=2)                                                                                                                              ttAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGa                          <  2:495777/148‑1 (MQ=255)                                                                                                                                  atatGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGGTGCCGTATACACA‑CTTTCCAGGCGTGCCCCTTCAGCCACTCGGACA‑CCTCACCAAAATGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGa                          >  1:216545/1‑144 (MQ=255)                                                                                                                                        cGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc             <  1:778323/151‑1 (MQ=255)                                                                                                                                          gTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACgg           >  1:471363/1‑151 (MQ=255)                                                                                                                                          gTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACgg           >  2:974625/1‑151 (MQ=255)                                                                                                                                            g‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc         <  2:237503/151‑1 (MQ=255)                                                                                                                                              a‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa        <  2:292312/151‑1 (MQ=255)                                                                                                                                                gggggggA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaa       <  2:348370/151‑1 (MQ=255)                                                                                                                                                 ggggggA‑TTCGCACCCCCGCTACGTTGCCGTATACACA‑CTTTTCAGGCGTGCTCCTTCAGCCCCTCGGACT‑CCTCACCAAATTGTGTTGCTACCAAACCTCATGGGGGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAACCGGTCaag      <  1:488228/151‑1 (MQ=255)                                                                                                                                                      gA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAagaaga  >  2:333371/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                        GTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTG‑A‑GGGGGGGA‑TTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGA  >  NC_000913/925753‑926048

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