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

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAA  >  NC_000913/925751‑926049                                                                                                                                                        |                                                                                                                                                          gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGa                                                                                                                                                           <  2:202855/150‑1 (MQ=255)   gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GT‑AGGGGGGAATc                                                                                                                                                        >  2:235773/1‑151 (MQ=255)     tCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCg                                                                                                                                                       <  1:46846/150‑1 (MQ=255)      cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGaa                                                                                                                                                     <  1:1186488/151‑1 (MQ=255)          gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAAcccc                                                                                                                                                 <  1:73405/151‑1 (MQ=255)           ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GA‑GGGGGGGGGTTGAAacccc                                                                                                                                                >  1:263965/1‑151 (MQ=255)             gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑GGGGGGGGGTCGAAacccc                                                                                                                                                >  1:494130/1‑149 (MQ=255)             gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTAAAAacccc                                                                                                                                                >  2:267191/1‑149 (MQ=255)             gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑AGGGGGGGGTTGAAacccc                                                                                                                                                >  2:949870/1‑149 (MQ=255)              tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTTGAAacccc                                                                                                                                                >  2:304245/1‑148 (MQ=255)              tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GA‑GGAGGGGAATCggaaacccgc                                                                                                                                              >  1:352085/1‑141 (MQ=38)               aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATa                                                                                                                                            <  1:799863/151‑1 (MQ=255)                cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGG‑GGGGGGGGTTGGAAccccc                                                                                                                                                >  2:1056523/1‑147 (MQ=255)                cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCCAAccccccctc                                                                                                                                            >  2:698860/1‑146 (MQ=255)                  ggCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGt                                                                                                                                         <  2:592065/151‑1 (MQ=255)                          aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCGAAACCCCCATAAGGTGcc                                                                                                                                     >  1:1162298/1‑147 (MQ=37)                          aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGta                                                                                                                                  <  1:197114/150‑1 (MQ=255)                                tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATacaca                                                                                                                            <  2:842925/150‑1 (MQ=255)                                       tAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGGACCCCCCCCACGTTTCCCTATACACA‑A‑TTTcc                                                                                                                    >  1:636901/1‑149 (MQ=25)                                          cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCATCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCGGACCCCCCATACGGTGCCGGATACACA‑C‑ttttccgg                                                                                                                 >  2:888554/1‑144 (MQ=21)                                                 cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCt                                                                                                        <  2:1038444/151‑1 (MQ=255)                                                  aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GG‑AGGGGGGGTTTGGACACCCCATACGTTTCCGTATACACA‑C‑ctt                                                                                                                      >  1:475687/1‑136 (MQ=21)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTGAACCCCCGGTACGTTGCCGGATACACA‑C‑TTTCCCGGCGTGCTcc                                                                                                         >  2:491007/1‑148 (MQ=255)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCGAACCCCCCATACGTTTCCGCATACACA‑C‑TTTCCCGGCGGGCTcc                                                                                                         >  1:275907/1‑148 (MQ=21)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCa                                                                                                     <  1:769504/151‑1 (MQ=255)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCa                                                                                                     <  2:950194/151‑1 (MQ=255)                                                      gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCGGAACCCCCCTAAGTTGCCGTATAAACACC‑TTTCCA‑GCGGGCTCCTTCAGc                                                                                                   >  2:92459/1‑151 (MQ=18)                                                      gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCGAACCCCCCCTTAGTTGCCGTGTACACA‑A‑ATTTCCGGCGGGCGCCTTCAGc                                                                                                   >  1:910200/1‑151 (MQ=18)                                                        ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGGACCCCCCACACGTTTTCGTATAAACA‑C‑CTTTCCGGCGTGCTCCTTCCGcc                                                                                                  >  2:192842/1‑150 (MQ=18)                                                        ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑a                                                                                                <  2:476780/151‑1 (MQ=255)                                                         caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑Ac                                                                                               <  2:83168/151‑1 (MQ=255)                                                          accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACt                                                                                              <  2:1199733/151‑1 (MQ=255)                                                              tGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCGAAACCCCCCCAACTTGTCCTTTAAACA‑C‑CTTTCCGGGGGGCCCCTTCAGCC‑ACTc                                                                                             >  1:478666/1‑148 (MQ=18)                                                                ttCCCTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGac                                                                                         <  2:73659/150‑1 (MQ=255)                                                                       ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑GGGGGGGGGTCGAAACCCCCCTACGGTGCCGCATAAACA‑A‑TTTTCCGGGGTGGTCCTTCTGCC‑CCCCGGACAccccccc                                                                                 >  1:980446/1‑146 (MQ=17)                                                                        cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGAAACCCCCATACGTTTTCGTATAAAAA‑C‑TTTCCCCGCGTGCTCCCTCCGCC‑ACCCCGACAcccccc                                                                                  >  2:905220/1‑145 (MQ=17)                                                                             agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAAtt                                                                            <  2:476629/150‑1 (MQ=255)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑TA‑GGGGGGGGTTCGAAACCCCCATACGTTGTCCGATACACA‑C‑TTTTCAGGCGGGGTCCCTCCACC‑ACCCCGACAACCCACCAAAATGtttt                                                                       >  2:641839/1‑151 (MQ=14)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGtttt                                                                       <  1:92497/151‑1 (MQ=255)                                                                                  gaTAAAATGCGAAAAAAAAGCTCACACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGATCGAACCCCCCAAAAGTTGCCGTATACAAA‑C‑CTTTCAGGCGTGGTCCTTCAGCC‑ACCCCGACAACTCCCCCAAATGtttt                                                                       >  2:764639/1‑150 (MQ=14)                                                                                     aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGAACCCCCCATACGTTGCCGTATAAACA‑C‑CTTCCAGGCGTGCTCCTTCCGCC‑ACCCGGACCCCCCCCCCAAAtttttt                                                                       >  1:752911/1‑147 (MQ=14)                                                                                      aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGAACCCCCCATACGGTGCCCGGTAAACC‑CATTTCCGGCGGGGCTCCTTC‑CCC‑CCTCGGGCACCCCACCAAAAtttttttttt                                                                   >  2:1164350/1‑146 (MQ=9)                                                                                          gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACC‑C‑TTTCCAGGCGTGCTCCTCCAGCC‑GCTCGGACACCTCACCAAATTGTTTTGCTACCaaa                                                              <  1:522272/151‑1 (MQ=255)                                                                                          gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCaaa                                                              <  1:269552/151‑1 (MQ=255)                                                                                            gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGGACCCCCCCTACGTTTCCGCATACACA‑C‑TTTCCAGGCGGGCTCCCTCAGCC‑ACTCGGACACCTCCCCCAAATGTTTTTCTAACaaa                                                              >  2:1106/1‑149 (MQ=14)                                                                                              aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑TT‑GGGGGGGGTTTGAAACCCCCATTCGTTTTCCTATAAACA‑C‑ATTCCAGGGGGGCGCCCTCAACC‑CCACCGGCACCTCACACAAAAGTTTTgtttacaaaac                                                            >  1:478327/1‑139 (MQ=2)                                                                                                aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCAt                                                        <  1:282147/151‑1 (MQ=255)                                                                                                  aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTGGACCCCCCATACGTTTCCGTATACACA‑A‑TTTCCCGGCGGGGTCCCTCAGCC‑CCCCCGGCACCCCCCCCAAATGTTTTTCTTCCCAACCCCATgg                                                      >  1:435038/1‑151 (MQ=9)                                                                                                     gCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑GGGGGGGGTTCGGACCCCCCCAACGGTGCCGTATAAACA‑C‑CTTCCCGGCGGGCGCCCTCAGCC‑ACTCCGGCACCCCCCCCAAATTTTTTTTTTCCAAAACCCATggggg                                                   >  1:491618/1‑149 (MQ=2)                                                                                                        cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCa                                               <  1:659285/151‑1 (MQ=255)                                                                                                         gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCaa                                              <  1:873072/151‑1 (MQ=255)                                                                                                          cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAAc                                             <  1:1103094/151‑1 (MQ=255)                                                                                                              ttCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACggg                                          <  1:670646/150‑1 (MQ=255)                                                                                                              ttCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACg                                            <  1:569176/148‑1 (MQ=255)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGgcg                                       <  2:936648/151‑1 (MQ=255)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGgc                                        <  2:980496/150‑1 (MQ=255)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGGTCGAAACCCCCAAACGTTGCCGTATACACA‑C‑CTTTCCGGCGTGCTCCTTTAACC‑ACTCGGACACCCCCCCAAAATTTTTTTTTTCCAAACCTCCTGGGT‑GGCAAAGGGgcg                                       >  2:633643/1‑151 (MQ=9)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTGGGACCCCCCGCACCTTGCCGGATAAAAA‑C‑TTTTCACGCGGGCCCCTTCACCC‑CCCCCGGCCCCCCACCAAAATTTTTTTTTTCCAAAACCCCTgggg‑gggcaccggg                                         >  1:1009799/1‑140 (MQ=2)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGgcgc                                      <  1:895393/151‑1 (MQ=255)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGCGGA‑GG‑GGGGGGGGTTTGAAACCCCCATCAGTTTCCGTATTAAAA‑C‑TTTTCCGGCGGGGTCCTTCTGCC‑ACCCGGGCACCCCCCCAAAATTTTTTTTTACCCAACCCCATGGGG‑GGCAACGGGgcg                                       >  1:1186391/1‑150 (MQ=9)                                                                                                                     gAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTGAAACCCCCAAACGTTGCCGTATACACA‑C‑CTTCCCGGCGTGGTCCTTCAGCC‑ACCCCGGAACCCCACCCAAATTTTTTTCTACCAAAACCCCTTGGG‑GGGAACGGGgggcg                                     >  2:416053/1‑147 (MQ=2)                                                                                                                        ctctTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTTAACCCCCCATACGTTGCCGTATACACA‑C‑TTTTCAGGCGTGCTCCTTCCACC‑ACCCCGACACCCCCCCCAATTTTTTTGCTTCCCAACCCCATTGGG‑GGCAACGGGGGGccact                                  >  2:893903/1‑148 (MQ=9)                                                                                                                          cttcttTAAATATGGCGGT‑GA‑GGGGGGGGTTTGAAACCCCCATACGTTTCCGTATACACA‑A‑ATTCCACGGGGGCTCCTTTAGCC‑ACCCGGGCACCTCACCAAAATGTTTTTTTAACAAAACCCACTGGT‑GGGAACgggggggtcat                                  >  2:818793/1‑139 (MQ=2)                                                                                                                            tcttTAAATATGGCGGT‑GA‑GGGGGGGGACTGCAACCCCGATAAGCTTCCGAATTCACA‑A‑ATTTCCGGTGCGCTCCTCCTGCCAACCGGGACACC‑CAACAAATTGTTTTTTTAATAAACCCCTTTGGTGGGGAACGGGG‑GCTACtatatg                             >  1:764882/1‑147 (MQ=255)                                                                                                                                   atatGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGGCGCTACTATAGGGAGTTGGAg                    <  2:180787/151‑1 (MQ=255)                                                                                                                                   atatGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GCCAACGGGGCGCTACTATAGGGAGTTGGAg                    <  2:275922/151‑1 (MQ=255)                                                                                                                                      tGGCGGT‑GA‑GGGGGGGATTCAAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGGCGCTACTATAGGGAGATGGAGTa                  >  2:1114349/1‑150 (MQ=255)                                                                                                                                       ggCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaa                <  2:609903/151‑1 (MQ=255)                                                                                                                                              gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt           <  1:177233/150‑1 (MQ=255)                                                                                                                                               a‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa         <  1:855034/151‑1 (MQ=255)                                                                                                                                                  ggggggATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaag       <  2:910245/151‑1 (MQ=255)                                                                                                                                                     gggATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag    <  2:614765/151‑1 (MQ=255)                                                                                                                                                       gATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa  <  1:1131086/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                          GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCC‑ACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGT‑GGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAA  >  NC_000913/925751‑926049

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