Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5160CA42.0% 109.8 / 37.6 81Q295K (CAA→AAA) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (38/9);  new base A (0/34);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 7.71e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436367‑436666
                                                                                                                                                       |                                                                                                                                                      
tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                        >  2:194122/1‑150 (MQ=255)
 ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                        <  2:757273/149‑1 (MQ=255)
 ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                        <  1:259935/149‑1 (MQ=255)
    tcGCCATGTTGAACCGGAACCGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGCGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATg                                                                                                                                                   >  1:571802/1‑151 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATg                                                                                                                                                   >  1:1194734/1‑150 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATg                                                                                                                                                   <  1:33177/150‑1 (MQ=255)
      gCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAc                                                                                                                                                 >  2:153585/1‑151 (MQ=255)
       ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAcc                                                                                                                                                <  1:29975/151‑1 (MQ=255)
         aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgc                                                                                                                                              >  2:74954/1‑151 (MQ=255)
         aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAAATGCACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgc                                                                                                                                              >  1:1143660/1‑151 (MQ=255)
           gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgccg                                                                                                                                            <  1:1142017/151‑1 (MQ=255)
             tGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGTCGCGCTGGATTTTGCTCTCTTCCCCCTTGTCTTACCGTTTACCTGTATCGGGCAAATGAccgccgct                                                                                                                                          >  2:511343/1‑149 (MQ=255)
              gAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGat                                                                                                                                          >  2:395402/1‑150 (MQ=255)
                aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATc                                                                                                                                       >  2:1332465/1‑151 (MQ=255)
                   ggATCAGGCGCTGGCCGGGCCTCTCTCGGGCGGTGAAGATAAGGAGTTGTTTTCCACGGTGCTGGAACTGAACCGGGGCGCTCTGGTTGTGGCTCTCGGACACTGGGACGTACGGTTTCCCGTTACGGGGCAAATGACCGCCGATATCGaa                                                                                                                                    <  1:865177/151‑1 (MQ=255)
                     aaCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAgg                                                                                                                                  >  2:86601/1‑151 (MQ=255)
                       cAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGAACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGata                                                                                                                                         >  1:351888/1‑142 (MQ=255)
                       cAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGCACCGTTTACCTGTATCGGGCAAATGGCCGCCGATATCGAAGGGc                                                                                                                                >  1:234467/1‑151 (MQ=255)
                         ggcgcggcgcGGGCCTCTCTCGGGGGTTGAAGATTAGGATTTTTTTTTCACTGTGCCGAAACTGACCCGTGGCGCGCGGGATGTGCCTCGCGGACCCCGGGGGGTCCTTTTTCCCTTTATCGGGAAAATGACCGCCGATATCGAAGGGCtt                                                                                                                              <  1:959486/151‑1 (MQ=255)
                            gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                           >  1:1116263/1‑151 (MQ=255)
                            gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                           >  1:851436/1‑151 (MQ=255)
                            gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                           >  2:1102771/1‑151 (MQ=255)
                             ctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                          >  1:1138776/1‑151 (MQ=255)
                               gggctgGGCGCGCGCGGGCGGGGAAGTTTACGGGTTGTTTTCCCCTGTGCCGGACCTGACCCGGGGCGCTGTGGATGTGGCTCGCGGCCACGGGGGGGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                         <  1:730932/148‑1 (MQ=255)
                               gcggggggCGCTTTCGGGGGGTGAATATTACGAGTTTTTTTTCCCTGTGCCGAAACTGACCCGGGGCGCGGTGGATGTGGCTCTCGCCCACGGGGGCGTACTTTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTg                                                                                                                            <  1:877624/142‑1 (MQ=255)
                                 gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                      >  1:766931/1‑151 (MQ=255)
                                   ggggCGCTCTCGGGCGGTGAAGTTTACGGTTTTTTTTTCCCTGGGCGGGAACTGACCCGTGGCGCGGGGGGTGTGGCCCCCGCCCCCGTGGGCGTACTTTTTCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                                    <  1:793341/150‑1 (MQ=255)
                                    gggCGCTCTCGGGGGGTGAAGATTACGTGTTGTGTTTCACGTTGCCGGAACTGAACCGTGGCGCGTTGGATGTGGCTCTGGGCCACTGGGGCGTCCTTTTTACTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                      <  1:655518/148‑1 (MQ=255)
                                    gggCGCTCTCGGGCGGTGAAGATTACGGGTTTTGTTCCCCTGTGCCGGAACTGACCGGTGGCGCGTTGGATGTGCCTCTCGGACCCTTGGGCGTCCCTTTTCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                      <  2:78459/148‑1 (MQ=255)
                                     gtcgcCCTCTGGCGGGGAAAATTACGTGTTGTTTTTCTCTGAGCCGGAACGGACCCGTGGCGCTTTGGTTGTGGCTCGCGGACACGTGGGCGTACCTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                                    <  2:762972/147‑1 (MQ=255)
                                     gtccctttctGGGGATAAAGATTAGGGGTTGTGTTTCAGTTTGCCGAAACGGAGCTGGGCGGCCCTGGTTTGGGTTCAGGCACCCTGGGTGCTACATTTTACTTTAATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                                   <  2:960873/143‑1 (MQ=255)
                                       cgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                                 >  2:1315640/1‑150 (MQ=255)
                                          tctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                             >  1:612835/1‑151 (MQ=255)
                                               gggggTGAATATTACGAGTGGTGTTTCCCTGGGCGGGAACTGACCCGTGGCGCGTTGGATGTGGCTCCCGGACACGGGGGCGTCCGTTTTACTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                             <  1:252524/143‑1 (MQ=255)
                                                  ggTGAAGATTACGTGTTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGCTGGTTGTGGCTCGCGGACACCGGGGCGTCCCGTTTCCCTGTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                     <  2:13208/151‑1 (MQ=255)
                                                      aaGTTTACGAGTTGTGTTTCCTTGTGCCGGAACTGACCGGTGGCGCGGTGGATGGGGCTCGGGAACCCTGGGGCGTCCCTTTTTCCTGTTGCGGGAAAATGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                                 <  2:463662/151‑1 (MQ=255)
                                                      aaGATTACGTGTTGTTTTTCCCTGTGCGGGAACGGACCGGTGGCGGGTTGGATGGGCCTCTCGACCCCCGGGGCGTCCCTTTTCCTTGTATCGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                                 <  1:967777/151‑1 (MQ=255)
                                                      aaGATTACGAGTTGTGTTTCCCTGTGCGGGAACTGACCCGTGGCGCGTTGGATGGGGCTCTCGGACCCCGGGGCGTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                                   <  2:1320135/149‑1 (MQ=255)
                                                          ttACGAGTTGTGTTTCACTGTGCGGAAACTGACCGGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCCTGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                              <  2:1349301/150‑1 (MQ=255)
                                                          ttACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                              >  2:630824/1‑150 (MQ=255)
                                                               aGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                        >  1:1153407/1‑151 (MQ=255)
                                                                    tttttCACTGTGCCGAAACGGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCTTACCTTTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                                                    <  2:850971/148‑1 (MQ=255)
                                                                      tttCCCTGTGCCGAAATTAACCCGTGGCGCGCTGGATGTGTCTCTCGGACCCGTGGGCGTCCCTTTTCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                                  <  2:454804/150‑1 (MQ=255)
                                                                      tttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                                  >  1:951/1‑150 (MQ=255)
                                                                         cccTGTGCGGGAACGAACCCGTGGCGCGCTGGTTGTGGCTCCGGGACCCCGGGGCGTACTTTTTCCCTGTATCGGACAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                              <  2:218387/149‑1 (MQ=255)
                                                                            tgtgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                           >  1:1059514/1‑151 (MQ=255)
                                                                            tgtgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                           >  1:635409/1‑151 (MQ=255)
                                                                               gCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                        >  1:1296826/1‑151 (MQ=255)
                                                                                 cGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                      >  2:111230/1‑151 (MQ=255)
                                                                                  ggAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                     >  2:273531/1‑151 (MQ=255)
                                                                                    aaCTGAACCGTGGCGCGTTGGATGTGGCTCTCGGACCCCGGGGCGTCCCGTTTCCCTTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                     <  2:1060557/149‑1 (MQ=255)
                                                                                    aaCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                   >  1:551331/1‑151 (MQ=255)
                                                                                     atagaaggggggggcgcTGGGTGTGGCTCTCGGCCCCCGGGGGGTCCTTTTTCCCTGTCTGGGGCAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                           <  2:1291800/129‑1 (MQ=255)
                                                                                        gacccGGGGGGCGCGGGTTGTGGCTCTCGGACCCCGGGGTGTCCCGTTTCCCTGTAGGGGGCAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCa                                                                 <  1:91639/146‑1 (MQ=255)
                                                                                          aCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                             >  2:11123/1‑151 (MQ=255)
                                                                                           cgggggCGCGCTGGATGGGGCTCTCGGCCCCCGGGGCGTACCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                               <  2:640173/144‑1 (MQ=255)
                                                                                           cggTGGCGCGTTGGATGTGGCTCTCGGACCCCTGGGGGTCCCTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa                                                            <  2:805603/149‑1 (MQ=255)
                                                                                             ggggCGCGGTGGTTGGGCCTCCCGCACCCGGGGGGCTACCTTTTTCCTGTAGGGGCCAAATGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                           <  2:743566/148‑1 (MQ=255)
                                                                                                   cgcTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                                    >  2:1187524/1‑151 (MQ=255)
                                                                                                      gggATGGGCCCCTCGGACACCTGGGCTTCCGTTTTCCCTGATCCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                                 <  1:38598/150‑1 (MQ=255)
                                                                                                             ggCTCTCGGACCCCGGGGCGTCCCTTTTCCTTGTATGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAg                                          <  2:1306689/151‑1 (MQ=255)
                                                                                                                ctatcggccccTGGGGGTTTCGTTTTCCTTTTATGGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGGTGTcg                                                 <  2:357133/133‑1 (MQ=255)
                                                                                                                 ctcGGATCCCGGGGCGTCCCTTTTCCCTGATTCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                       <  1:1331595/150‑1 (MQ=255)
                                                                                                                 cggggcccccggggCGTCCCTTTTCCCTGTAGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGt                                         <  1:728238/137‑1 (MQ=255)
                                                                                                                 cccatcccctggggCGTACCTTTTTCCTTGATGGACCACATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAGTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                        <  1:62581/138‑1 (MQ=255)
                                                                                                                    ggccacTTGGGGGTCCTGTTTCCCTATATGGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                       <  2:988559/144‑1 (MQ=255)
                                                                                                                    ggACCCCGGGGCGTCCCTTTTCCTTTTAGCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGcc                                     <  2:859413/149‑1 (MQ=255)
                                                                                                                      acacCTGGGCGTACCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc                                      <  2:340184/146‑1 (MQ=255)
                                                                                                                       cacCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                                  >  1:915753/1‑149 (MQ=255)
                                                                                                                        cctggggCGTCCTTTTTCCCGGTATGGGGAAAGGGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                                  <  1:1272323/144‑1 (MQ=255)
                                                                                                                                gTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAAt                       >  1:42811/1‑151 (MQ=255)
                                                                                                                                gTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGACTGAAGATGAGTAAt                       >  1:584343/1‑151 (MQ=255)
                                                                                                                                  tACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCCTGACGGCGAACCTGTTACATTAGACTGGAAAGGCTATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc                      >  1:233973/1‑151 (MQ=255)
                                                                                                                                   aCCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                     >  2:1371528/1‑151 (MQ=255)
                                                                                                                                   aCCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATAACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                     >  1:327913/1‑151 (MQ=255)
                                                                                                                                   aCCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAACTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                     >  1:1019481/1‑151 (MQ=255)
                                                                                                                                       tttCCCTTTGTCGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc                      <  2:933955/146‑1 (MQ=255)
                                                                                                                                        tttccttGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg                  <  2:1323599/144‑1 (MQ=255)
                                                                                                                                        ttcccTGTATCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTgg                 <  2:1302633/147‑1 (MQ=255)
                                                                                                                                           ccTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATc             >  1:910634/1‑151 (MQ=255)
                                                                                                                                                    gggCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCt     >  2:1262614/1‑151 (MQ=255)
                                                                                                                                                    gggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTAc          <  2:1046917/146‑1 (MQ=255)
                                                                                                                                                    gggAAAATGACCGCCGATATCGAAGGGCTTTGGTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt        <  2:118724/148‑1 (MQ=255)
                                                                                                                                                    gggAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt         <  2:59770/147‑1 (MQ=255)
                                                                                                                                                       cAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATCAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTc  >  1:306538/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436367‑436666

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