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

CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTT  >  NC_000913/436366‑436660
                                                                                                                                                        |                                                                                                                                                
cTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                  >  2:426551/1‑151 (MQ=255)
      cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCaaa                                                                                                                                               <  2:69198/148‑1 (MQ=255)
      cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATg                                                                                                                                             >  1:20571/1‑150 (MQ=255)
        ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAcc                                                                                                                                          <  2:114062/151‑1 (MQ=255)
           tGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgcc                                                                                                                                       >  1:477548/1‑151 (MQ=255)
                 aCCGGAACGGGCGCGGCGCGGGGCGCTCTCGGGCGGTGAAGATTACGAGTGGTGTTTCACTGGGCCGGAACTAACCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCTTTTCCCTTTAGCGGGCAAATGACCGCCGatat                                                                                                                                  <  1:445505/150‑1 (MQ=255)
                       aCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAggg                                                                                                                           >  2:82197/1‑151 (MQ=255)
                         aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCt                                                                                                                         >  2:525746/1‑151 (MQ=255)
                          ggCGCTGCGCTGGGCGCTCTCTGGGGGTGAAGATTACGAGTTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTACCGTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCtt                                                                                                                        <  2:75312/151‑1 (MQ=255)
                          ggCGCTGCGCGGGGCGCTCTCGGGGGGTGAAGATTACGAGTTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCt                                                                                                                         <  1:24014/150‑1 (MQ=255)
                              cggcggtgGGCGCTCTCTGGCGGTGAATTTTACGTGTTGTGTTTCTCTGGGCGGGAACTAACCCGGGGCGCGTTGGAGTTGGCTCGGGGACCCTGGGGCTTCCGTTTTCCCGGTCGGGGCCAAAGGACCGCCGATATCGAAGGGCttt                                                                                                                       <  1:532158/146‑1 (MQ=255)
                               cgcgttgGGGGCTTTCGGGGGGTGAAGTTTACGATTTGTTTTTCACTGTGCCGAAACTGACCGGTGGCGCGCTGGATGTGGCTCTGGGACCCCGGGGCGTCCTTTTTCCCTGTACGGGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                   <  2:72640/150‑1 (MQ=255)
                                     gggCGTCCTCGGGGGGTGAAGATTACGGGTTGTTTTCCCCTGTGCCGAATCTGACCCGGGGCGCGGTGGTGGTGCCTCCCGGCCACCGGGGCGTCCCGTTTCCTTGTGTCGGGAAAATAACCGCCGATATCGAAGGGCttt                                                                                                                       <  1:533497/141‑1 (MQ=255)
                                     gggCGCTCTCGGGGGGGGAAGATTACGAGTTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGTGGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCTTTTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                             <  1:648604/151‑1 (MQ=255)
                                      ggCGCTCTCTGGGGGTAAAGTTTACGTGTTGTTTTTCCCGGTGCGGGAACTGACCCGGGGCGCGTGGGATGTGGCTCCGGGACCCTGGGGCGTCCTTTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                                            <  1:938494/151‑1 (MQ=255)
                                      ggCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTACCTTTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                                            <  2:217318/151‑1 (MQ=255)
                                       gcgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                           >  1:518190/1‑151 (MQ=255)
                                          ctctctGGCGGTGAAGTTTACGAGTTGTTTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCTTTTCCCTGTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                           <  1:909419/148‑1 (MQ=255)
                                           tctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                       >  2:68057/1‑151 (MQ=255)
                                             tcggggggTGAAGATTAGGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGGTGGTTGTGGCTCTCGGACACCGGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                        <  2:428534/142‑1 (MQ=255)
                                                gggggTGAAGATTACGGGTTGTGTTTCACTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTTTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAc                                                                                                  <  2:775409/148‑1 (MQ=255)
                                                ggCGGTGAAGATTACGAGTTGTTTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCCCGGACACCGGGGCGTACCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                        <  1:436036/145‑1 (MQ=255)
                                                ggCGGGGAAGATAACGGGTTGTTTTTCCCTGTGCCGGATCTGAACCGGGGCGCGTTGGATGGGCCTCGCGCACCCCGGGGCTTACCGTTTCCCTTTACCGGGAAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                      <  1:462867/147‑1 (MQ=255)
                                                 gCGGTGAAGATTACGGGTTGTGTTTCACTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                      <  2:476624/146‑1 (MQ=255)
                                                            tacgttttGTGTTTCTCTGTGCCGGAACTGACCCGTGGCGCGGGGGATGTGGCTCTCGGACACCTGGGCGTACCTTTTACTTGTGTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                           <  1:647959/140‑1 (MQ=255)
                                                               gttttGTGTTCCCCTGTCCGGAACCGGACCCGGGGCGCGCGGGTGGTGGCTCCCGGACCCCTGGGCGTACCTTTTCCTTTTAGCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa                                                                                     <  1:938911/146‑1 (MQ=255)
                                                                 gttgttTTTCACTGTGCCGGAACGAAACCGGGGGGCGCTGGATGGGGCTCTCGGACCCCGGGGCGTACTTTTTACCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                                      <  1:745309/146‑1 (MQ=255)
                                                                  ttGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                 >  1:725575/1‑150 (MQ=255)
                                                                   tgtgTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                >  1:700852/1‑150 (MQ=255)
                                                                     ttttccccTGTGCCGAAACTAACCCGGGGCGCGGTGGATGTGGCTCCCGGACCCGTGGGGGTCCCTTTTCCCTGTATCGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                  <  1:813190/139‑1 (MQ=255)
                                                                       tttCCCTGTGCCGGAACTGAACGGTGGCGCGCTGGATGTGGCTCTCGGACACGGGGGCGTCCCGTTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                            <  2:142055/150‑1 (MQ=255)
                                                                       tttCACTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCTTTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                            <  2:413192/150‑1 (MQ=255)
                                                                       ttcccttgtgCCGAAACTGACCGGTGGCGCGCGGGATGGGGCTCGCGGACCCGTGGGCGTCCCTTTTCCTTGAATCGGGCAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                            <  1:248562/144‑1 (MQ=255)
                                                                          cacggggCCGGAACGGAACGGGGGCGCGGGGGATGTGGCTCCGGGACACCGGGGCGTCCCTTTTCCCTGTATGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                             <  2:359280/140‑1 (MQ=255)
                                                                           aCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                       >  1:425397/1‑151 (MQ=255)
                                                                                gCCGGAACTGAACCGGGGCGCGCGGGATGTGGCTCTCGGCCCCCGGGGCGTACCTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCa                                                                   <  2:372408/150‑1 (MQ=255)
                                                                                 cggaaaCGGACCCGTGGCGCGTGGTATGGGCCTCTCGGCCCCCGGGGCGTCCCTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAca                                                                    <  2:372390/144‑2 (MQ=255)
                                                                                 ccGGAACTGACCCGCGGGGCTCTGGATGTGCCTCTCGGACACCGGGGCGTACTTTTTCCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                  <  1:793631/150‑1 (MQ=255)
                                                                                     aaCTGACCCGGGGCGCGTTGGATGTGGCCCGCGGACACCTGGGCTTTCCTTTTACCTGTTTCGGGAAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                             <  1:115147/151‑1 (MQ=255)
                                                                                     aaCTAACCCGTGGCGCGTTGGTGGGGCCTCTCGGACCCCGGGGCGTCCCTTTTCCTTTTATCGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                             <  1:909806/151‑1 (MQ=255)
                                                                                        cTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc                                                            >  2:897582/1‑150 (MQ=255)
                                                                                          gAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                         >  1:161737/1‑151 (MQ=255)
                                                                                           agctggggCGCGTGGGTTGGGGCTCCCGGCCCCTTGGGCGTCCCGTTTCCTTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                          <  2:937551/143‑1 (MQ=255)
                                                                                           aaCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                        >  1:863315/1‑151 (MQ=255)
                                                                                           aaCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCTTGTTGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCa                                                           <  2:767644/148‑1 (MQ=255)
                                                                                                    gcgcTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                   >  1:412068/1‑147 (MQ=255)
                                                                                                      ggtGGATGTGGCTCTCGGCCACGGGGGCGTCCCTTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAg                                                 <  2:336129/145‑1 (MQ=255)
                                                                                                           aTGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgcg                                         >  2:430222/1‑150 (MQ=255)
                                                                                                               ggCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa                                      >  2:668097/1‑149 (MQ=255)
                                                                                                               ggCTCTCGGACACCGGGGCGTCCCTTTTCCTTTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaa                                       <  2:781084/148‑1 (MQ=255)
                                                                                                                     cGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                              >  2:633000/1‑151 (MQ=255)
                                                                                                                      ggccacgggggCGTCCCGTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                              <  2:509139/142‑1 (MQ=255)
                                                                                                                      ggACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                             >  2:3223/1‑151 (MQ=255)
                                                                                                                          acCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                            >  2:552764/1‑148 (MQ=255)
                                                                                                                          acCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                            >  2:637881/1‑148 (MQ=255)
                                                                                                                           ccggggCGTCCCTTTTCCTTGTTTCGGCCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAg                          <  2:379071/146‑1 (MQ=255)
                                                                                                                           ccTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                            >  1:545231/1‑147 (MQ=255)
                                                                                                                              gggCGTCCCTTTTCCTTGTATGGGAAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGaa                           <  2:524983/145‑1 (MQ=255)
                                                                                                                              gggCGTACCTTTTCCCTTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGAt                        <  2:518242/148‑1 (MQ=255)
                                                                                                                               ggCGTCCCTTTTCCCTTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                    <  1:426509/151‑1 (MQ=255)
                                                                                                                               ggCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                    >  2:185906/1‑151 (MQ=255)
                                                                                                                                    tcccttttCCCTGTAGGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTa                   <  1:856256/143‑1 (MQ=255)
                                                                                                                                    aCCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc               >  1:356342/1‑151 (MQ=255)
                                                                                                                                       gTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg            >  1:611423/1‑151 (MQ=255)
                                                                                                                                       gTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg            >  1:937146/1‑151 (MQ=255)
                                                                                                                                         tttACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTgg           >  1:25114/1‑151 (MQ=255)
                                                                                                                                                  aTCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt  >  2:412075/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                
CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTT  >  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: