Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,456,659 G→A 5.2% R292H (CGC→CAC)  mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,456,6590GA 5.2% 247.5 / 2.2 77R292H (CGC→CAC) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base G (33/39);  new base A (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGCCCATCACCCGACGGCGATTCTGG  >  NC_000913/4456513‑4456799
                                                                                                                                                  |                                                                                                                                            
gCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGc                                                                                                                                             >  8:29261/1‑148 (MQ=255)
    ggTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGc                                                                                                                                             >  3:165583/1‑144 (MQ=255)
    ggTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGc                                                                                                                                             <  4:165583/144‑1 (MQ=255)
    ggTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGt                                                                                                                                        >  5:305114/1‑149 (MQ=255)
    ggTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGt                                                                                                                                        <  2:147310/149‑1 (MQ=255)
      tCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGTCGATTGCAGCGGCTCGCCATGTTg                                                                                                                                      <  6:122078/149‑1 (MQ=255)
        aCTGGCAGGCGAAAAAGCTGACCACCGATGCTTCAGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGGGAACATAATTTGCACAATGGCCTGATCTCGAATGCAGCGGCTCGCCATGTTGgt                                                                                                                                    >  3:31105/1‑149 (MQ=255)
            gCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAg                                                                                                                                >  5:257082/1‑149 (MQ=255)
                gCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGcc                                                                                                                            <  2:75773/149‑1 (MQ=255)
                gCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGcc                                                                                                                            >  2:229246/1‑149 (MQ=255)
                  gAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCgg                                                                                                                          >  8:18204/1‑149 (MQ=255)
                    aaaaGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCG‑AAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAg                                                                                                                       <  4:219041/149‑1 (MQ=255)
                      aaGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGa                                                                                                                      <  3:260845/149‑1 (MQ=255)
                           gACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCg                                                                                                                 >  1:19935/1‑149 (MQ=255)
                               accGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAgc                                                                                                                               <  2:129742/131‑1 (MQ=255)
                               accGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAgc                                                                                                                               >  1:129742/1‑131 (MQ=255)
                               accGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTaa                                                                                                             <  8:140199/149‑1 (MQ=255)
                                 cGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCACCATGTTGGTGTAGCGCCGGCAGATGCCGCTAAcg                                                                                                           >  5:108709/1‑149 (MQ=255)
                                 cGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAACTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCACCATGTTGGTGTAGCGCCGGCAGATGCCGCTAAcg                                                                                                           >  1:288192/1‑149 (MQ=255)
                                       ttCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCa                                                                                                                        <  7:201535/130‑1 (MQ=255)
                                       ttCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCa                                                                                                                        >  8:201535/1‑130 (MQ=255)
                                          cGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGtt                                                                                                  >  3:274910/1‑149 (MQ=255)
                                             aTGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGt                                                                                               >  6:75748/1‑149 (MQ=255)
                                                 gAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTttat                                                                                           <  1:229246/149‑1 (MQ=255)
                                                  aaGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCACCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTttatt                                                                                          <  2:288192/149‑1 (MQ=255)
                                                  aaGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCACCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTttatt                                                                                          <  6:108709/149‑1 (MQ=255)
                                                      tttGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGa                                                                                                                      <  2:159589/117‑1 (MQ=255)
                                                      tttGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGa                                                                                                                      >  1:159589/1‑117 (MQ=255)
                                                         gCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATg                                                                                      <  7:111802/146‑1 (MQ=255)
                                                         gCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATg                                                                                      >  8:111802/1‑146 (MQ=255)
                                                          cTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCtc                                                                                   <  6:257082/148‑1 (MQ=255)
                                                           tGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCtcgt                                                                                 >  8:44770/1‑149 (MQ=255)
                                                           tGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCtcgt                                                                                 >  3:63510/1‑149 (MQ=255)
                                                                   gAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCgg                                                                                                                          >  2:402140/1‑100 (MQ=255)
                                                                   gAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCgg                                                                                                                          <  1:402140/100‑1 (MQ=255)
                                                                         gTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCTGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGtt                                                                   <  8:392645/149‑1 (MQ=255)
                                                                         gTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGtt                                                                   >  4:256263/1‑149 (MQ=255)
                                                                         gTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGtt                                                                   <  7:44770/149‑1 (MQ=255)
                                                                             gCGATGTGATATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGTGATTTCAGCGGCTCTCCATGTTGGTTTAGCGCCGGAAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCgt                                                               <  5:329780/149‑1 (MQ=255)
                                                                              cGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCgtg                                                              >  1:125040/1‑149 (MQ=255)
                                                                              cGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCgtg                                                              <  7:18204/149‑1 (MQ=255)
                                                                               gAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCgtgg                                                             >  3:374416/1‑149 (MQ=255)
                                                                                aaGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCgtggt                                                            <  6:305114/149‑1 (MQ=255)
                                                                                aaGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCgtggt                                                            <  6:181466/149‑1 (MQ=255)
                                                                                     aaaTGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGc                                                       <  4:374416/149‑1 (MQ=255)
                                                                                     aaaTGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGc                                                       >  5:172645/1‑149 (MQ=255)
                                                                                       aTGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGaa                                                         >  6:89662/1‑145 (MQ=255)
                                                                                       aTGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGaa                                                         <  5:89662/145‑1 (MQ=255)
                                                                                         ggTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTATCGCCGGCAGATGACGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAAt                                                   >  4:248878/1‑149 (MQ=255)
                                                                                          gTCGCTGGTAGGCGATCATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTACCGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATg                                                  >  1:97935/1‑149 (MQ=255)
                                                                                          gTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAAt                                                   <  7:218600/148‑1 (MQ=255)
                                                                                          gTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAAt                                                   <  7:29261/148‑1 (MQ=255)
                                                                                            cGCTGGAAGGCGTACATAATCTGCACAATGGCCTGATGGCGTTTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGc                                                <  8:372421/149‑1 (MQ=255)
                                                                                             gCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGc                                                <  7:352272/148‑1 (MQ=255)
                                                                                             gCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGc                                                <  6:172645/148‑1 (MQ=255)
                                                                                                 gTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGGCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCAc                                           <  4:151885/149‑1 (MQ=255)
                                                                                                       gAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGtata                                     <  4:63510/149‑1 (MQ=255)
                                                                                                       gAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGAGAATGGCGTCACGGtata                                     >  5:353119/1‑149 (MQ=255)
                                                                                                             aataTGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGa                               >  5:99304/1‑149 (MQ=255)
                                                                                                                 tGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGAtttt                           <  3:256263/149‑1 (MQ=255)
                                                                                                                    acaATGGTCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGc                                                       <  8:362836/118‑1 (MQ=255)
                                                                                                                    acaATGGTCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGc                                                       >  7:362836/1‑118 (MQ=255)
                                                                                                                       aTGTCCTTATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGTCAGATGCCGCTAACGCGCTGGGTTAGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGc                         <  2:270025/145‑1 (MQ=255)
                                                                                                                       aTGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTATCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGc                         >  1:270025/1‑145 (MQ=255)
                                                                                                                        tGGCCTGATGGCGATTTCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTTACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATTGCTTCACGGTATATGACGATTTTGCCCATc                    <  4:31105/149‑1 (MQ=255)
                                                                                                                        tGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGCCCATc                    <  5:78585/149‑1 (MQ=255)
                                                                                                                          gCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGCCCATCa                   <  5:75748/148‑1 (MQ=255)
                                                                                                                          gCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGACGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGCCCATCAc                  <  3:248878/149‑1 (MQ=255)
                                                                                                                            cTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAg                                                                     >  7:30081/1‑96 (MQ=255)
                                                                                                                            cTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAg                                                                     <  8:30081/96‑1 (MQ=255)
                                                                                                                             tGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAg                                                                     >  7:309223/1‑95 (MQ=255)
                                                                                                                             tGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAg                                                                     <  8:309223/95‑1 (MQ=255)
                                                                                                                             tGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGCCCATCAccc                >  6:338565/1‑148 (MQ=255)
                                                                                                                                   cGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGCCCATCACCCGACGGc          >  1:314433/1‑148 (MQ=255)
                                                                                                                                    gATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGGTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGCCCATCACCCGACGGCGa        <  2:19935/149‑1 (MQ=255)
                                                                                                                                         cAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGCCCATCACCCGACGGCGATTCTg   >  1:113573/1‑149 (MQ=255)
                                                                                                                                          aGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGCCCATCACCCGACGGCGATTCTgg  <  6:99304/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGATTGCAGCGGCTCGCCATGTTGGTGTAGCGCCGGCAGATGCCGCTAACGCGCTGGGTTCGTTTATTAATGCTCGTCGCCGTCTGGAGTTGCGTGGTGAAGCGAATGGCGTCACGGTATATGACGATTTTGCCCATCACCCGACGGCGATTCTGG  >  NC_000913/4456513‑4456799

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