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

TTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGG  >  NC_000913/436355‑436651
                                                                                                                                                       |                                                                                                                                                   
tttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTa                                                                                                                                                     <  1:602701/150‑1 (MQ=255)
tttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTa                                                                                                                                                     <  2:665964/150‑1 (MQ=255)
 ttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc                                                                                                                                                   >  1:39898/1‑151 (MQ=255)
     gATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTg                                                                                                                                                 <  2:606415/149‑1 (MQ=255)
      aTGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTAt                                                                                                                                              >  2:303228/1‑151 (MQ=255)
          gcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCggg                                                                                                                                          >  2:217177/1‑151 (MQ=255)
          gcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCggg                                                                                                                                          >  2:557018/1‑151 (MQ=255)
           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)
                                                                                                                                                       |                                                                                                                                                   
TTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGG  >  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: