Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5160CA41.5% 66.0 / 27.7 54Q295K (CAA→AAA) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (27/4);  new base A (0/22);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 3.23e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGC  >  NC_000913/436366‑436662
                                                                                                                                                            |                                                                                                                                                    
cTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGc                                                                                                                                                      >  1:418175/1‑151 (MQ=255)
   tctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGAAACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑aa                                                                                                                                                   >  2:220474/1‑150 (MQ=255)
      cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑Ga                                                                                                                                              <  2:28845/151‑1 (MQ=255)
        ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GAcc                                                                                                                                            <  2:53046/151‑1 (MQ=255)
            gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GAccgcc                                                                                                                                         <  1:813288/150‑1 (MQ=255)
                   cGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCg                                                                                                                                  >  2:740228/1‑150 (MQ=255)
                    ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGa                                                                                                                                 >  1:654476/1‑150 (MQ=255)
                    ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑GACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGa                                                                                                                                 >  2:692103/1‑150 (MQ=255)
                          ggggcggcggggggCGCTCTCGGGGGGAGAAGATTGCGGTTTGTTTTTCCCTGTCCCGA‑AACTGACCCGGGGCGCGTT‑GGGGGTGCCTCCCGCCCCCTG‑GGGGG‑TCCT‑T‑TTTACCTTTTGGGGAA‑AAAT‑GACCGCCGATATCGAAGGGCtt                                                                                                                          <  2:1180320/140‑1 (MQ=255)
                           gcgctgcgctGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCttt                                                                                                                         >  1:402729/1‑151 (MQ=255)
                                 cgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                    >  1:979128/1‑150 (MQ=255)
                                 cgcggGGCGCCCTCGGGGGGTGAAGATTCGGGTTTGTTTTTCCCTGTGCCGG‑AACTGACCCGGGGCGCGCT‑GGATGTGGCTCCCGGACCCTG‑GGGCG‑TCCT‑T‑TTTTCTTGTACCGGAA‑AATT‑GACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                    <  2:1106924/150‑1 (MQ=255)
                                  ggtgGGCCTTCTCGGGCGGTGAAGTTTCCGATTTGTGTTTCCCTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCCCTCGGCCACCT‑GGGGG‑TCCC‑T‑TTTCCTTTTATGGGGA‑AAAA‑GACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                  <  1:817194/149‑1 (MQ=255)
                                     gggcctctgcggggggTGAAGATTAGGGTTTTTGTTTCCCTGGGCCGG‑AACTGAACGGTGGCGCGTT‑GGTTGTGCCCCGCGCAC‑CCG‑GGGGG‑TTCCTT‑TTTCCCTTTAGCGGAA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                 <  2:1078737/135‑1 (MQ=255)
                                     gggCGCTCTGTGGCGGGAAAGTTTACGTGTTTTTTTTCTGTGTGCGGG‑AACTGACCGGGGGCGCGCT‑GGGTGTGGCTCTCGGACCCTG‑GGGCG‑TCCT‑G‑TTTCCCTGTGGCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCttt                                                                                                                         <  2:254210/141‑1 (MQ=255)
                                        cgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                             >  1:1168517/1‑150 (MQ=255)
                                              ctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                        >  1:804978/1‑149 (MQ=255)
                                                ggggttGAAGATTCGGATTTGTTTTGTATGGTCCCGG‑AACTGCTCCGGGGCGCGCG‑GGTTGTCTCCTTCGG‑CGACG‑GGGCT‑TTCC‑GTTTTCCCTTTGTCGGAA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                                     <  2:784792/145‑1 (MQ=255)
                                                ggCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGACCCGGGGCGCGTT‑GGTTGTGGCTCTCGGACACTG‑GGGCG‑TACC‑G‑TTTCCTTTTATGGGGA‑AAAT‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                        <  1:1032022/147‑1 (MQ=255)
                                                ggCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                                     >  2:203604/1‑150 (MQ=255)
                                                   ggtgaagtttcgggtttGTGTTTCACTGTGCGGG‑AACTGACCGGGGGCGCGTG‑GGATGGGTCTCTCGGACCCCG‑GGGCG‑TCCT‑T‑TTTCCTTGTAGCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                 <  2:815475/136‑1 (MQ=255)
                                                   ggTAAAAATTTCGATTTGTTTTTCACTGTGCGAG‑AACTGACCTGCGGCCCGTG‑GGGTGGGGCTCGCGGCCACGT‑GGGCG‑TACT‑T‑TTTACCTGTACCGGGA‑AAGT‑CACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                  <  2:349134/150‑1 (MQ=255)
                                                        aGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                             >  2:791576/1‑150 (MQ=255)
                                                            taggAGTTTTTTTTCCCTGGGCCGG‑AACTGACCGGGGGCGCGTT‑GGTTGTGGCTCTGGGACCCCG‑GGGCG‑TCCC‑T‑TTTCCTTGTATCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                                         <  1:1050362/147‑1 (MQ=255)
                                                            tacgttttGTGTTTCCCTGTGCCGG‑AACGGACCGGGGGCGCGTT‑GGTGGTGCCTCCGGGACACCT‑GGGCT‑GCCC‑T‑TTTCCCTTTGCGGGAA‑AAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                                         <  2:435075/144‑1 (MQ=255)
                                                            tACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGCCACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                          >  2:2904/1‑149 (MQ=255)
                                                              cGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc                                                                                      >  2:140274/1‑151 (MQ=255)
                                                                ggttgttTTTCCTTGTGCCGA‑AACGGACCCGTGGCGCGGT‑GGGTGGGGCTC‑CGGGAACCTGGGGCG‑TCCC‑T‑TTTTCCTGTACCGGGA‑AAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                                        <  1:634064/146‑1 (MQ=255)
                                                                   tgtgTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                  >  1:415923/1‑150 (MQ=255)
                                                                    tttttttcccgggCCGGAAACTGACCCGGGGCGCGTT‑GGTTGTGCCTCTCGGCCCCCG‑GGGCG‑TCCC‑T‑TTTCCTTGTTTCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                  <  1:1086262/138‑1 (MQ=255)
                                                                          ctcTGTGCGGG‑AAC‑GAACCGGGGGGCGGTGGGGTGTGCCTC‑CGGGCCCCT‑GGGGGTTCCC‑T‑TTTACCTGTATGGGAA‑AAAG‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                           <  2:189651/148‑1 (MQ=255)
                                                                          cACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                          >  2:788386/1‑151 (MQ=255)
                                                                               tgCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                                      >  1:380946/1‑150 (MQ=255)
                                                                                    gaaaCTGACCCTTGGCGGGTT‑GGTTGTGTCTCCGGAACCCCT‑GGGCG‑T‑CC‑C‑TTTTCCTTTTCCGGGG‑AAATGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                                   <  1:396090/147‑1 (MQ=255)
                                                                                      accggcccggTGGCGCGTT‑GGATGTGGCTCCCGGCCACCT‑GGGCG‑TACC‑T‑TTTCCCTGTATGGGGA‑AAAT‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                 <  2:262218/140‑1 (MQ=255)
                                                                                      aaCTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                               >  1:468079/1‑151 (MQ=255)
                                                                                         tGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                            >  2:200081/1‑151 (MQ=255)
                                                                                          gAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                           >  2:261005/1‑151 (MQ=255)
                                                                                              cgtgcggcgcT‑GGAGGTGCTCCTCGTACCCTG‑GGGCG‑TACT‑T‑TTTCCTTTTGTCGGAA‑AAAT‑AACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                           <  2:1222531/141‑1 (MQ=255)
                                                                                                 ggCGGGTT‑GGGTGTGCCTCTCGTCCACCT‑GGGTT‑TCCC‑G‑TTTCCTTGTATCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                            <  1:703850/143‑1 (MQ=255)
                                                                                                     cgcT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGTCGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                                >  2:322051/1‑151 (MQ=255)
                                                                                                      ggt‑GGATGTGTCTCTCGGCCCCTG‑GGGCG‑TCCT‑T‑TTTCCTTGTAGCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                                <  2:243903/148‑1 (MQ=255)
                                                                                                      gcT‑GGATGTGGCTATCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                                >  2:651959/1‑150 (MQ=255)
                                                                                                        t‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                             >  1:1188573/1‑151 (MQ=255)
                                                                                                           gttggtgcctccGGACCCCTGGGGCG‑TCCC‑T‑TTTACCTTTA‑CGGGCAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                              <  2:95360/138‑1 (MQ=255)
                                                                                                               tgGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa                                        >  2:937685/1‑150 (MQ=255)
                                                                                                                ggCTCTGGGCCCCCG‑GGGCG‑TCCC‑T‑TTTCCTTTTAGGGGGA‑AAAT‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaaa                                       <  1:284147/150‑1 (MQ=255)
                                                                                                                    ctcGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc                                  >  1:851242/1‑151 (MQ=255)
                                                                                                                             cgggggCG‑TCCC‑T‑TTTCCTTGTAGGGGAA‑AAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGa                           <  1:1132430/147‑1 (MQ=255)
                                                                                                                                             tttCCCTGTATCGGAA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg                <  1:920961/148‑1 (MQ=255)
                                                                                                                                              ttACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc            >  1:433743/1‑151 (MQ=255)
                                                                                                                                                  cTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt        >  1:1273077/1‑151 (MQ=255)
                                                                                                                                                     tATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt     >  1:637848/1‑151 (MQ=255)
                                                                                                                                                        cGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGc  >  2:1178118/1‑151 (MQ=255)
                                                                                                                                                         gggA‑AAAT‑GCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt    <  1:140276/148‑1 (MQ=255)
                                                                                                                                                            |                                                                                                                                                    
CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGG‑AACTGAACCGTGGCGCGCT‑GGATGTGGCTCTCGGACACCT‑GGGCG‑TACC‑G‑TTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGC  >  NC_000913/436366‑436662

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑