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

TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436367‑436666
                                                                                                                                                       |                                                                                                                                                      
tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGc                                                                                                                                                        >  2:194122/1‑150 (MQ=255)
 ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGc                                                                                                                                                        <  2:757273/149‑1 (MQ=255)
 ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGc                                                                                                                                                        <  1:259935/149‑1 (MQ=255)
    tcGCCATGTTGAACCGGAACCGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGCGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATg                                                                                                                                                   >  1:571802/1‑151 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATg                                                                                                                                                   >  1:1194734/1‑150 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATg                                                                                                                                                   <  1:33177/150‑1 (MQ=255)
      gCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAc                                                                                                                                                 >  2:153585/1‑151 (MQ=255)
       ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAcc                                                                                                                                                <  1:29975/151‑1 (MQ=255)
         aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccgc                                                                                                                                              >  2:74954/1‑151 (MQ=255)
         aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAAATGCACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccgc                                                                                                                                              >  1:1143660/1‑151 (MQ=255)
           gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccgccg                                                                                                                                            <  1:1142017/151‑1 (MQ=255)
             tGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGTCGCGCTGGATTTTGCTCTCTTCCCCCTTGTCT‑TACCGTTTACC‑TGTATCGGGCAAATGAccgccgct                                                                                                                                          >  2:511343/1‑149 (MQ=255)
              gAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGat                                                                                                                                          >  2:395402/1‑150 (MQ=255)
                aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATc                                                                                                                                       >  2:1332465/1‑151 (MQ=255)
                   ggATCAGGCGCTGGCCGGGCCTCTCTCGGGCGGTGAAGATAAGGAGTTGTTTTCCACGGTGCTGGAACTGAACCGGGGCGCTCTGGTTGTGGCTCTCGGACACTGGGACG‑TACGGTTTCCC‑GTTACGGGGCAAATGACCGCCGATATCGaa                                                                                                                                    <  1:865177/151‑1 (MQ=255)
                     aaCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAgg                                                                                                                                  >  2:86601/1‑151 (MQ=255)
                       cAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGAACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGata                                                                                                                                         >  1:351888/1‑142 (MQ=255)
                       cAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑CACCGTTTACC‑TGTATCGGGCAAATGGCCGCCGATATCGAAGGGc                                                                                                                                >  1:234467/1‑151 (MQ=255)
                         ggcgcggcgcGGGCCTCTCTCGGGGGTTGAAGATTAGGATTTTTTTTTCACTGTGCCGAAACTGACCCGTGGCGCGCGGGATGTGCCTCGCGGACCCCGGGGGG‑TCCTTTTTCCC‑TTTATCGGGAAAATGACCGCCGATATCGAAGGGCtt                                                                                                                              <  1:959486/151‑1 (MQ=255)
                            gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                           >  1:1116263/1‑151 (MQ=255)
                            gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                           >  1:851436/1‑151 (MQ=255)
                            gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                           >  2:1102771/1‑151 (MQ=255)
                             ctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                          >  1:1138776/1‑151 (MQ=255)
                               gggctgGGCGCGCGCGGGCGGGGAAGTTTACGGGTTGTTTTCCCCTGTGCCGGACCTGACCCGGGGCGCTGTGGATGTGGCTCGCGGCCACGGGGGGG‑TCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                         <  1:730932/148‑1 (MQ=255)
                               gcggggggCGCTTTCGGGGGGTGAATATTACGAGTTTTTTTTCCCTGTGCCGAAACTGACCCGGGGCGCGGTGGATGTGGCTCTCGCCCACGGGGGCG‑TACTTTTTCCT‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTg                                                                                                                            <  1:877624/142‑1 (MQ=255)
                                 gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                      >  1:766931/1‑151 (MQ=255)
                                   ggggCGCTCTCGGGCGGTGAAGTTTACGGTTTTTTTTTCCCTGGGCGGGAACTGACCCGTGGCGCGGGGGGTGTGGCCCCCGCCCCCGTGGGCG‑TACTTTTTCCC‑TGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                                    <  1:793341/150‑1 (MQ=255)
                                    gggCGCTCTCGGGGGGTGAAGATTACGTGTTGTGTTTCACGTTGCCGGAACTGAACCGTGGCGCGTTGGATGTGGCTCTGGGCCACTGGGGCG‑TCCTTTTTACT‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                      <  1:655518/148‑1 (MQ=255)
                                    gggCGCTCTCGGGCGGTGAAGATTACGGGTTTTGTTCCCCTGTGCCGGAACTGACCGGTGGCGCGTTGGATGTGCCTCTCGGACCCTTGGGCG‑TCCCTTTTCCC‑TGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                      <  2:78459/148‑1 (MQ=255)
                                     gtcgcCCTCTGGCGGGGAAAATTACGTGTTGTTTTTCTCTGAGCCGGAACGGACCCGTGGCGCTTTGGTTGTGGCTCGCGGACACGTGGGCG‑TACCTTTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                                    <  2:762972/147‑1 (MQ=255)
                                     gtccctttctGGGGATAAAGATTAGGGGTTGTGTTTCAGTTTGCCGAAACGGAGCTGGGCGGCCCTGGTTTGGGTTCAGGCAC‑CCTGGGTGCTACATTTTACT‑TTAATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                                   <  2:960873/143‑1 (MQ=255)
                                       cgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                                 >  2:1315640/1‑150 (MQ=255)
                                          tctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                             >  1:612835/1‑151 (MQ=255)
                                               gggggTGAATATTACGAGTGGTGTTTCCCTGGGCGGGAACTGACCCGTGGCGCGTTGGATGTGGCTCCCGGACACGGGGGCG‑TCCGTTTTACT‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                             <  1:252524/143‑1 (MQ=255)
                                                  ggTGAAGATTACGTGTTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGCTGGTTGTGGCTCGCGGACACCGGGGCG‑TCCCGTTTCCC‑TGTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                     <  2:13208/151‑1 (MQ=255)
                                                      aaGTTTACGAGTTGTGTTTCCTTGTGCCGGAACTGACCGGTGGCGCGGTGGATGGGGCTCGGGAACCCTGGGGCG‑TCCCTTTTTCC‑TGTTGCGGGAAAATGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                                 <  2:463662/151‑1 (MQ=255)
                                                      aaGATTACGTGTTGTTTTTCCCTGTGCGGGAACGGACCGGTGGCGGGTTGGATGGGCCTCTCGACCCCCGGGGCG‑TCCCTTTTCCT‑TGTATCGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                                 <  1:967777/151‑1 (MQ=255)
                                                      aaGATTACGAGTTGTGTTTCCCTGTGCGGGAACTGACCCGTGGCGCGTTGGATGGGGCTCTCGGACCCCGGGGCG‑TCCCGTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                                   <  2:1320135/149‑1 (MQ=255)
                                                          ttACGAGTTGTGTTTCACTGTGCGGAAACTGACCGGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCG‑TCCCTTTTCCC‑TGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                              <  2:1349301/150‑1 (MQ=255)
                                                          ttACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                              >  2:630824/1‑150 (MQ=255)
                                                               aGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                        >  1:1153407/1‑151 (MQ=255)
                                                                    tttttCACTGTGCCGAAACGGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCT‑TACCTTTTACC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                                                    <  2:850971/148‑1 (MQ=255)
                                                                      tttCCCTGTGCCGAAATTAACCCGTGGCGCGCTGGATGTGTCTCTCGGACCCGTGGGCG‑TCCCTTTTCCC‑TGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                                  <  2:454804/150‑1 (MQ=255)
                                                                      tttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                                  >  1:951/1‑150 (MQ=255)
                                                                         cccTGTGCGGGAACGAACCCGTGGCGCGCTGGTTGTGGCTCCGGGACCCCGGGGCG‑TACTTTTTCCC‑TGTATCGGACAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                              <  2:218387/149‑1 (MQ=255)
                                                                            tgtgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                           >  1:1059514/1‑151 (MQ=255)
                                                                            tgtgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                           >  1:635409/1‑151 (MQ=255)
                                                                               gCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                        >  1:1296826/1‑151 (MQ=255)
                                                                                 cGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                      >  2:111230/1‑151 (MQ=255)
                                                                                  ggAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                     >  2:273531/1‑151 (MQ=255)
                                                                                    aaCTGAACCGTGGCGCGTTGGATGTGGCTCTCGGACCCCGGGGCG‑TCCCGTTTCCC‑TTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                     <  2:1060557/149‑1 (MQ=255)
                                                                                    aaCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                   >  1:551331/1‑151 (MQ=255)
                                                                                     atagaaggggggggcgcTGGGTGTGGCTCTCGGCCCCCGGGGGG‑TCCTTTTTCCC‑TGTCTGGGGCAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                           <  2:1291800/129‑1 (MQ=255)
                                                                                        gacccGGGGGGCGCGGGTTGTGGCTCTCGGACCCCGGGGTG‑TCCCGTTTCCC‑TGTAGGGGGCAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCa                                                                 <  1:91639/146‑1 (MQ=255)
                                                                                          aCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                             >  2:11123/1‑151 (MQ=255)
                                                                                           cgggggCGCGCTGGATGGGGCTCTCGGCCCCCGGGGCG‑TACCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                               <  2:640173/144‑1 (MQ=255)
                                                                                           cggTGGCGCGTTGGATGTGGCTCTCGGACCCCTGGGGG‑TCCCTTTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa                                                            <  2:805603/149‑1 (MQ=255)
                                                                                             ggggCGCGGTGGTTGGGCCTCCCGCACCCGGGGGGC‑TACCTTTTTCC‑TGTAGGGGCCAAATGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                           <  2:743566/148‑1 (MQ=255)
                                                                                                   cgcTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                                    >  2:1187524/1‑151 (MQ=255)
                                                                                                      gggATGGGCCCCTCGGACACCTGGGCT‑TCCGTTTTCCC‑TGATCCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                                 <  1:38598/150‑1 (MQ=255)
                                                                                                             ggCTCTCGGACCCCGGGGCG‑TCCCTTTTCCT‑TGTATGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAg                                          <  2:1306689/151‑1 (MQ=255)
                                                                                                                ctatcggccccTGGGGG‑TTTCGTTTTCCTTTTATGGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGGTGTcg                                                 <  2:357133/133‑1 (MQ=255)
                                                                                                                 ctcGGATCCCGGGGCG‑TCCCTTTTCCC‑TGATTCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                       <  1:1331595/150‑1 (MQ=255)
                                                                                                                 cggggcccccggggCG‑TCCCTTTTCCC‑TGTAGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGt                                         <  1:728238/137‑1 (MQ=255)
                                                                                                                 cccatcccctggggCG‑TACCTTTTTCC‑TTGATGGACCACATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAGTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                        <  1:62581/138‑1 (MQ=255)
                                                                                                                    ggccacTTGGGGG‑TCCTGTTTCCC‑TATATGGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                       <  2:988559/144‑1 (MQ=255)
                                                                                                                    ggACCCCGGGGCG‑TCCCTTTTCCT‑TTTAGCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGcc                                     <  2:859413/149‑1 (MQ=255)
                                                                                                                      acacCTGGGCG‑TACCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc                                      <  2:340184/146‑1 (MQ=255)
                                                                                                                       cacCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                                  >  1:915753/1‑149 (MQ=255)
                                                                                                                        cctggggCG‑TCCTTTTTCCC‑GGTATGGGGAAAGGGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                                  <  1:1272323/144‑1 (MQ=255)
                                                                                                                                g‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAAt                       >  1:42811/1‑151 (MQ=255)
                                                                                                                                g‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGACTGAAGATGAGTAAt                       >  1:584343/1‑151 (MQ=255)
                                                                                                                                  tACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCCTGACGGCGAACCTGTTACATTAGACTGGAAAGGCTATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc                      >  1:233973/1‑151 (MQ=255)
                                                                                                                                   aCCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                     >  2:1371528/1‑151 (MQ=255)
                                                                                                                                   aCCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATAACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                     >  1:327913/1‑151 (MQ=255)
                                                                                                                                   aCCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAACTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                     >  1:1019481/1‑151 (MQ=255)
                                                                                                                                       tttCCC‑TTTGTCGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc                      <  2:933955/146‑1 (MQ=255)
                                                                                                                                        tttccttGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg                  <  2:1323599/144‑1 (MQ=255)
                                                                                                                                        ttccc‑TGTATCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTgg                 <  2:1302633/147‑1 (MQ=255)
                                                                                                                                           cc‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATc             >  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)
                                                                                                                                                       |                                                                                                                                                      
TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  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: ‑