Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,373,622 A→C intergenic (‑84/+26) pmrD ← / ← menE inactive two‑component system connector protein/O‑succinylbenzoate‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,373,6220AC98.5% 445.6 / ‑1.8 130intergenic (‑84/+26)pmrD/menEinactive two‑component system connector protein/O‑succinylbenzoate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (75/53);  total (77/53)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

Only 100 of 132 total aligned reads displayed.

CCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCC  >  NC_000913/2373526‑2373717
                                                                                                |                                                                                               
gcAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAAt                                                                                              <  2:2574198/99‑1 (MQ=255)
ccAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAAt                                                                                              >  2:74869/1‑100 (MQ=255)
ccAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAAt                                                                                              >  1:1332932/1‑100 (MQ=255)
 cAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATg                                                                                             >  2:1547793/1‑100 (MQ=255)
  aGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATga                                                                                            <  1:3925586/100‑1 (MQ=255)
  aGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATga                                                                                            >  2:3293394/1‑100 (MQ=255)
    ccattccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGCAATGGCAGCAATgaga                                                                                          <  2:3870388/100‑1 (MQ=255)
     cattccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGaa                                                                                         >  1:556374/1‑100 (MQ=255)
     cattccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGaa                                                                                         <  2:1128639/100‑1 (MQ=255)
      attccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAg                                                                                        >  2:3496387/1‑100 (MQ=255)
      attccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAg                                                                                        >  1:1530394/1‑100 (MQ=255)
      attccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCCTTATCAACGGCAATGGAATGGCAGCAATGAGAAg                                                                                        <  2:3693576/100‑1 (MQ=255)
       ttccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGt                                                                                       <  1:1560746/100‑1 (MQ=255)
       ttccattGCATTATCCTGTTTACTACGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGt                                                                                       <  1:3797130/100‑1 (MQ=255)
        tccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTg                                                                                      >  2:2095457/1‑100 (MQ=255)
         ccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGa                                                                                     >  2:22840/1‑100 (MQ=255)
         ccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGa                                                                                     >  1:1219977/1‑100 (MQ=255)
             tGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACg                                                                                 >  1:3843476/1‑100 (MQ=255)
             tGCATTATCCTGTTTGCTAAGAGGTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATTGCAGCAATGCGAAGTGAAACg                                                                                 <  2:3425067/100‑1 (MQ=255)
              gCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACgg                                                                                >  1:2827518/1‑100 (MQ=255)
                 ttATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAgg                                                                             >  2:1647189/1‑100 (MQ=255)
                   aTCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGt                                                                           >  1:1286098/1‑100 (MQ=255)
                   aTCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGt                                                                           <  1:3759553/100‑1 (MQ=255)
                   aTCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGt                                                                           >  1:882426/1‑100 (MQ=255)
                   aTCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGt                                                                           <  1:3161755/100‑1 (MQ=255)
                    tCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGtt                                                                          <  1:3111498/100‑1 (MQ=255)
                     ccTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTa                                                                         <  2:1118464/100‑1 (MQ=255)
                     ccTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTa                                                                         >  1:1692080/1‑100 (MQ=255)
                      cTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGACTGGCAGCAATGAGAAGTGAAACGGAGGGTTaa                                                                        >  1:2417411/1‑100 (MQ=255)
                       tGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAat                                                                       <  2:2598545/100‑1 (MQ=255)
                        gTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAatt                                                                      >  2:2529120/1‑100 (MQ=255)
                        gTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAatt                                                                      <  1:400429/100‑1 (MQ=255)
                           tGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAattatt                                                                   >  1:363398/1‑100 (MQ=255)
                           tGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAattatt                                                                   >  1:1334308/1‑100 (MQ=255)
                             cTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGc                                                                 <  1:2900449/100‑1 (MQ=255)
                             cTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGc                                                                 <  1:1464415/100‑1 (MQ=255)
                              tAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCt                                                                >  2:280704/1‑100 (MQ=255)
                                agagTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGa                                                              >  2:3108581/1‑100 (MQ=255)
                                   gTTTTCACATCCATTTACCGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGt                                                           <  2:513400/100‑1 (MQ=255)
                                   gTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGt                                                           >  2:2136850/1‑100 (MQ=255)
                                   cttgtCACATCAATTCCCAGATTAACAGATAAACTCTTTAACAACGGCCATGGAATGTGAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGt                                                           <  2:3616957/96‑1 (MQ=255)
                                    ttttCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGtt                                                          >  2:159293/1‑100 (MQ=255)
                                     tttCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTg                                                         >  2:3927719/1‑100 (MQ=255)
                                      ttCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGc                                                        <  1:2546372/100‑1 (MQ=255)
                                        cacaTCAATTTACGGATTAACAGATAAGGTCTTTATCAACGGCAATGGAATGGCAGCAAAGAGAAGTGAAACGGAGGGTTAATTATTGATGACGTTGccc                                                      >  1:3594788/1‑98 (MQ=255)
                                        cacaTCAATTTACGGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCAc                                                      >  1:2460271/1‑100 (MQ=255)
                                         acaTCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCAcc                                                     <  2:1942703/100‑1 (MQ=255)
                                          caTCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCAccc                                                    >  2:1843157/1‑100 (MQ=255)
                                          caTCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCAccc                                                    >  2:1317447/1‑100 (MQ=255)
                                          caTCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCAccc                                                    >  2:3259655/1‑100 (MQ=255)
                                           aTCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCa                                                   >  2:1120816/1‑100 (MQ=255)
                                            tCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCAc                                                  <  1:2728377/100‑1 (MQ=255)
                                            tCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCAc                                                  <  1:2568450/100‑1 (MQ=255)
                                               aTTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCActct                                               <  1:511857/100‑1 (MQ=255)
                                                 ttACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCttt                                             >  2:800373/1‑100 (MQ=255)
                                                    cAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAg                                          >  1:663697/1‑100 (MQ=255)
                                                     aGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAgc                                         >  2:2294735/1‑100 (MQ=255)
                                                      gATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAgcg                                        >  1:1763829/1‑100 (MQ=255)
                                                      gATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAgcg                                        <  1:144776/100‑1 (MQ=255)
                                                      gATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAgcg                                        >  2:2142834/1‑100 (MQ=255)
                                                        ttAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGcc                                      <  2:702998/100‑1 (MQ=255)
                                                        ttAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGcc                                      <  1:3670274/100‑1 (MQ=255)
                                                        ctAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGcc                                      <  2:946116/99‑1 (MQ=255)
                                                         tAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCt                                     <  1:1168698/100‑1 (MQ=255)
                                                          aaCAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGccag                                    >  1:966344/1‑98 (MQ=255)
                                                           aCAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAGCGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGa                                   <  1:1987277/100‑1 (MQ=255)
                                                            cAGATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGAc                                  <  2:2750428/100‑1 (MQ=255)
                                                              gATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGt                                >  1:1111306/1‑100 (MQ=255)
                                                              gATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGt                                >  2:905518/1‑100 (MQ=255)
                                                              gATAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGt                                <  1:4010680/100‑1 (MQ=255)
                                                               aTAAACTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTg                               >  2:935226/1‑100 (MQ=255)
                                                                 aaaCTCGTGATTAACGGCAATAAAGTGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTTCTGACGTTGCAACAACTGTTTTAGCTCCTGACTTGaa                             >  2:39994/1‑100 (MQ=255)
                                                                  aaCTCTTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGaaa                            >  2:338878/1‑100 (MQ=255)
                                                                    ctctTAATAAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGACGGCTAATTATTGCCCCCGTTGCACTCACCCCTTTACCCCCACACGTGAAAtc                          >  2:453488/1‑99 (MQ=255)
                                                                     tctTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAAttt                         >  1:854921/1‑100 (MQ=255)
                                                                     tctTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAAttt                         >  2:2165775/1‑100 (MQ=255)
                                                                      ctTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAAtttt                        <  1:1063481/100‑1 (MQ=255)
                                                                      ctTTATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAAtttt                        >  1:3387984/1‑100 (MQ=255)
                                                                        ttATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTaa                      <  2:1545687/100‑1 (MQ=255)
                                                                         tATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAAt                     <  1:3520704/100‑1 (MQ=255)
                                                                         tATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAAt                     >  1:316671/1‑100 (MQ=255)
                                                                         tATCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTACTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAAt                     <  1:4010924/100‑1 (MQ=255)
                                                                          aTCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATa                    <  2:1612901/100‑1 (MQ=255)
                                                                          aTCAACGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATa                    >  2:1241812/1‑100 (MQ=255)
                                                                            cAACGGCAATGGACTGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACCCTTTTAGCGCCTGCCGTGAAATTTTGATAcc                  >  2:2371202/1‑100 (MQ=255)
                                                                             aaCGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATAccg                 >  2:2778892/1‑100 (MQ=255)
                                                                             aaCGGCAATGGAACGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGTCGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATAccg                 <  1:2821499/100‑1 (MQ=255)
                                                                               cGGCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATAccgcc               >  2:1531759/1‑100 (MQ=255)
                                                                                 gCAATGGAATGGCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGt             >  1:2157679/1‑100 (MQ=255)
                                                                                    atggaatggCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGtttt          >  2:694094/1‑100 (MQ=255)
                                                                                    atggaatcgCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAAGACCGCCGttct          <  2:3705810/100‑3 (MQ=255)
                                                                                     tggaatggCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGttttt         >  1:2166659/1‑100 (MQ=255)
                                                                                     tggaatggCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGttttt         >  1:1957219/1‑100 (MQ=255)
                                                                                       gaatggCAGCGATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCa       <  1:2599764/100‑1 (MQ=255)
                                                                                       gaatggCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCa       >  1:1898461/1‑100 (MQ=255)
                                                                                       gaatggCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCa       <  2:2749044/100‑1 (MQ=255)
                                                                                          tggCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCt    <  1:921638/100‑1 (MQ=255)
                                                                                           ggCAGCAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTc   >  1:2056325/1‑100 (MQ=255)
                                                                                            gcagcaATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCTGCTcc  >  2:1229146/1‑100 (MQ=255)
                                                                                            gcagcaATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTcc  <  1:2525971/100‑1 (MQ=255)
                                                                                                |                                                                                               
CCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCC  >  NC_000913/2373526‑2373717

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