Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913925,9000AG45.6% 105.7 / 17.4 69noncoding (72/88 nt)serWtRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (7/30);  new base G (31/0);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 7.44e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925751‑926050
                                                                                                                                                       |                                                                                                                                                          
gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGAt                                                                                                                                                           <  1:924063/151‑1 (MQ=255)
  gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTc                                                                                                                                                         <  1:875089/151‑1 (MQ=255)
  gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTc                                                                                                                                                         <  1:317331/151‑1 (MQ=255)
    tCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGtttg                                                                                                                                                        >  2:166757/1‑148 (MQ=255)
     cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGaa                                                                                                                                                      <  2:243969/151‑1 (MQ=255)
       cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCggaac                                                                                                                                                    >  1:721122/1‑147 (MQ=255)
           cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCg                                                                                                                                                <  2:313536/151‑1 (MQ=255)
           cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCg                                                                                                                                                <  1:756039/151‑1 (MQ=255)
                aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGGAAcccccaca                                                                                                                                              >  2:851883/1‑146 (MQ=255)
                aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGGACCCCCCAt                                                                                                                                              >  1:1045330/1‑148 (MQ=255)
                  gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGA‑GGGGGGGATTGGAcccccccgacact                                                                                                                                          >  2:351066/1‑144 (MQ=37)
                  gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑AAGtt                                                                                                                                        <  2:269272/151‑1 (MQ=255)
                       ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAAACCCCGAT‑ACGTTg                                                                                                                                       >  2:498727/1‑147 (MQ=255)
                        tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGt                                                                                                                                   <  1:699045/150‑1 (MQ=255)
                         aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTTGAACCCCCCCC‑AAGTTGCCgca                                                                                                                                  >  2:606391/1‑148 (MQ=37)
                           aaaTTCACAGGTAACATACTCCACCGGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGtata                                                                                                                                <  2:767189/150‑1 (MQ=255)
                                  caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTTGAAACCCCCAT‑ACGTTGCCGGATACACACctt                                                                                                                        >  1:750096/1‑151 (MQ=25)
                                        aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTTGAACCCCCCAC‑ACGTTGCCGTATACACACTTTcc                                                                                                                      >  2:914615/1‑147 (MQ=255)
                                          aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAAACCCCGGTAAGGTGGCCGT‑TACACACCTTTCAGGCg                                                                                                                 >  1:257798/1‑150 (MQ=21)
                                          aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGGGA‑GGGGGGGGGTTGAACCCCCCCT‑ACCTTGCCCGATAAACACctt                                                                                                                        >  2:819476/1‑143 (MQ=21)
                                             cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGG‑GGGGGGGGTTTGAACCCCCCCT‑ACGTTGCCCTATACACACCTTTCAgg                                                                                                                   >  1:853255/1‑145 (MQ=21)
                                              tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGGACCCCCGCC‑AAGGTGCCCTATAAACACTTTCCAGGCg                                                                                                                 >  2:760727/1‑146 (MQ=21)
                                              tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGGTCGAACCCCCCCC‑ACGTTTGCGTATAAACACCTTCCatg                                                                                                                   >  2:361734/1‑142 (MQ=21)
                                                 aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAACCCCCCAT‑ACGTTGGCGTATACACACCTTCCAGGCGTGcccc                                                                                                           >  1:286294/1‑149 (MQ=21)
                                                      cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAACCCCCCCA‑ACAGTGCCGGATACACACCTTTCCGGCGGGCTCCTTCCG‑cc                                                                                                   >  2:830178/1‑151 (MQ=18)
                                                            aTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGGGTGCTCCTTCAG‑CCCCTCgg                                                                                             >  1:928779/1‑151 (MQ=255)
                                                                       cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAACCCCCCAT‑ACGTTGCCGTATACACACCTTTCCGGCGTGCTCCTTCCC‑CCACCCGGGCA‑CCCCAcc                                                                                  >  2:861269/1‑150 (MQ=17)
                                                                        gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCaa                                                                                <  1:1113637/151‑1 (MQ=255)
                                                                         tACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGGTA‑AGGGGGGTTTAGTAAACCCGAT‑ACGTTGCCGTATTAACACTTTCCGGGCGCGATCCTTTCG‑CCACCCGGCCA‑CCCCACCaa                                                                                >  1:670023/1‑150 (MQ=17)
                                                                          acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAAt                                                                              <  2:13948/151‑1 (MQ=255)
                                                                              acagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGtt                                                                          <  2:504234/151‑1 (MQ=255)
                                                                                agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGtttt                                                                        <  2:94468/151‑1 (MQ=255)
                                                                                agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGGGA‑GGGGGGGATTTGAACCCCCCAC‑ACGTTGCCGTATACAAACTTTTCCGGCGTGGTCCTTCAG‑CCACCCGGGCC‑CCTCACCAAATTGtttt                                                                        >  1:175244/1‑151 (MQ=14)
                                                                                 gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGAAACCCCCAT‑ACGTTGCCGTATACACACTTTTCAGGCGTGCTCCTTCAG‑CCACCCGGACA‑CCCCCCCAAATTGtttt                                                                        >  1:112485/1‑150 (MQ=14)
                                                                                   tAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAACCCCCCAT‑AAGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACCCGGGCA‑CCTCACCAAAATGtttt                                                                        >  1:257723/1‑148 (MQ=14)
                                                                                    aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCt                                                                     <  1:914637/150‑1 (MQ=255)
                                                                                    aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTa                                                                    <  2:512746/151‑1 (MQ=255)
                                                                                          cGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTCGAACCCCCCAT‑ACGTTTCCGTATACACACTTTTCCGGCGTGCTCCTTCAC‑CCCCCCCGACA‑ACTCACCAAATTGTTTTTCTACCcaa                                                               >  2:870700/1‑150 (MQ=14)
                                                                                           gAAAAAAAAGCCCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGGGG‑GGGGGGGGTTCGAACCCCCCAC‑ACGGTGCCGGATAAACACTTTTCCGGCGGGCTCCTTCAG‑CCCCCCCGACA‑CCCCACCAAAAAGTTTTTTTTCCaaa                                                               >  1:171707/1‑149 (MQ=2)
                                                                                              aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGGGA‑GGGGGGGGTTCGAACCCCCCGT‑ACGTTTGCGGTTACAAACCTTCCACGCGTGCTCCTTCAG‑CCACTCGGGCA‑CCCCACCCAAATGTTTTGCTTCCAAAACTc                                                           >  1:945702/1‑150 (MQ=14)
                                                                                               aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGGAACCCCCAG‑TCGTTTCCGTATAAACACTTTTCAGGCGGGGTCCCTCCG‑CCCCCCGGGCA‑CCCCCCCCAAATTTTTTTTTACCaa                                                                >  2:968927/1‑144 (MQ=2)
                                                                                                aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAACCCCCGAT‑ACGTTGCCGTATAAACACTTTTCCGGCGGGCCCCTTCAG‑CCACCCCGACA‑CCTCACCAAAATGTTTTTCTTCCAAAcc                                                             >  1:333286/1‑146 (MQ=14)
                                                                                                aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑g                                                       <  1:361742/151‑1 (MQ=255)
                                                                                                 aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGGACCCCCCAT‑ACGTTTTCGTATACACACTTTTCACGCGTGCTCCCTCCC‑CCCCCCGCGCA‑CCCCACCCAAATGTTTTGCTACCAAAACTCa                                                          >  2:88062/1‑148 (MQ=9)
                                                                                                  aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑ggg                                                     <  2:747730/151‑1 (MQ=255)
                                                                                                    gCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGAACCCCCCAT‑AAGGTGGCGTATAAACACTTTCCAGGCGGGCTCCTTCAG‑CCACCCCGACA‑CCCCACCAAAATTTTTTTCTACCAAACCTCAT‑ggggg                                                   >  1:766588/1‑149 (MQ=14)
                                                                                                       cGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGc                                                 <  1:1775/150‑1 (MQ=255)
                                                                                                          aCTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGCGG‑GGGGGGGGTTCGAACCCCCCCA‑ACGGTTCCCTTTACAAACCCTCCAAGCGTGCTTCCTCCC‑CCACCCCGGCA‑CCTCACCCAAATGCTTTTCTCCCAAAACCTAT‑GGGTGGCaa                                               >  1:535022/1‑149 (MQ=2)
                                                                                                              tCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACgggg                                          <  2:151568/150‑1 (MQ=255)
                                                                                                              tCGTACGAGCTCTTCTTTAAATATGGC‑GGGTA‑GGGGGGGGATCGGAACCCCCAT‑ACGGTTCCGTATACACACTTTTCAGGGGGGGTCCTTCAGCCCCCCCGGCAA‑CCCCCCCAAATTTTTTTG‑TTCCACAACTCCTGGGGGGGCAACgg                                            >  2:529770/1‑149 (MQ=2)
                                                                                                               cGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGgc                                         <  1:396610/150‑1 (MQ=255)
                                                                                                                  aCGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGGTTTGGAACCCCCAC‑ACCGTTCCGTATAAAAACATTTCCGGGGGGCTCCCTCCG‑CCCCCCCGGCACCCCCACCAAATTTTTTTCTTCCCAAACCCATT‑GGGGGGAAAGggggg                                         >  2:130928/1‑147 (MQ=2)
                                                                                                                  aCGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGct                                      <  2:1052389/150‑1 (MQ=255)
                                                                                                                  aCGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGct                                      <  1:789817/150‑1 (MQ=255)
                                                                                                                     aGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAAACCCCCAT‑ACGTTGCCGTATAAACACCTTCCAGGCGTGCTCCTTCAA‑CCACTCCGGCA‑CCCCCCCAAATTTTTTTTTTACCAAACCTCAT‑TGGGGGCAACGGGGGGccac                                    >  1:982/1‑149 (MQ=14)
                                                                                                                      gCTCTTCTTTAAATATGGC‑GGTGG‑GGGGGGGGATTGAAACCCCCCC‑ACCTTTGCGTATTGACCCTTTTCCGGCGGGGTCCTTCCGCCCACCGGGCCA‑CCCCCCCCAA‑TGTTTTGCTAACCAACCCCAC‑GGGGGGGGACGGGGGGccac                                    >  2:69227/1‑148 (MQ=255)
                                                                                                                         cttcttTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACtata                                <  2:633592/149‑1 (MQ=255)
                                                                                                                            cttTAAATATGGC‑GGTGA‑GGGGGGGGGTCGGACCCCCCCT‑ACGTTTCCGTATACACACTTTTCAGGCGGGCTCCTTCCG‑CCCCCCCGGCA‑CCCCCCCCAAATGTTTTTTTACCAAAACCCCT‑CGGGGGGAACGGGGGGCCACTTTAggggg                           >  2:722507/1‑149 (MQ=2)
                                                                                                                              ttAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAg                           <  1:606409/149‑1 (MQ=255)
                                                                                                                              ttAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGt                          <  2:497927/150‑1 (MQ=255)
                                                                                                                                        c‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               >  1:776408/1‑151 (MQ=255)
                                                                                                                                            tGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt            >  2:1036764/1‑151 (MQ=255)
                                                                                                                                             gA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           <  2:66838/151‑1 (MQ=255)
                                                                                                                                              a‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          <  1:1040580/151‑1 (MQ=255)
                                                                                                                                                gggggggATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          >  1:40076/1‑150 (MQ=255)
                                                                                                                                                gggggggATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          <  1:851902/150‑1 (MQ=255)
                                                                                                                                                     ggATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAagaaga    <  2:1136905/151‑1 (MQ=255)
                                                                                                                                                      gATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAagaaga    >  2:965308/1‑150 (MQ=255)
                                                                                                                                                      gATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa   <  1:830198/151‑1 (MQ=255)
                                                                                                                                                       aTTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  >  2:1152024/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                          
GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGA‑GGGGGGGATTCGAACCCCCGAT‑ACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAG‑CCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925751‑926050

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