Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913925,9000AG39.0% 105.6 / 11.0 59noncoding (72/88 nt)serWtRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (5/31);  new base G (23/0);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 6.82e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925755‑926050
                                                                                                                                                   |                                                                                                                                                             
tCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTc                                                                                                                                                            <  1:212510/149‑1 (MQ=255)
      ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAAccccc                                                                                                                                                    <  2:700126/151‑1 (MQ=255)
          aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGG‑GGGGGGGATTTGAAacccc                                                                                                                                                    >  1:721074/1‑147 (MQ=255)
                   ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑TGGGGGGGATCAAACCCCCGACACGtttc                                                                                                                                          >  2:102054/1‑146 (MQ=37)
                   ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTTTAACCCCCCCTACGTTggc                                                                                                                                         >  2:271099/1‑147 (MQ=37)
                    tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGG‑GGGGGGGGTTTGAAtcccccgcacgg                                                                                                                                             >  2:456964/1‑137 (MQ=37)
                    tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCg                                                                                                                                        <  2:8723/149‑1 (MQ=255)
                     aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGCGA‑GGGGGGGGTTTGGACACCCCCTACgctttc                                                                                                                                         >  1:1135234/1‑142 (MQ=37)
                        aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑Taca                                                                                                                                 <  1:1111436/151‑1 (MQ=255)
                            cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑Tacaca‑c                                                                                                                             <  1:105294/150‑1 (MQ=255)
                               aGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑Cttt                                                                                                                          <  1:636746/150‑1 (MQ=255)
                                   aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAg                                                                                                                      <  2:523279/150‑1 (MQ=255)
                                     cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTCGAACCCCCCCTACGGTTGCGTA‑TAAACA‑CCTTTCaag                                                                                                                     >  2:903561/1‑147 (MQ=21)
                                        aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTTGAACCCCCCATACGTTGCCGGA‑TACACA‑CTTTTCAGGCGTg                                                                                                                 >  2:1006331/1‑150 (MQ=21)
                                          tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAACCCCCCATACGTTGCCGTA‑TAAACA‑CCTTTCAAGCggg                                                                                                                 >  2:916639/1‑146 (MQ=21)
                                           ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATTCGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTcc                                                                                                             <  1:385725/151‑1 (MQ=255)
                                             aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTCGAAACCCCCATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCtt                                                                                                           >  2:132977/1‑151 (MQ=255)
                                             aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGAGG‑AGGGGGGGTTCGAACCCCCCCCCCGTTTTCGTA‑TACACC‑CCTTTCCGGCGGGGTCCtt                                                                                                           >  1:138620/1‑151 (MQ=21)
                                               ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGAACCCCCCAAACGTTGCCGTA‑TAAACA‑CTTTTCCGGCGGGCTcc                                                                                                             >  2:122973/1‑147 (MQ=21)
                                               ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCa                                                                                                         <  2:224928/151‑1 (MQ=255)
                                                cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGATCGGAACCCCCCCAAGTTGCCGTA‑TTAACA‑CATTTCAGGCGTggccc                                                                                                             >  1:72622/1‑142 (MQ=18)
                                                 gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑c                                                                                                      <  1:299905/151‑1 (MQ=255)
                                                                          acagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGtt                                                                            <  1:839831/151‑1 (MQ=255)
                                                                             gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGAGATTCGAAACCCCCATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAAATGtttt                                                                          >  2:413779/1‑150 (MQ=255)
                                                                              aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGc                                                                        <  2:367955/151‑1 (MQ=255)
                                                                               tAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGAGA‑GGGGGGGGGTTGAAACCCCCCAACGTTTGCGTA‑TTCCCA‑CCTTTCTCGCGTGCTCCTTCCC‑ACC‑CTCCGGCACCCCCCCCAAATGTTTTTCt                                                                       >  1:104504/1‑151 (MQ=9)
                                                                                    tGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGAACCCCCCGTTCCTTTTCGTATTACA‑A‑CATTTTAGGGGGGGGCCTTCAA‑CCC‑CCCGGACAACTCACCAAATTGTTTTGTTA‑ACa                                                                  >  1:1103032/1‑150 (MQ=9)
                                                                                            aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGAGAGGGGGGGGGTTCG‑ACCCCCCATACGTTGCCGTA‑TACACA‑CTTTCCCGGCGGGCTCCCTCAG‑CCA‑CTCCGACCCCCCCCCCAAAATTTTTTCTA‑ACaaa                                                                >  1:859949/1‑144 (MQ=2)
                                                                                                 cTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCCGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCCCACCCAATTGTTTTTCTA‑CCAAACCTC‑AT‑gggggg                                                  >  2:436349/1‑151 (MQ=14)
                                                                                                   cGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTgg                                                  <  1:321567/149‑1 (MQ=255)
                                                                                                   cGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGc                                                 <  2:195364/150‑1 (MQ=255)
                                                                                                    gCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCaa                                               <  2:315006/151‑1 (MQ=255)
                                                                                                     cACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAAc                                              <  1:829692/151‑1 (MQ=255)
                                                                                                      aCTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGAACCCCCCATACGTTGCCGGA‑TACACA‑CTTTTCAGGCGGGCTCCTTCAG‑CCC‑CTCGGGCACCCCACCAAAATGTTTTTTTA‑CCAAAACTC‑CT‑gggggg                                                  >  2:1103660/1‑146 (MQ=11)
                                                                                                       cTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGT‑GGGGGGGGTTTGGACCCCCCCCACGTTGCCCTA‑TAAACACCTTTCCAGG‑GGGGTCCTTCCA‑CCC‑CCCCGGCACCCCACCACAATGTTTTTTTA‑CCAAAACCC‑CT‑gggggg                                                  >  1:866961/1‑145 (MQ=2)
                                                                                                       cTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGC‑GCGGGGGGCCCCAACCCCCGCGTTGTCGTCTTA‑TACAAA‑CTTTCTGCGCCGGCCCCCTCTC‑CCA‑ACCCGACCCCTCAACAAACCATTTTCCTA‑CCAAACCTC‑AC‑GGGGGGCa                                                >  1:63972/1‑147 (MQ=2)
                                                                                                          tCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGAACCCCCCATACGTTTCCGTA‑TACACA‑CTTTTCCAGCGTGCTCCCTCCG‑CCC‑CTCGGGCACCCCCCCCAAATTTTTTTCTT‑CCCAAACTC‑CT‑TGGGGGCAACgggg                                          >  2:368815/1‑150 (MQ=2)
                                                                                                           cGTACGAGCTCTTCTTTAAATATGGCGGGGGA‑GGGGGGGGATCGAACCCCGCAAAC‑TTCCCGTA‑TAAAGC‑CGTTTTCGGGGTGCTCCTTCCGCCCACCCCTGACACCAC‑CCAAATAGTTTT‑CTATCCAAACCAC‑TT‑GGTTGGCAACggg                                           >  1:386816/1‑149 (MQ=2)
                                                                                                            gTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAAACCCCCATACGGTGCCGGA‑TACAAA‑ATTTTCCGGCGGGGTCCTTCAA‑CCC‑CCCGGACACCCCACCAAAATGTTTTTCTA‑CCAAACCTC‑AT‑GGGGGGCAAAgggggg                                        >  2:331161/1‑148 (MQ=12)
                                                                                                            gTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGgcgc                                       <  2:957560/151‑1 (MQ=255)
                                                                                                            gTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGgcgc                                       <  1:903202/151‑1 (MQ=255)
                                                                                                                  gCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACtat                                 <  2:1150221/151‑1 (MQ=255)
                                                                                                                  gCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCCACACGTTGCCGTA‑TATACA‑CTTTCCAGGCGTGCTCCTTCAA‑CCA‑CCCCGAGACCCCACCAAAATTTTTTTCTA‑CCAAAACTC‑ATGGGGGGGCAAC‑GGGCGCGACtat                                 >  2:494250/1‑151 (MQ=14)
                                                                                                                   ctctTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTTTAAACCCCCACACCTTTCCGTA‑TACACA‑CATTCCCGGCGGGCTCCTTCAG‑CCA‑CCCCGCCCCCCCACCCAAATGTTTTTTTT‑ACAAACCCC‑AT‑GTGGGGGAACgggggg                                        >  1:652655/1‑141 (MQ=2)
                                                                                                                     cttcttTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAAACCCCCATAACTTTTCGTA‑TAAACA‑CTTTCCAGGGGGGCTCCCTCCA‑CCC‑CTCGGGCACCCCCCCAAAATTTTTTTCTA‑CCAAAACTCAAT‑GGGGGGAAAGgggggggac                                     >  1:822117/1‑140 (MQ=2)
                                                                                                                     cttcttTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACtata                                <  1:405558/149‑1 (MQ=255)
                                                                                                                     cttcttTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACtata                                <  2:175925/149‑1 (MQ=255)
                                                                                                                      ttcttTAAATATGGC‑GGTGA‑GGGGGTGATTCGATCCCCCTATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCTCCAAATAATTTTGCTG‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACtata                                <  1:435286/148‑1 (MQ=18)
                                                                                                                           tAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGtt                         <  2:615447/150‑1 (MQ=255)
                                                                                                                                attgC‑GGTGA‑GGAGGGTATTCGACCCCCCGATACGTAGCCGTA‑TACACC‑CTTTCCATGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTATTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAgggagtgggagt                    <  1:918850/147‑1 (MQ=21)
                                                                                                                                atGGC‑GGTGA‑GGGGGGGATTCGAAACCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCAACAAAATGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAAAGGGGCGCTACTATAGGGAGATGGa                      >  1:160474/1‑148 (MQ=255)
                                                                                                                                      ggTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg              <  2:47056/151‑1 (MQ=255)
                                                                                                                                         gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt            <  2:701339/150‑1 (MQ=255)
                                                                                                                                         gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           <  2:492908/151‑1 (MQ=255)
                                                                                                                                         gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           <  1:921936/151‑1 (MQ=255)
                                                                                                                                          a‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           >  1:956348/1‑150 (MQ=255)
                                                                                                                                          a‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          <  2:408062/151‑1 (MQ=255)
                                                                                                                                          a‑GGGGGAGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          <  1:413535/151‑1 (MQ=255)
                                                                                                                                               ggggATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa      <  1:282226/151‑1 (MQ=255)
                                                                                                                                                   aTTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  1:1114564/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                             
TCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTA‑TACACA‑CTTTCCAGGCGTGCTCCTTCAG‑CCA‑CTCGGACACCTCACCAAATTGTTTTGCTA‑CCAAACCTC‑AT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925755‑926050

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