Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,163,678 A→C 9.2% D268A (GAT→GCT)  baeS → sensory histidine kinase BaeS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,163,6780AC 9.2% 308.6 / 24.9 87D268A (GAT→GCT) baeSsensory histidine kinase BaeS
Reads supporting (aligned to +/- strand):  ref base A (45/34);  new base C (6/2);  total (51/36)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GAAAAACCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATG  >  NC_000913/2163580‑2163775
                                                                                                  |                                                                                                 
gaAAAACCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGa                                                                                                   <  1:136111/99‑1 (MQ=255)
gaAAAACCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGa                                                                                                   >  2:371459/1‑99 (MQ=255)
  aaaaCCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGCTg                                                                                                 <  2:1708594/99‑1 (MQ=255)
  aaaaCCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATg                                                                                                 >  2:160425/1‑99 (MQ=255)
  aaaaCCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATg                                                                                                 >  1:1433324/1‑99 (MQ=255)
   aaaCCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATgg                                                                                                >  1:2026322/1‑99 (MQ=255)
   aaaCCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATgg                                                                                                <  2:1595896/99‑1 (MQ=255)
    aaCCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGCTGgc                                                                                               >  2:148637/1‑99 (MQ=255)
    aaCCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGgc                                                                                               >  2:1630258/1‑99 (MQ=255)
     aCCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGCTGgcg                                                                                              >  1:1180069/1‑99 (MQ=255)
     aCCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGgcg                                                                                              <  2:1900652/99‑1 (MQ=255)
      ccAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGgcgt                                                                                             <  1:665915/99‑1 (MQ=255)
       cAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGgcgtg                                                                                            <  1:889946/99‑1 (MQ=255)
         gCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGgcgtgcg                                                                                          >  2:869054/1‑99 (MQ=255)
             aTGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTaaa                                                                                      >  2:2394442/1‑99 (MQ=255)
               gcggcgCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAAtt                                                                                    <  2:1312723/99‑1 (MQ=255)
                 ggcgCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGCTGGCGTGCGTAAATTCa                                                                                  <  1:1082781/99‑1 (MQ=255)
                 ggcgCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCa                                                                                  <  1:674962/99‑1 (MQ=255)
                  gcgcGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCAc                                                                                 >  1:659961/1‑99 (MQ=255)
                    gcgATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGc                                                                               >  1:1411415/1‑99 (MQ=255)
                     cgATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGcc                                                                              >  2:2180294/1‑99 (MQ=255)
                      gATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCg                                                                             >  2:504786/1‑99 (MQ=255)
                      gATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCg                                                                             >  2:1348537/1‑99 (MQ=255)
                       aTTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGCTGGCGTGCGTAAATTCACGCCgg                                                                            >  1:2152269/1‑99 (MQ=255)
                       aTTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCgg                                                                            >  2:2412558/1‑99 (MQ=255)
                       aTTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCgg                                                                            >  1:1159479/1‑99 (MQ=255)
                        ttttATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGga                                                                           <  2:1506420/99‑1 (MQ=255)
                          ttATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGgaga                                                                         <  2:2310882/99‑1 (MQ=255)
                           tATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGAc                                                                        <  2:767529/99‑1 (MQ=255)
                            aTGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACg                                                                       >  2:1440715/1‑99 (MQ=255)
                             tGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACgg                                                                      <  2:197922/99‑1 (MQ=255)
                                    tatTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTc                                                               >  2:1180296/1‑99 (MQ=255)
                                         ctcACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTAc                                                          <  1:837468/99‑1 (MQ=255)
                                          tcACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACa                                                         >  1:286545/1‑99 (MQ=255)
                                             cGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGc                                                      <  2:1910341/99‑1 (MQ=255)
                                             cGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGc                                                      <  1:1970955/99‑1 (MQ=255)
                                              gAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCg                                                     >  2:713812/1‑99 (MQ=255)
                                              gAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCg                                                     <  2:1671093/99‑1 (MQ=255)
                                               aaCTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGCTGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCgg                                                    >  2:1291298/1‑99 (MQ=255)
                                               aaCTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCgg                                                    >  2:422228/1‑99 (MQ=255)
                                                aCTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGa                                                   >  2:671806/1‑99 (MQ=255)
                                                  tGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAgg                                                 <  1:895051/99‑1 (MQ=255)
                                                   gCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGt                                                <  2:2212923/99‑1 (MQ=255)
                                                   gCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGt                                                <  1:2323889/99‑1 (MQ=255)
                                                     gTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCg                                              >  1:1839706/1‑99 (MQ=255)
                                                     gTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCg                                              <  2:1233538/99‑1 (MQ=255)
                                                      tACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGCTGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCgg                                             >  2:1783121/1‑99 (MQ=255)
                                                      tACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCgg                                             >  2:341127/1‑99 (MQ=255)
                                                      tACGCCAGTAGCGGTGTCGCGCGGTGAACTGGAGGCCAATCAGGATGGCGTGCGTAAGTTCCCGCTGGAGAAGGTGGCTTCGTTACAGGATGAGGTCgg                                             <  2:1854101/99‑1 (MQ=255)
                                                       aCGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGt                                            >  2:1749509/1‑99 (MQ=255)
                                                        cGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTa                                           <  2:210929/99‑1 (MQ=255)
                                                        cGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTa                                           <  1:528709/99‑1 (MQ=255)
                                                          ccATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTaca                                         <  1:407630/99‑1 (MQ=255)
                                                          ccATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTaca                                         >  1:1726313/1‑99 (MQ=255)
                                                           cATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTacac                                        <  2:1913718/99‑1 (MQ=255)
                                                            aTTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACt                                       <  2:2348788/99‑1 (MQ=255)
                                                             ttAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTg                                      >  1:1834561/1‑99 (MQ=255)
                                                                gCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGAcc                                   >  2:86110/1‑99 (MQ=255)
                                                                   gTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCaaa                                >  2:1526158/1‑99 (MQ=255)
                                                                   gTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCaaa                                >  1:953818/1‑99 (MQ=255)
                                                                    tgctgcGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAAc                               >  2:995915/1‑99 (MQ=255)
                                                                      ctgcGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTg                             <  2:522313/99‑1 (MQ=255)
                                                                       tgcGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTgg                            <  2:915335/99‑1 (MQ=255)
                                                                       tgcGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTgg                            >  1:710295/1‑99 (MQ=255)
                                                                        gcgcGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGt                           <  2:738131/99‑1 (MQ=255)
                                                                         cgcgGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGtt                          >  2:1427446/1‑99 (MQ=255)
                                                                         cgcgGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGtt                          <  1:592136/99‑1 (MQ=255)
                                                                           cgGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGa                        >  1:1277241/1‑99 (MQ=255)
                                                                               gAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGAt                    >  1:324440/1‑99 (MQ=255)
                                                                               gAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGAt                    >  1:1161676/1‑99 (MQ=255)
                                                                                  cTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGAtctc                 <  1:348986/99‑1 (MQ=255)
                                                                                  cTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGAtctc                 >  1:426580/1‑99 (MQ=255)
                                                                                  cTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGAtctc                 >  1:2034555/1‑99 (MQ=255)
                                                                                     gAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCAt              <  1:183911/99‑1 (MQ=255)
                                                                                       aGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCa            <  1:929968/99‑1 (MQ=255)
                                                                                       aGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCa            >  1:606523/1‑99 (MQ=255)
                                                                                       aGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCa            >  1:671331/1‑99 (MQ=255)
                                                                                       aGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCa            >  1:1060277/1‑99 (MQ=255)
                                                                                        gCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAg           <  1:485732/99‑1 (MQ=255)
                                                                                          cATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGtt         >  2:1307434/1‑99 (MQ=255)
                                                                                           aTTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTg        >  1:881045/1‑99 (MQ=255)
                                                                                            ttCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGt       >  2:446177/1‑99 (MQ=255)
                                                                                            ttCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGt       >  1:1765761/1‑99 (MQ=255)
                                                                                            ttCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGt       <  2:2155079/99‑1 (MQ=255)
                                                                                            ttCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGt       >  1:2169631/1‑99 (MQ=255)
                                                                                             tCAGGCTGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTc      >  1:872201/1‑99 (MQ=255)
                                                                                             tCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTc      <  1:1871522/99‑1 (MQ=255)
                                                                                               aGGATGGCGGGCGTAAATTCACGCCGGAGACGGAGGCGTCTTTACAGGCGGAGCTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGa    >  1:1376216/1‑99 (MQ=255)
                                                                                                 gATGGCGTGCTTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGTCGGAGGCCGGTACACTGACCAAACCGGTTTACGATCTCCATCAGTTGTCGAtt  >  2:1038286/1‑98 (MQ=255)
                                                                                                  |                                                                                                 
GAAAAACCAGCAAATGCGGCGCGATTTTATGGCCGATATTTCTCACGAACTGCGTACGCCATTAGCGGTGCTGCGCGGTGAACTGGAAGCCATTCAGGATGGCGTGCGTAAATTCACGCCGGAGACGGTGGCGTCTTTACAGGCGGAGGTCGGTACACTGACCAAACTGGTTGACGATCTCCATCAGTTGTCGATG  >  NC_000913/2163580‑2163775

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: