Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5040AC46.5% 101.7 / 25.9 71T291P (ACC→CCC) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base A (30/8);  new base C (0/33);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 5.27e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGG  >  NC_000913/436355‑436651
                                                                                                                                                       |                                                                                                                                                   
tttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTa                                                                                                                                                     <  1:602701/150‑1 (MQ=255)
tttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTa                                                                                                                                                     <  2:665964/150‑1 (MQ=255)
 ttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTAcc                                                                                                                                                   >  1:39898/1‑151 (MQ=255)
     gATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTg                                                                                                                                                 <  2:606415/149‑1 (MQ=255)
      aTGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTAt                                                                                                                                              >  2:303228/1‑151 (MQ=255)
          gcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCggg                                                                                                                                          >  2:217177/1‑151 (MQ=255)
          gcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCggg                                                                                                                                          >  2:557018/1‑151 (MQ=255)
           cTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGc                                                                                                                                         >  2:426551/1‑151 (MQ=255)
                 cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCaaa                                                                                                                                      <  2:69198/148‑1 (MQ=255)
                 cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATg                                                                                                                                    >  1:20571/1‑150 (MQ=255)
                   ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGAcc                                                                                                                                 <  2:114062/151‑1 (MQ=255)
                      tGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGAccgcc                                                                                                                              >  1:477548/1‑151 (MQ=255)
                            aCCGGAACGGGCGCGGCGCGGGGCGCTCTCGGGCGGTGAAGATTACGAGTGGTGTTTCACTGGGCCGGAA‑CTAACCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCT‑TTTCCCTTTAGCGGGCAAATGACCGCCGatat                                                                                                                         <  1:445505/150‑1 (MQ=255)
                                  aCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAggg                                                                                                                  >  2:82197/1‑151 (MQ=255)
                                    aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCt                                                                                                                >  2:525746/1‑151 (MQ=255)
                                     ggCGCTGCGCTGGGCGCTCTCTGGGGGTGAAGATTACGAGTTGTTTTTCACTGTGCCGGAA‑CTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTACCG‑TTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCtt                                                                                                               <  2:75312/151‑1 (MQ=255)
                                     ggCGCTGCGCGGGGCGCTCTCGGGGGGTGAAGATTACGAGTTGTTTTTCACTGTGCCGGAA‑CTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCT‑TTTCCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCt                                                                                                                <  1:24014/150‑1 (MQ=255)
                                         cggcggtgGGCGCTCTCTGGCGGTGAATTTTACGTGTTGTGTTTCTCTGGGCGGGAA‑CTAACCCGGGGCGCGTTGGAGTTGGCTCGGGGAC‑CCTGGGGCTTCCGTTTTCCCGGTCGGGGCCAAAGGACCGCCGATATCGAAGGGCttt                                                                                                              <  1:532158/146‑1 (MQ=255)
                                          cgcgttgGGGGCTTTCGGGGGGTGAAGTTTACGATTTGTTTTTCACTGTGCCGAAA‑CTGACCGGTGGCGCGCTGGATGTGGCTCTGGGACCCCGGGGCGTCCTT‑TTTCCCTGTACGGGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                          <  2:72640/150‑1 (MQ=255)
                                                gggCGTCCTCGGGGGGTGAAGATTACGGGTTGTTTTCCCCTGTGCCGAAT‑CTGACCCGGGGCGCGGTGGTGGTGCCTCCCGGCCACCGGGGCGTCCCG‑TTTCCTTGTGTCGGGAAAATAACCGCCGATATCGAAGGGCttt                                                                                                              <  1:533497/141‑1 (MQ=255)
                                                gggCGCTCTCGGGGGGGGAAGATTACGAGTTGTTTTTCACTGTGCCGGAA‑CTGACCCGGGGCGCGTGGGATGTGGCTCTCGGACCCCTGGGCGTCCCT‑TTTCCTTTTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                    <  1:648604/151‑1 (MQ=255)
                                                 ggCGCTCTCTGGGGGTAAAGTTTACGTGTTGTTTTTCCCGGTGCGGGAA‑CTGACCCGGGGCGCGTGGGATGTGGCTCCGGGACCCTGGGGCGTCCTT‑TTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                                   <  1:938494/151‑1 (MQ=255)
                                                 ggCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTACCT‑TTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                                   <  2:217318/151‑1 (MQ=255)
                                                  gcgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                  >  1:518190/1‑151 (MQ=255)
                                                     ctctctGGCGGTGAAGTTTACGAGTTGTTTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCT‑TTTCCCTGTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                  <  1:909419/148‑1 (MQ=255)
                                                      tctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                              >  2:68057/1‑151 (MQ=255)
                                                        tcggggggTGAAGATTAGGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGGTGGTTGTGGCTCTCGGACACCGGGGCGTCCCT‑TTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                               <  2:428534/142‑1 (MQ=255)
                                                           gggggTGAAGATTACGGGTTGTGTTTCACTGTGCCGGAA‑CTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCT‑TTTCCCTTTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAc                                                                                         <  2:775409/148‑1 (MQ=255)
                                                           ggCGGTGAAGATTACGAGTTGTTTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCCCGGACACCGGGGCGTACCG‑TTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                               <  1:436036/145‑1 (MQ=255)
                                                           ggCGGGGAAGATAACGGGTTGTTTTTCCCTGTGCCGGAT‑CTGAACCGGGGCGCGTTGGATGGGCCTCGCGCACCCCGGGGCTTACCG‑TTTCCCTTTACCGGGAAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                             <  1:462867/147‑1 (MQ=255)
                                                            gCGGTGAAGATTACGGGTTGTGTTTCACTGTGCCGGAA‑CTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTCCCT‑TTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                             <  2:476624/146‑1 (MQ=255)
                                                                       tacgttttGTGTTTCTCTGTGCCGGAA‑CTGACCCGTGGCGCGGGGGATGTGGCTCTCGGACACCTGGGCGTACCT‑TTTACTTGTGTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                  <  1:647959/140‑1 (MQ=255)
                                                                          gttttGTGTTCCCCTGT‑CCGGAACCGGACCCGGGGCGCGCGGGTGGTGGCTCCCGGACCCCTGGGCGTACCT‑TTTCCTTTTAGCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa                                                                            <  1:938911/146‑1 (MQ=255)
                                                                            gttgttTTTCACTGTGCCGGAA‑CGAAACCGGGGGGCGCTGGATGGGGCTCTCGGACCCCGGGGCGTACTT‑TTTACCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                             <  1:745309/146‑1 (MQ=255)
                                                                             ttGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                        >  1:725575/1‑150 (MQ=255)
                                                                              tgtgTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                       >  1:700852/1‑150 (MQ=255)
                                                                                ttttccccTGTGCCGAAA‑CTAACCCGGGGCGCGGTGGATGTGGCTCCCGGACCCGTGGGGGTCCCT‑TTTCCCTGTATCGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                         <  1:813190/139‑1 (MQ=255)
                                                                                  tttCCCTGTGCCGGAA‑CTGAACGGTGGCGCGCTGGATGTGGCTCTCGGACACGGGGGCGTCCCG‑TTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                   <  2:142055/150‑1 (MQ=255)
                                                                                  tttCACTGTGCCGGAA‑CTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCT‑TTTCCTTTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                   <  2:413192/150‑1 (MQ=255)
                                                                                  ttcccttgtgCCGAAA‑CTGACCGGTGGCGCGCGGGATGGGGCTCGCGGACCCGTGGGCGTCCCT‑TTTCCTTGAATCGGGCAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                   <  1:248562/144‑1 (MQ=255)
                                                                                     cacggggCCGGAA‑CGGAACGGGGGCGCGGGGGATGTGGCTCCGGGACACCGGGGCGTCCCT‑TTTCCCTGTATGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                    <  2:359280/140‑1 (MQ=255)
                                                                                      aCTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                              >  1:425397/1‑151 (MQ=255)
                                                                                           gCCGGAA‑CTGAACCGGGGCGCGCGGGATGTGGCTCTCGGCCCCCGGGGCGTACCT‑TTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCa                                                          <  2:372408/150‑1 (MQ=255)
                                                                                            cggaaa‑CGGACCCGTGGCGCGTGGTATGGGCCTCTCGGCCCCCGGGGCGTCCCT‑TTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAca                                                           <  2:372390/144‑2 (MQ=255)
                                                                                            ccGGAA‑CTGACCCGCGGGGCTCTGGATGTGCCTCTCGGACACCGGGGCGTACTT‑TTTCCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                         <  1:793631/150‑1 (MQ=255)
                                                                                                aa‑CTGACCCGGGGCGCGTTGGATGTGGCCCGCGGACACCTGGGCTTTCCT‑TTTACCTGTTTCGGGAAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                    <  1:115147/151‑1 (MQ=255)
                                                                                                aa‑CTAACCCGTGGCGCGTTGGTGGGGCCTCTCGGACCCCGGGGCGTCCCT‑TTTCCTTTTATCGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                    <  1:909806/151‑1 (MQ=255)
                                                                                                   cTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc                                                   >  2:897582/1‑150 (MQ=255)
                                                                                                     gAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                >  1:161737/1‑151 (MQ=255)
                                                                                                      agctggggCGCGTGGGTTGGGGCTCCCGGCCCCTTGGGCGTCCCG‑TTTCCTTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                 <  2:937551/143‑1 (MQ=255)
                                                                                                      aaCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                               >  1:863315/1‑151 (MQ=255)
                                                                                                      aaCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCT‑TTTCCTTGTTGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCa                                                  <  2:767644/148‑1 (MQ=255)
                                                                                                               gcgcTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                          >  1:412068/1‑147 (MQ=255)
                                                                                                                 ggtGGATGTGGCTCTCGGCCACGGGGGCGTCCCT‑TTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAg                                        <  2:336129/145‑1 (MQ=255)
                                                                                                                      aTGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgcg                                >  2:430222/1‑150 (MQ=255)
                                                                                                                          ggCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa                             >  2:668097/1‑149 (MQ=255)
                                                                                                                          ggCTCTCGGACACCGGGGCGTCCCT‑TTTCCTTTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaa                              <  2:781084/148‑1 (MQ=255)
                                                                                                                                cGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                     >  2:633000/1‑151 (MQ=255)
                                                                                                                                 ggccacgggggCGTCCCG‑TTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                     <  2:509139/142‑1 (MQ=255)
                                                                                                                                 ggACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                    >  2:3223/1‑151 (MQ=255)
                                                                                                                                     acCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                   >  2:552764/1‑148 (MQ=255)
                                                                                                                                     acCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                   >  2:637881/1‑148 (MQ=255)
                                                                                                                                      ccggggCGTCCCT‑TTTCCTTGTTTCGGCCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAg                 <  2:379071/146‑1 (MQ=255)
                                                                                                                                      ccTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                   >  1:545231/1‑147 (MQ=255)
                                                                                                                                         gggCGTCCCT‑TTTCCTTGTATGGGAAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGaa                  <  2:524983/145‑1 (MQ=255)
                                                                                                                                         gggCGTACCT‑TTTCCCTTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGAt               <  2:518242/148‑1 (MQ=255)
                                                                                                                                          ggCGTCCCT‑TTTCCCTTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt           <  1:426509/151‑1 (MQ=255)
                                                                                                                                          ggCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt           >  2:185906/1‑151 (MQ=255)
                                                                                                                                               tcccttttCCCTGTAGGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTa          <  1:856256/143‑1 (MQ=255)
                                                                                                                                               aCCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc      >  1:356342/1‑151 (MQ=255)
                                                                                                                                                  g‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg   >  1:611423/1‑151 (MQ=255)
                                                                                                                                                  g‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg   >  1:937146/1‑151 (MQ=255)
                                                                                                                                                    tttACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTgg  >  1:25114/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
TTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGG  >  NC_000913/436355‑436651

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