Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,473,927 C→A A43S (GCC→TCC)  bdcA ← c‑di‑GMP‑binding biofilm dispersal mediator protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,473,9270CA98.2% 409.3 / ‑1.6 114A43S (GCC→TCC) bdcAc‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (58/54);  total (60/54)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 114 total aligned reads displayed.

GCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGAT  >  NC_000913/4473830‑4474023
                                                                                                 |                                                                                                
gcgcCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAAt                                                                                                <  1:1358594/100‑1 (MQ=255)
 cgccgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACGGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATc                                                                                               >  1:1705688/1‑100 (MQ=255)
    cgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCttt                                                                                            <  2:1517497/100‑1 (MQ=255)
     gcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTc                                                                                           >  1:2997903/1‑100 (MQ=255)
      ctctTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCg                                                                                          <  1:2482190/100‑1 (MQ=255)
      ctctTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCg                                                                                          <  2:1488389/100‑1 (MQ=255)
      ctctTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCg                                                                                          <  2:1487328/100‑1 (MQ=255)
         ttACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAcc                                                                                       <  1:2218656/100‑1 (MQ=255)
         ttACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAcc                                                                                       >  2:2897747/1‑100 (MQ=255)
          tACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAccc                                                                                      >  2:1767637/1‑100 (MQ=255)
          tACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAccc                                                                                      >  2:150681/1‑100 (MQ=255)
             gaacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGc                                                                                   <  1:2781558/100‑1 (MQ=255)
             gaacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGc                                                                                   <  2:2623602/100‑1 (MQ=255)
             gaacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGc                                                                                   <  1:2781611/100‑1 (MQ=255)
              aacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCa                                                                                  >  2:2435977/1‑100 (MQ=255)
               acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCAt                                                                                 <  1:346539/100‑1 (MQ=255)
                cgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATa                                                                                <  1:2931293/100‑1 (MQ=255)
                   cATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGt                                                                             >  2:2636722/1‑100 (MQ=255)
                    aTCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTg                                                                            >  1:1929666/1‑100 (MQ=255)
                      tgctGACAGAGACGCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGaa                                                                          >  1:1027322/4‑100 (MQ=255)
                      cAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGaa                                                                          >  2:357342/1‑100 (MQ=255)
                        aTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATc                                                                        >  1:1020864/1‑100 (MQ=255)
                        aTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATc                                                                        >  1:1021025/1‑100 (MQ=255)
                         tGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCg                                                                       <  1:128415/100‑1 (MQ=255)
                         tGACAGCGTCTCTGTAAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCg                                                                       <  2:2200855/100‑1 (MQ=255)
                          gACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGt                                                                      >  1:2034746/1‑100 (MQ=255)
                            cAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                    <  1:1333604/100‑1 (MQ=255)
                            cAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                    <  2:1131189/100‑1 (MQ=255)
                            cAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                    <  1:256865/100‑1 (MQ=255)
                            cAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCCCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                    <  1:2127363/100‑1 (MQ=255)
                            cAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCCCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                    <  1:2127381/100‑1 (MQ=255)
                               cGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACAtt                                                                 >  1:2451808/1‑100 (MQ=255)
                               cGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACAtt                                                                 <  1:167598/100‑1 (MQ=255)
                               cGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACAtt                                                                 <  1:168052/100‑1 (MQ=255)
                                 tctcTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTgg                                                               <  1:1542008/100‑1 (MQ=255)
                                   tctGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                             <  1:822853/100‑1 (MQ=255)
                                   tctGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                             <  1:1897894/100‑1 (MQ=255)
                                   tctGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                             <  2:2116182/100‑1 (MQ=255)
                                   tctGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                             <  1:1591837/100‑1 (MQ=255)
                                   tctGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                             >  2:1237910/1‑100 (MQ=255)
                                     tGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcccc                                                           >  1:2237677/1‑100 (MQ=255)
                                       tCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAt                                                         <  2:1166012/100‑1 (MQ=255)
                                          gCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCgg                                                      >  2:1471922/1‑100 (MQ=255)
                                          gCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCgg                                                      >  2:364347/1‑100 (MQ=255)
                                          gCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCgg                                                      >  1:3015666/1‑100 (MQ=255)
                                            aCTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTg                                                    >  1:245645/1‑100 (MQ=255)
                                             cTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCCTAGGTGCATCGTACATTGGCCCCATCCGTGa                                                   >  2:633458/1‑100 (MQ=255)
                                                  tgtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGaaa                                              >  2:2689257/1‑100 (MQ=255)
                                                  tgtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGaaa                                              >  1:1236869/1‑100 (MQ=255)
                                                  tgtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGaaa                                              >  1:1658549/1‑100 (MQ=255)
                                                  tgtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGaaa                                              >  1:1966986/1‑100 (MQ=255)
                                                   gtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAac                                             >  2:7815/1‑100 (MQ=255)
                                                   gtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAac                                             <  2:1104162/100‑1 (MQ=255)
                                                   gtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAac                                             <  2:1103964/100‑1 (MQ=255)
                                                   gtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAac                                             >  1:708381/1‑100 (MQ=255)
                                                   gtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAac                                             >  1:2483410/1‑100 (MQ=255)
                                                     gAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacga                                           <  2:1256621/100‑1 (MQ=255)
                                                      aaTACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacgac                                          <  1:490576/100‑1 (MQ=255)
                                                        tACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGt                                        >  2:390839/1‑100 (MQ=255)
                                                         aCTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTa                                       <  1:1767637/100‑1 (MQ=255)
                                                          ctgctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTTGCCCCATCGGTGACGAAACGACGTAc                                      >  2:2266500/1‑100 (MQ=255)
                                                           tgctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACg                                     <  1:1130428/100‑1 (MQ=255)
                                                            gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGa                                    >  2:263073/1‑100 (MQ=255)
                                                               gTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAg                                 <  1:2725759/100‑1 (MQ=255)
                                                               gTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAg                                 >  2:2524106/1‑100 (MQ=255)
                                                                 cGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCg                               <  1:2624925/100‑1 (MQ=255)
                                                                  gCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCgg                              <  1:1257612/100‑1 (MQ=255)
                                                                  gCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCgg                              >  1:692673/1‑100 (MQ=255)
                                                                  gCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCgg                              >  1:147900/1‑100 (MQ=255)
                                                                  gCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCgg                              >  1:1730505/1‑100 (MQ=255)
                                                                      cAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCAcc                          <  2:1522925/100‑1 (MQ=255)
                                                                         tctcTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGAt                       <  1:1079429/100‑1 (MQ=255)
                                                                         tctcTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGAt                       <  1:2992881/100‑1 (MQ=255)
                                                                          ctctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATa                      <  2:1979256/100‑1 (MQ=255)
                                                                           tctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAc                     >  2:1920620/1‑100 (MQ=255)
                                                                           tctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAc                     >  1:2795039/1‑100 (MQ=255)
                                                                            ctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAcc                    >  1:3120563/1‑100 (MQ=255)
                                                                              tgtgCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGCCGTACGATAGCGGCACCGATACCAc                  >  2:700559/1‑100 (MQ=255)
                                                                              tgtgCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCAc                  >  2:700939/1‑100 (MQ=255)
                                                                               gtgCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACg                 >  2:762928/1‑100 (MQ=255)
                                                                                tgCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGa                >  2:1342888/1‑100 (MQ=255)
                                                                                 gCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGAc               >  1:2620123/1‑100 (MQ=255)
                                                                                 gCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGAc               <  1:167198/100‑1 (MQ=255)
                                                                                  ccAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACt              >  2:1748885/1‑100 (MQ=255)
                                                                                   cAGGCGTTTAGCGGAATCTTTCGACCCCGCATCGGTGAATCGTACATTGGCCCCATCGGTGACGCAACGACGTACGATAGCGGCACCGATACCACGTCTg             >  2:3039790/1‑100 (MQ=255)
                                                                                   cAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTg             >  1:959384/1‑100 (MQ=255)
                                                                                   cAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTg             <  2:182082/100‑1 (MQ=255)
                                                                                    aGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGc            <  1:2000124/100‑1 (MQ=255)
                                                                                       cGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGCTAGCGGCACCGATACCACGACTGCCAc         >  2:3096010/1‑100 (MQ=255)
                                                                                       cGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCAc         >  2:645124/1‑100 (MQ=255)
                                                                                       cGTTTAGCGGAATCTTTCGACCCCGAATAGGCGACTCGTACATTGGTCCCATCGGTTACGAAACGACGTACGATAGCGGCACCGATACCATGACTGCCAc         >  2:3095929/1‑100 (MQ=255)
                                                                                        gTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCAcc        <  2:606548/100‑1 (MQ=255)
                                                                                        gTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCAcc        >  1:695668/1‑100 (MQ=255)
                                                                                         tttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCg       <  1:1063839/100‑1 (MQ=255)
                                                                                         tttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCg       <  1:1063820/100‑1 (MQ=255)
                                                                                          ttAGCGGAATCGTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAGCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCgn      <  2:712032/100‑2 (MQ=255)
                                                                                            aGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAgg    <  2:1756309/100‑1 (MQ=255)
                                                                                             gCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGa   <  2:2290906/100‑1 (MQ=255)
                                                                                             gCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGa   >  1:399507/1‑100 (MQ=255)
                                                                                              cGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGAt  >  1:2627264/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGAT  >  NC_000913/4473830‑4474023

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