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

CCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376361‑1376648
                                                                                                                                                      |                                                                                                                                               
ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                             <  1:553320/149‑1 (MQ=255)
ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCg                                                                                                                                           >  3:591385/1‑151 (MQ=255)
 tgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAATGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                          >  1:906950/2‑151 (MQ=255)
  gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                         <  2:149850/151‑1 (MQ=255)
  gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                         >  3:205840/1‑151 (MQ=255)
     gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                      <  3:397121/151‑1 (MQ=255)
      aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                     >  1:864589/1‑151 (MQ=255)
        gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGcc                                                                                                                                   <  3:79156/151‑1 (MQ=255)
        gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGcc                                                                                                                                   <  1:286491/151‑1 (MQ=255)
          ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccg                                                                                                                                 >  1:103219/1‑151 (MQ=255)
             cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgccg                                                                                                                              >  3:1273010/1‑151 (MQ=255)
             cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                               <  1:583857/150‑1 (MQ=255)
                   gTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAt                                                                                                                        >  2:136560/1‑151 (MQ=255)
                       ccTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                    >  1:546847/1‑151 (MQ=255)
                          ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                  >  3:864658/1‑150 (MQ=255)
                          ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                  >  3:350499/1‑150 (MQ=255)
                          ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                  >  1:1227617/1‑150 (MQ=255)
                            aCGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                                >  2:261959/1‑150 (MQ=255)
                              ggCGCCCGAAAAGGGGGGAGTTTTCCCCTGCGCGTACGGGAAAGGGCCGCCGCGCGTCCCGGACCGGTGGTTTACCCCGACCCTTTAACCCGCTATCAGGTTAATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                             <  1:421954/151‑1 (MQ=255)
                              ggCGCCCGAAAAAGGGTGAGTGTTCCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTCCGGGAACCGGGGTGTATCTCGACCCTTAAAACCGCTATCAGGATAATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGaaaa                                                                                                               <  3:703457/149‑1 (MQ=255)
                               gCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                             >  3:149775/1‑150 (MQ=255)
                                cGACCGAAAAATGGTGAGTGATCCCCGGCGCGTCCTGAAACAGGTCGCCGCGCGTCCCGGACCCGTGGTTTATCTCGACCCTTAAACCCCCTATCAGGATCATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                            <  3:651319/150‑1 (MQ=255)
                                   ccgaaaaggggtgggtgTTCCCCTGCGCGTACTGAAAGGGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCCCGACCCGTTAACCCGCTAACAGGATCATATAACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                        <  1:871098/142‑1 (MQ=255)
                                    cggaaaaTTGTGGAGTGCTCCCCTGGGCGTAGGGAAACAGCGCGCCGGGGGAACCGAAGCGGTGGTTTATCTGGAGCCGTAAAACCGCTATCAGGATCATATGACAACGCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                            <  3:1068389/144‑1 (MQ=255)
                                     gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                       >  3:772572/1‑150 (MQ=255)
                                     gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                       >  1:1403187/1‑150 (MQ=255)
                                       aaatTGGTGATGTATCCCCTGGGCTTCCGGAAACGGGGCGCCGCGCGCCCCGGACCCGTGTTGTTTCTCGATCCTTTAAACCGTTACCGGGATCATATAATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                           <  1:754955/144‑1 (MQ=255)
                                       aaaaTGGTGAGTGATCCCCTGCGCGTACTGGAACGGGTCGCCGCGCGTACCGGACCCGTGGTGTATCTCGAACCGTTAACCCGCTATCAGGATAATATGACCAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                      <  3:1244607/149‑1 (MQ=255)
                                       aaaaGGGGGAGTGTTTCCCTGCGCGTACGGAAACGGGCCGCCGCGCGTACCGAAACGGTGGTGTACCTCAACCCGTTAACCCGCTATAAGGATAATATAAAAAACCCCCCCCCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                     <  1:530781/150‑1 (MQ=255)
                                        aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                    >  4:132462/1‑150 (MQ=255)
                                            ggTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                >  1:150329/1‑150 (MQ=255)
                                            ggTGAGTGATTCCCTGCGCGTACTGAAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGACCCGTTAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                <  3:1201923/150‑1 (MQ=255)
                                             ggggagtgaTCCCCGGCGCTTACTGGACCGGGCCGCCGCGCTTACCGGACCCGTGTTTTATCTCAACCCTTAAACCCCCTATCAGAATCATAGAACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                <  2:361726/147‑1 (MQ=255)
                                                 agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                              >  3:1364145/1‑148 (MQ=255)
                                                     aTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                       >  1:1273826/1‑151 (MQ=255)
                                                      ttCCCTGCGCTTAGTGAAACAGGCCGCCGCGCGTCCCGGAACGGTGGTTTATCTCGACCCTTAAACCCGCTATAGGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                           <  1:486096/146‑1 (MQ=255)
                                                      ttCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                      >  1:258337/1‑151 (MQ=255)
                                                        cccGGCGCGTACTGAAACGGGTCGCCGCGCGTCCCGGACCGGGGGTGTATCTCGAACCGTTAAACCGCTATCAGATTAATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                     <  4:437034/150‑1 (MQ=255)
                                                          cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                  >  3:52634/1‑151 (MQ=255)
                                                            gcgcGTAGGGGAACGGGCCCCCGCGCGTCCCGGACCGGTGGTTTATCTCGACCCGTTAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTa                                                                                        <  1:718933/143‑1 (MQ=255)
                                                            gcgcGTACTGAAAGAGGTCGCCGCGCGTACCGGACCCGTGTTGTATCCCGAACCTTAAACCCTCTACAAGGACAAAATGACCAACCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                                 <  1:480435/150‑1 (MQ=255)
                                                               cgtccggaacctgttCGCCGCGCTACCCGAACCTGTGGTTTACCTCAACCCTTAACCCCGTAACCAGATCAATTTAATAACCACCCCCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                              <  4:153053/136‑1 (MQ=255)
                                                                  aCTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                          >  1:413346/1‑151 (MQ=255)
                                                                   cggaaaaagtccgccgcGCGTCCCGGAACCGTGGTCTGTCCCGACCCTTTAAGCCGCTATCAGATTAATTTATTCACCCCCCCCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         <  1:1114982/140‑1 (MQ=255)
                                                                    gggAACAGGTCGCCGCGCTAACGGGACCCGTGTTGTCTCCCAACCCGTAAACCCCCAAGAAGTATAATATAATCAACCCCCCCGCAAAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                        <  4:54453/150‑1 (MQ=255)
                                                                     ggAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                       >  2:340574/1‑151 (MQ=255)
                                                                     gaaaCGGGTCGCCGCGCGTACCGGACCGGTGGTGTATCTCGAACCGTAAAACCGCTATCAGGATCATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         <  2:22723/147‑1 (MQ=255)
                                                                     gaaaCAGGTCGCCGCGCGTACCGGAACCGGGGTGTATCTCAACCCGTTAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                          <  2:457031/146‑1 (MQ=255)
                                                                       aaCGGGCCGCCGCGCGTACCGGAACGGTGGTGTATCTCAACCCTTTAACCCGCTACCAGGATCATATGATCACCCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                     <  1:746826/151‑1 (MQ=255)
                                                                         cAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                    >  3:933369/1‑150 (MQ=255)
                                                                             tCGCGGCGCGTACCGAAACCGTGGTGTACCTCGACCCTTTAAACCGCTATAAGGATCATATGACCAACCCCCTCGCCGTTCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                <  3:251053/150‑1 (MQ=255)
                                                                             tCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                               >  3:1399904/1‑151 (MQ=255)
                                                                              cgccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                              >  1:1255444/1‑151 (MQ=255)
                                                                                ccgcGCTTCCGGGACCCGGGGTGTATCTCAAACGTTAAACCCGCTATCAGGATCATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                             <  3:519825/150‑1 (MQ=255)
                                                                                ccgcGCGTCCCGGACCCGTGGTGTATCTCGACCCGTTAACCCGCTATCAGAATAATATGATCAACCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                <  4:13148/147‑1 (MQ=255)
                                                                                ccgcGCGTCCCGGACCCGGGTTGTATCTCAACCCTTAAACCCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                              <  1:308785/149‑1 (MQ=255)
                                                                                ccgcGCGCCCCGGACCGGTGGTGTATCCCGAACCGTAAACCCGCTAAAAGGATAATATGACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                             <  3:127073/150‑1 (MQ=255)
                                                                                ccccgcgTACCGGACCGGGGGTTTTCCCCGACCCTTAAACCCGCTATAAGGATCATGTGAACACCCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                               <  1:731098/145‑1 (MQ=255)
                                                                                 cgcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                            >  2:400840/1‑150 (MQ=255)
                                                                                  gcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGACTATCGGCGATTTCTATCACATGAATATCga                                                            >  3:473798/1‑149 (MQ=255)
                                                                                    gcgTCCCGGACCCGTGGTGTTTCCGGACCCTTTAACCCGCTATCAGAATCTTATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         <  2:313381/150‑1 (MQ=255)
                                                                                      gccccGAAACCGGGGTGTTTCCCAAACCGTTAAACCGATATCAGAATCATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                      <  1:293113/148‑1 (MQ=255)
                                                                                         ccGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                    >  3:1374699/1‑150 (MQ=255)
                                                                                          cGGAAACGTGGGTGTCTCTCAAACCTTAACACCGCTATCAGGATAATATTATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                   <  1:52668/151‑1 (MQ=255)
                                                                                           ggAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                  >  3:976636/1‑150 (MQ=255)
                                                                                            ggaccggtggtgTATCTCGACCCGTAAACCCGCTAACAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                  <  3:473548/144‑1 (MQ=255)
                                                                                            ggaccggtggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATAATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                  <  2:548287/144‑1 (MQ=255)
                                                                                            gAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                >  1:1393086/1‑151 (MQ=255)
                                                                                              accggtggtgTACCTCGACCCTTAAACCCGCTATCAGAATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                <  4:330373/146‑1 (MQ=255)
                                                                                                cGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                            >  2:173219/1‑151 (MQ=255)
                                                                                                   tggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                           >  3:302573/1‑150 (MQ=255)
                                                                                                      gtgtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                        >  4:545175/1‑151 (MQ=255)
                                                                                                       tgtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                       >  1:1259013/1‑151 (MQ=255)
                                                                                                       tgtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                       >  2:181259/1‑151 (MQ=255)
                                                                                                          aTCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                    >  3:1124702/1‑151 (MQ=255)
                                                                                                          aTCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                    >  3:352922/1‑151 (MQ=255)
                                                                                                               cccccctttaccccccTATCAGTTTCATATCATCCCCCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                <  1:933245/136‑1 (MQ=255)
                                                                                                                        aacccGTTATAAGGATCATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACaa                           <  1:660998/144‑1 (MQ=255)
                                                                                                                        aaCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                      >  1:1233110/1‑151 (MQ=255)
                                                                                                                           ccGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                    >  3:186232/1‑151 (MQ=255)
                                                                                                                           ccGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                    >  2:404186/1‑151 (MQ=255)
                                                                                                                            cgttaacAGGATAATTTAACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGa                     <  3:56925/143‑1 (MQ=255)
                                                                                                                            cGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                   >  3:49971/1‑151 (MQ=255)
                                                                                                                            cGCTATAAGGTTCATATCATCACCCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                      <  1:752836/148‑1 (MQ=255)
                                                                                                                            cGCTACCAGAATCAAATGAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                   <  3:580085/151‑1 (MQ=255)
                                                                                                                             gCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                  >  1:1302852/1‑151 (MQ=255)
                                                                                                                               tATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                >  1:87383/1‑151 (MQ=255)
                                                                                                                                aTCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt               >  3:962564/1‑151 (MQ=255)
                                                                                                                                            atGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg     <  3:258357/150‑1 (MQ=255)
                                                                                                                                            atGATCAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 <  4:277501/138‑1 (MQ=255)
                                                                                                                                              gATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  3:39930/1‑151 (MQ=255)
                                                                                                                                               accAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCa   <  3:413411/147‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
CCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  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: