Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5160CA43.9% 78.9 / 27.8 66Q295K (CAA→AAA) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (26/11);  new base A (0/29);  total (26/40)
Fisher's exact test for biased strand distribution p-value = 5.65e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTT  >  NC_000913/436366‑436660
                                                                                                                                                        |                                                                                                                                                
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)
                                                                                                                                                  aTCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt  >  2:412075/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                
CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAA‑CTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCG‑TTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTT  >  NC_000913/436366‑436660

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