Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC43.5% 111.8 / 12.4 62T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (31/4);  new base C (0/27);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 2.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655
                                                                                                                                                           |                                                                                                                                                      
gtttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACATGTCGCCGCGCGTACCGGAT‑A‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAacac                                                                                                                                                       <  1:282733/149‑1 (MQ=255)
 cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACAcc                                                                                                                                                      >  1:552235/1‑151 (MQ=255)
  ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACAccc                                                                                                                                                     >  2:816017/1‑151 (MQ=255)
   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCt                                                                                                                                                    >  1:453431/1‑151 (MQ=255)
   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCt                                                                                                                                                    >  1:446461/1‑151 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCg                                                                                                                                                  <  1:929970/151‑1 (MQ=255)
      cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                                 <  2:382406/151‑1 (MQ=255)
               ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCccg                                                                                                                                        <  1:323576/151‑1 (MQ=255)
                 tCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                                      >  2:957034/1‑151 (MQ=255)
                     cgcgTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                                   >  2:1109170/1‑150 (MQ=255)
                         tAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtc                                                                                                                              >  1:3333/1‑151 (MQ=255)
                            ccTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           >  1:810612/1‑151 (MQ=255)
                                   ggCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                    >  2:1048165/1‑151 (MQ=255)
                                        ccGCAAAAGTGGTTAGAGTTCCCCGGCTCCTACGGTAAACGGCCGCCGCGCGT‑CCGGAA‑CGCTTGT‑TGTATC‑CCAACCCTTAAA‑CCGTCTACGAGGCAAATATGCAATCCCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                                 <  1:581595/150‑1 (MQ=255)
                                        ccGAAAAA‑TGGTGAGTGTTCCCCGGCGCGTACGGGAACAGTTCCCCGCGCGTACCGAAC‑C‑CGTGG‑TGTATCTCAACCCTTTAAA‑CCG‑CTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                                <  2:883595/150‑1 (MQ=255)
                                         cGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                               >  2:293826/1‑150 (MQ=255)
                                            aaaa‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                           >  1:768043/1‑151 (MQ=255)
                                            aaaa‑GGGTGAGTGTTCCCCTGCGCGTACGGAAACGGGTCGCCGCGCGTCCCGGAC‑C‑CGTGG‑TGTATCTCGACCCTTTAAC‑CCG‑TTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                             <  2:494454/149‑1 (MQ=255)
                                             aaa‑GGGGGAGTGATCCCCGGCGCGTACTGAACCGGGTCCCCGCGCGTACCGGAA‑C‑GGGGG‑TGTATCTCGACCCGTAAAC‑CCG‑CTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                          <  2:1067830/151‑1 (MQ=255)
                                                       gtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     >  2:1012051/1‑147 (MQ=255)
                                                           tccccTGCGCGTACTGGAACGGGTCGCCGCGCGTCCGGGAC‑C‑CGGGG‑TGTATCTCGACCCTTAAAC‑CCG‑CTACAAGGATCAAATAATAAACCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                              <  2:840604/148‑1 (MQ=255)
                                                             cccgccgcgTACGGGACCGGGCCCCCGCGCGTACCGGAC‑C‑GGGGG‑TTTACCTCGAACCGTTAAA‑CCG‑TTATCAGGACAAAATGATAACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                             <  1:145853/144‑1 (MQ=255)
                                                              ccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                          >  2:354114/1‑151 (MQ=255)
                                                                 gcgcGTACTAAAACGGGCCCCCGCGCGCCCCGAAA‑C‑GGGGG‑TTTATCCCAACCTGTTAAC‑CCT‑CTAACGGGATAAAATACACAACCCCCTCGGCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                              <  1:116017/144‑1 (MQ=255)
                                                                     gTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                   >  1:249019/1‑151 (MQ=255)
                                                                      tacggaaagaGGCCGCCGCGCGT‑CCCGAA‑C‑CGGGGTTTTATCTAGAACCTTTAAC‑CCC‑TTACAGGGATAAAATAATCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                      <  1:15691/138‑1 (MQ=255)
                                                                          ggacaaGGTCCCGGCGCGTACGGAAC‑C‑CGGGT‑TTTATTACGAACCTTAAAA‑CCC‑TTATCAGAATAAAATATTCACCCCCCGCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                                        <  2:106141/136‑1 (MQ=255)
                                                                            accgggccgccgcGCGT‑CCGGAA‑C‑CGGGGTTTTATCTCGACCCTTAAAC‑CCG‑CTACAAGAATAAAATGATAAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                            <  1:336094/144‑1 (MQ=255)
                                                                             acgggTCGCCGCGCG‑ACCGGAACC‑CGGGG‑TGTA‑CCCCAACCCTTAAACCCG‑TTATCAGAATAATATAAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                               <  1:259881/144‑1 (MQ=255)
                                                                                ggTCCCCGCGGGTACGGGAA‑C‑GGGGG‑TGTATCTCGACCCGTTAAC‑CCG‑CTACAAGAATCATATGAACAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                             <  1:35613/146‑1 (MQ=255)
                                                                                    gccgcGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    >  1:200251/1‑151 (MQ=255)
                                                                                    gccgcGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    >  2:1040992/1‑151 (MQ=255)
                                                                                     ccgcGGGTCCCGAAC‑C‑GGTGG‑TTTTTCCCAACCTGTAAAA‑CCG‑CTATAGGGATAATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   <  1:406997/151‑1 (MQ=255)
                                                                                     ccgcGCGTCCCGGAC‑C‑CGTGG‑TTTATCTCGACCCTTAAAC‑CCG‑CTATCAGGATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                                     <  2:90448/149‑1 (MQ=255)
                                                                                         gcgTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  2:39326/1‑150 (MQ=255)
                                                                                          cgTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  2:219000/1‑149 (MQ=255)
                                                                                          cgTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  2:565268/1‑149 (MQ=255)
                                                                                             aCCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            >  2:418623/1‑150 (MQ=255)
                                                                                                 gAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                       >  1:1019727/1‑151 (MQ=255)
                                                                                                   acc‑cGTGG‑TGTCTCTCGACCCGTTAAC‑CCG‑CTATCAGAATCATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                           <  2:1073713/144‑1 (MQ=255)
                                                                                                        gtgg‑ttTATCTCGACCCTTTAAC‑CCG‑CAAAAAGGATAATATTATAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                     <  2:108698/148‑1 (MQ=255)
                                                                                                        gtgg‑tgTATCCCGAACCGTAAAC‑CCG‑CTATCAGGATCATATAACAAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                   <  2:1019809/150‑1 (MQ=255)
                                                                                                          gg‑tgTATCTCGACCCGTAAAC‑CCG‑CTATCAGGATCATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                                <  2:458652/151‑1 (MQ=255)
                                                                                                            ggtttacccagaCCCTTAAACCCG‑CTATCAGGATAATTTGAACACCCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                    <  2:977726/136‑1 (MQ=255)
                                                                                                               tATCTC‑AACCCTTAAACCCG‑CTATAAGAATCATATAATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                              <  2:21516/149‑1 (MQ=255)
                                                                                                                  cccaaaCCGTTAAC‑CCG‑CTATCAGGATAATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                          <  2:623568/146‑1 (MQ=255)
                                                                                                                      aaCCGTAAAC‑CCG‑CTATCAGGATCATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                       <  2:442079/149‑1 (MQ=255)
                                                                                                                        ccGTAAAC‑CCG‑TTATCAGGATAATATGACCAACCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                                   <  2:15110/151‑1 (MQ=255)
                                                                                                                            tAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAActg                               >  1:669903/1‑149 (MQ=255)
                                                                                                                                 ccG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                           >  1:1050509/1‑151 (MQ=255)
                                                                                                                                 ccG‑CTATCAGAATAAAAGAATAACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             <  1:560523/149‑1 (MQ=255)
                                                                                                                                   g‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                         >  2:410898/1‑151 (MQ=255)
                                                                                                                                     cTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                        >  2:1089737/1‑151 (MQ=255)
                                                                                                                                       aTCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt                      >  2:1044512/1‑151 (MQ=255)
                                                                                                                                       aTCAGGATAATATGAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                        <  2:329945/149‑1 (MQ=255)
                                                                                                                                           ggacaatataaaaaaCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt                  <  1:1090332/138‑1 (MQ=255)
                                                                                                                                           ggATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt                  >  1:827976/1‑151 (MQ=255)
                                                                                                                                              tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt               >  1:1132444/1‑151 (MQ=255)
                                                                                                                                                 tataaccAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                 <  2:453459/140‑1 (MQ=255)
                                                                                                                                                  atGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc           >  1:719503/1‑151 (MQ=255)
                                                                                                                                                  atGACCACCCCCCCCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                    <  2:1018537/142‑1 (MQ=255)
                                                                                                                                                   tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc           >  1:1062395/1‑150 (MQ=255)
                                                                                                                                                      tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat       >  1:844966/1‑151 (MQ=255)
                                                                                                                                                         acacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGc    >  2:142369/1‑151 (MQ=255)
                                                                                                                                                           acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  2:65516/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                      
GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAA‑TGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑C‑CGTGG‑TGTATCTCGAACCGTTAAA‑CCG‑CTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655

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