Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC41.3% 178.3 / 28.9 92T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (46/8);  new base C (0/38);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 5.07e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376361‑1376648
                                                                                                                                                      |                                                                                                                                               
ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCt                                                                                                                                             <  1:553320/149‑1 (MQ=255)
ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCg                                                                                                                                           >  3:591385/1‑151 (MQ=255)
 tgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAATGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGc                                                                                                                                          >  1:906950/2‑151 (MQ=255)
  gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGcc                                                                                                                                         <  2:149850/151‑1 (MQ=255)
  gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGcc                                                                                                                                         >  3:205840/1‑151 (MQ=255)
     gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGAt                                                                                                                                      <  3:397121/151‑1 (MQ=255)
      aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATg                                                                                                                                     >  1:864589/1‑151 (MQ=255)
        gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGcc                                                                                                                                   <  3:79156/151‑1 (MQ=255)
        gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGcc                                                                                                                                   <  1:286491/151‑1 (MQ=255)
          ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCccg                                                                                                                                 >  1:103219/1‑151 (MQ=255)
             cgccgcGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCccgccg                                                                                                                              >  3:1273010/1‑151 (MQ=255)
             cgccgcGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCccgcc                                                                                                                               <  1:583857/150‑1 (MQ=255)
                   gTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACAt                                                                                                                        >  2:136560/1‑151 (MQ=255)
                       ccTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                    >  1:546847/1‑151 (MQ=255)
                          ggACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                  >  3:864658/1‑150 (MQ=255)
                          ggACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                  >  3:350499/1‑150 (MQ=255)
                          ggACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                  >  1:1227617/1‑150 (MQ=255)
                            aCGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                                >  2:261959/1‑150 (MQ=255)
                              ggCGCCCGAAAAGGGGG‑GAGTTTTCCCCTGCGCGTACGGGAAAGGGCCGCCGCGCGTCCCGGACCG‑GT‑GGTTTACCCCGACCCTTTAAC‑CCGCTATCAGGTTA‑ATATGACCAAC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                             <  1:421954/151‑1 (MQ=255)
                              ggCGCCCGAAAAAGGGT‑GAGTGTTCCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTCCGGGAACC‑GG‑GGTGTATCTCGACCCTTAAAA‑CCGCTATCAGGATA‑ATATGATAAAC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGaaaa                                                                                                               <  3:703457/149‑1 (MQ=255)
                               gCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                             >  3:149775/1‑150 (MQ=255)
                                cGACCGAAAAATGGT‑GAGTGATCCCCGGCGCGTCCTGAAACAGGTCGCCGCGCGTCCCGGACCC‑GT‑GGTTTATCTCGACCCTTAAAC‑CCCCTATCAGGATC‑ATATAATCAAC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                            <  3:651319/150‑1 (MQ=255)
                                   ccgaaaaggggt‑gggtgTTCCCCTGCGCGTACTGAAAGGGGTCGCCGCGCGTCCCGGACCC‑GT‑GGTGTATCCCGACCCGTTAAC‑CCGCTAACAGGATC‑ATATAACCACC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                        <  1:871098/142‑1 (MQ=255)
                                    cggaaaaTTGTGGAGTGCTCCCCTGGGCGTAGGGAAACAGCGCGCCGGGGGAACCGAAGCG‑GT‑GGTTTATCTGGAGCCGTAAAA‑CCGCTATCAGGATC‑ATATGA‑CAACGCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                            <  3:1068389/144‑1 (MQ=255)
                                     gCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                       >  3:772572/1‑150 (MQ=255)
                                     gCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                       >  1:1403187/1‑150 (MQ=255)
                                       aaatTGGT‑GATGTATCCCCTGGGCTTCCGGAAACGGGGCGCCGCGCGCCCCGGACCC‑GT‑GTTGTTTCTCGATCCTTTAAA‑CCGTTACCGGGATC‑ATATAATAAAC‑CCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                           <  1:754955/144‑1 (MQ=255)
                                       aaaaTGGT‑GAGTGATCCCCTGCGCGTACTGGAACGGGTCGCCGCGCGTACCGGACCC‑GT‑GGTGTATCTCGAACCGTTAAC‑CCGCTATCAGGATA‑ATATGACCAAC‑CCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                      <  3:1244607/149‑1 (MQ=255)
                                       aaaaGGGG‑GAGTGTTTCCCTGCGCGTACGGAAACGGGCCGCCGCGCGTACCGAAACG‑GT‑GGTGTACCTCAACCCGTTAAC‑CCGCTATAAGGATA‑ATATAAAAAAC‑CCCCCCCCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                     <  1:530781/150‑1 (MQ=255)
                                        aaaTGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                    >  4:132462/1‑150 (MQ=255)
                                            ggT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                >  1:150329/1‑150 (MQ=255)
                                            ggT‑GAGTGATTCCCTGCGCGTACTGAAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGACCCGTTAAC‑CCGCTATCAGGATC‑ATATGATCAAC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                <  3:1201923/150‑1 (MQ=255)
                                             ggggagtgaTCCCCGGCGCTTACTGGACCGGGCCGCCGCGCTTACCGGACCC‑GT‑GTTTTATCTC‑AACCCTTAAACCCCCTATCAGAATC‑ATAGAACCACC‑CCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                <  2:361726/147‑1 (MQ=255)
                                                 agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                              >  3:1364145/1‑148 (MQ=255)
                                                     aTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                       >  1:1273826/1‑151 (MQ=255)
                                                      ttCCCTGCGCTTAGTGAAACAGGCCGCCGCGCGTCCCGGAACG‑GT‑GGTTTATCTCGACCCTTAAAC‑CCGCTATAGGGATC‑ATATGATCAAC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                           <  1:486096/146‑1 (MQ=255)
                                                      ttCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                      >  1:258337/1‑151 (MQ=255)
                                                        cccGGCGCGTACTGAAACGGGTCGCCGCGCGTCCCGGACCG‑GG‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGATTA‑ATATAATCACC‑CCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                     <  4:437034/150‑1 (MQ=255)
                                                          cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                  >  3:52634/1‑151 (MQ=255)
                                                            gcgcGTAGGGGAACGGGCCCCCGCGCGTCCCGGACCG‑GT‑GGTTTATCTCGACCCGTTAAC‑CCGCTATCAGGATC‑ATATGATCAAC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTa                                                                                        <  1:718933/143‑1 (MQ=255)
                                                            gcgcGTACTGAAAGAGGTCGCCGCGCGTACCGGACCC‑GT‑GTTGTATCCCGAACCTTAAAC‑CCTCTACAAGGACA‑AAATGACCAAC‑CCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                                 <  1:480435/150‑1 (MQ=255)
                                                               cgtccggaacctgttCGCCGCGC‑TACCCGAACCTGT‑GGTTTACCTC‑AACCCTTAAC‑CCCGTAACCAGATCAATTTAATAACC‑ACCCCCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                              <  4:153053/136‑1 (MQ=255)
                                                                  aCTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                          >  1:413346/1‑151 (MQ=255)
                                                                   cggaaaaagtccgccgcGCGTCCCGGAACC‑GT‑GGTCTGTCCCGACCCTTTAAG‑CCGCTATCAGATTA‑ATTTATTCACC‑CCCCCCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         <  1:1114982/140‑1 (MQ=255)
                                                                    gggAACAGGTCGCCGCGCTAACGGGACCC‑GT‑GTTGTCTCCCAACCCGTAAAC‑CCCCAAGAAGTATA‑ATATAATCAAC‑CCCCCCGCAAAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                        <  4:54453/150‑1 (MQ=255)
                                                                     ggAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                       >  2:340574/1‑151 (MQ=255)
                                                                     gaaaCGGGTCGCCGCGCGTACCGGACCG‑GT‑GGTGTATCTCGAACCGTAAAA‑CCGCTATCAGGATC‑ATATGATAAAC‑CCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         <  2:22723/147‑1 (MQ=255)
                                                                     gaaaCAGGTCGCCGCGCGTACCGGAACC‑GG‑GGTGTATCTCAACCCGTTAAC‑CCGCTATCAGGATC‑ATATGATCAAC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                          <  2:457031/146‑1 (MQ=255)
                                                                       aaCGGGCCGCCGCGCGTACCGGAACG‑GT‑GGTGTATCTCAACCCTTTAAC‑CCGCTACCAGGATC‑ATATGATCACC‑CCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                     <  1:746826/151‑1 (MQ=255)
                                                                         cAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                    >  3:933369/1‑150 (MQ=255)
                                                                             tCGCGGCGCGTACCGAAACC‑GT‑GGTGTACCTCGACCCTTTAAA‑CCGCTATAAGGATC‑ATATGACCAAC‑CCCCTCGCCGTTCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                <  3:251053/150‑1 (MQ=255)
                                                                             tCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                               >  3:1399904/1‑151 (MQ=255)
                                                                              cgccgcGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                              >  1:1255444/1‑151 (MQ=255)
                                                                                ccgcGCTTCCGGGACCC‑GG‑GGTGTATCTC‑AAACGTTAAACCCGCTATCAGGATC‑ATATGATAAAC‑CCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                             <  3:519825/150‑1 (MQ=255)
                                                                                ccgcGCGTCCCGGACCC‑GT‑GGTGTATCTCGACCCGTTAAC‑CCGCTATCAGAATA‑ATATGATCAAC‑CCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                <  4:13148/147‑1 (MQ=255)
                                                                                ccgcGCGTCCCGGACCC‑GG‑GTTGTATCTC‑AACCCTTAAACCCGCTATCAGGATA‑ATATGATCAAC‑CCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                              <  1:308785/149‑1 (MQ=255)
                                                                                ccgcGCGCCCCGGACCG‑GT‑GGTGTATCCCGAACCGTAAAC‑CCGCTAAAAGGATA‑ATATGACAACC‑CCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                             <  3:127073/150‑1 (MQ=255)
                                                                                ccccgcgTACCGGACCG‑GG‑GGTTTTCCCCGACCCTTAAAC‑CCGCTATAAGGATC‑ATGTGAACACC‑CCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                               <  1:731098/145‑1 (MQ=255)
                                                                                 cgcgcgTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                            >  2:400840/1‑150 (MQ=255)
                                                                                  gcgcgTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGACTATCGGCGATTTCTATCACATGAATATCga                                                            >  3:473798/1‑149 (MQ=255)
                                                                                    gcgTCCCGGACCC‑GT‑GGTGTTTCCGGACCCTTTAAC‑CCGCTATCAGAATC‑TTATGATCACC‑CCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         <  2:313381/150‑1 (MQ=255)
                                                                                      gccccGAAACC‑GG‑GGTGTTTCCCAAACCGTTAAA‑CCGATATCAGAATC‑ATATGATAAAC‑CCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                      <  1:293113/148‑1 (MQ=255)
                                                                                         ccGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                    >  3:1374699/1‑150 (MQ=255)
                                                                                          cGGAAAC‑GTGGGTGTCTCTCAAACCTTAACA‑CCGCTATCAGGATA‑ATATTATCACC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                   <  1:52668/151‑1 (MQ=255)
                                                                                           ggAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                  >  3:976636/1‑150 (MQ=255)
                                                                                            ggaccggt‑ggtgTATCTCGACCCGTAAAC‑CCGCTAACAGGATC‑ATATGATCAAC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                  <  3:473548/144‑1 (MQ=255)
                                                                                            ggaccggt‑ggtgTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATAATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                  <  2:548287/144‑1 (MQ=255)
                                                                                            gAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                >  1:1393086/1‑151 (MQ=255)
                                                                                              accggt‑ggtgTACCTCGACCCTTAAAC‑CCGCTATCAGAATA‑ATATAATCAAC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                <  4:330373/146‑1 (MQ=255)
                                                                                                c‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                            >  2:173219/1‑151 (MQ=255)
                                                                                                   t‑ggtgTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                           >  3:302573/1‑150 (MQ=255)
                                                                                                      gtgtATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                        >  4:545175/1‑151 (MQ=255)
                                                                                                       tgtATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                       >  1:1259013/1‑151 (MQ=255)
                                                                                                       tgtATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                       >  2:181259/1‑151 (MQ=255)
                                                                                                          aTCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                    >  3:1124702/1‑151 (MQ=255)
                                                                                                          aTCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                    >  3:352922/1‑151 (MQ=255)
                                                                                                               cccccctttaccccccTATCAGTTTC‑ATATCATCCCC‑CCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                <  1:933245/136‑1 (MQ=255)
                                                                                                                        aacccGTTATAAGGATC‑ATATAACCACC‑CCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACaa                           <  1:660998/144‑1 (MQ=255)
                                                                                                                        aa‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                      >  1:1233110/1‑151 (MQ=255)
                                                                                                                           ccGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                    >  3:186232/1‑151 (MQ=255)
                                                                                                                           ccGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                    >  2:404186/1‑151 (MQ=255)
                                                                                                                            cgttaacAGGATA‑ATTTAACCACC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGa                     <  3:56925/143‑1 (MQ=255)
                                                                                                                            cGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                   >  3:49971/1‑151 (MQ=255)
                                                                                                                            cGCTATAAGGTTC‑ATATCATCACC‑CCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                      <  1:752836/148‑1 (MQ=255)
                                                                                                                            cGCTACCAGAATC‑AAATGAAAAAC‑CCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                   <  3:580085/151‑1 (MQ=255)
                                                                                                                             gCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                  >  1:1302852/1‑151 (MQ=255)
                                                                                                                               tATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                >  1:87383/1‑151 (MQ=255)
                                                                                                                                aTCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt               >  3:962564/1‑151 (MQ=255)
                                                                                                                                            atGATCACC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg     <  3:258357/150‑1 (MQ=255)
                                                                                                                                            atGATCAAC‑CCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 <  4:277501/138‑1 (MQ=255)
                                                                                                                                              gATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  3:39930/1‑151 (MQ=255)
                                                                                                                                               accAAC‑CCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCa   <  3:413411/147‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
CCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGT‑GAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAAC‑ACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376361‑1376648

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: ‑