Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 348,519 T→C V205A (GTG→GCG)  prpB → putative phosphonomutase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913348,5190TC84.0% 507.6 / 85.0 214V205A (GTG→GCG) prpBputative phosphonomutase 2
Reads supporting (aligned to +/- strand):  ref base T (20/14);  new base C (92/87);  total (113/101)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

Only 100 of 233 total aligned reads displayed.

GCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTA  >  NC_000913/348427‑348609
                                                                                            |                                                                                          
  aGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGgcgc                                                                                          <  2:5920285/93‑1 (MQ=255)
  aGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGgcgc                                                                                          >  2:4255835/1‑93 (MQ=255)
  aGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGgcgc                                                                                          <  1:83122/93‑1 (MQ=255)
  aGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGgcgc                                                                                          <  1:3165275/93‑1 (MQ=255)
    gCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCg                                                                                        <  2:4499237/93‑1 (MQ=255)
    gCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCg                                                                                        >  1:1912715/1‑93 (MQ=255)
     ccTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGa                                                                                       <  1:1207175/93‑1 (MQ=255)
     ccTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGa                                                                                       >  1:5169999/1‑93 (MQ=255)
     ccTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGa                                                                                       >  1:2337065/1‑93 (MQ=255)
         tGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCt                                                                                   >  2:6186015/1‑93 (MQ=255)
         tGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCt                                                                                   >  2:241557/1‑93 (MQ=255)
         tGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCt                                                                                   <  1:370499/93‑1 (MQ=255)
          gTTGAAGTGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTc                                                                                  <  2:5872625/93‑1 (MQ=255)
           ttGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCg                                                                                 <  2:1221445/93‑1 (MQ=255)
             gAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGcc                                                                               >  1:4241221/1‑93 (MQ=255)
             gAAGCGGGTGCCGAGATGTTGTTCCCGGAGGAGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGcc                                                                               >  1:4806052/1‑93 (MQ=255)
               aGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCaa                                                                             >  2:2867105/1‑93 (MQ=255)
               aGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCaa                                                                             >  2:1104803/1‑93 (MQ=255)
                gCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGACCTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAAc                                                                            >  2:3328319/1‑93 (MQ=255)
                   ggTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATc                                                                         <  1:3927792/93‑1 (MQ=255)
                   ggTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATc                                                                         >  1:4367908/1‑93 (MQ=255)
                   ggTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATc                                                                         >  1:4432122/1‑93 (MQ=255)
                    gTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCa                                                                        <  1:5878134/93‑1 (MQ=255)
                     tGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCGGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCAc                                                                       <  1:4903628/93‑1 (MQ=255)
                      gCCTAGATGTTTTTCCCGGAGGCGATTACCGACCTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCAcc                                                                      <  1:6422542/93‑1 (MQ=255)
                        cGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGa                                                                    >  1:5623483/1‑93 (MQ=255)
                         gagaTGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGaa                                                                   >  2:1702988/1‑93 (MQ=255)
                          agaTGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAAt                                                                  <  2:3998678/93‑1 (MQ=255)
                          agaTGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAAt                                                                  >  1:1335322/1‑93 (MQ=255)
                          agaTGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAAt                                                                  >  1:5379041/1‑93 (MQ=255)
                           gaTGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAAtt                                                                 >  2:4210345/1‑93 (MQ=255)
                              gttgttCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTgg                                                              <  1:401281/93‑1 (MQ=255)
                              gttgttCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTgg                                                              >  1:2792956/1‑93 (MQ=255)
                               ttgttCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGt                                                             >  1:6317545/1‑93 (MQ=255)
                               ttgttCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGt                                                             <  2:3798127/93‑1 (MQ=255)
                                tgttCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTg                                                            >  1:5522684/1‑93 (MQ=255)
                                  ttCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGcc                                                          >  2:4128232/1‑93 (MQ=255)
                                  ttCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGcc                                                          <  1:2911358/93‑1 (MQ=255)
                                   tCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCa                                                         >  2:3654622/1‑93 (MQ=255)
                                    cccGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGGCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAc                                                        >  1:2089239/1‑93 (MQ=255)
                                      cGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAcgc                                                      >  2:6134200/1‑93 (MQ=255)
                                      cGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAcgc                                                      >  1:1993968/1‑93 (MQ=255)
                                       ggAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAcgcc                                                     >  1:4959664/1‑93 (MQ=255)
                                         aGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAcgccgc                                                   <  2:236878/93‑1 (MQ=255)
                                         aGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAcgccgc                                                   <  1:2625674/93‑1 (MQ=255)
                                         aGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAcgccgc                                                   >  2:4657688/1‑93 (MQ=255)
                                         aGGCGATGCCCGAACTCGCCATGTTCCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCAcgccgc                                                   <  1:168484/93‑1 (MQ=255)
                                          ggCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCt                                                  >  2:2023401/1‑93 (MQ=255)
                                          ggCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCt                                                  <  1:3770431/93‑1 (MQ=255)
                                            cGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGt                                                <  1:1893274/93‑1 (MQ=255)
                                             gATTACCGAACTCGCCATGTCCCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGtt                                               <  1:2172944/93‑1 (MQ=255)
                                             gATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGtt                                               <  2:1607171/93‑1 (MQ=255)
                                             gATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGtt                                               <  2:2116313/93‑1 (MQ=255)
                                              aTTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGttt                                              >  1:449239/1‑93 (MQ=255)
                                               ttACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTa                                             >  1:1060515/1‑93 (MQ=255)
                                               ttACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTa                                             >  1:4579707/1‑93 (MQ=255)
                                                 aCCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTacc                                           >  1:5738366/1‑93 (MQ=255)
                                                  ccGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTacca                                          <  2:1889621/93‑1 (MQ=255)
                                                   cGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTaccac                                         <  1:4293118/93‑1 (MQ=255)
                                                    gAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTaccacc                                        <  2:1319332/93‑1 (MQ=255)
                                                      aCTCGCCATGTACTGCCAGTTTGCCGATGCGTTGCAGGCGCTGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACcga                                      <  1:2972311/93‑1 (MQ=255)
                                                      aCTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACcga                                      <  2:4815712/93‑1 (MQ=255)
                                                        tCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACcgacg                                    >  2:5283292/1‑93 (MQ=255)
                                                        tCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACcgacg                                    <  1:6521532/93‑1 (MQ=255)
                                                         cGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACcgacga                                   >  1:6486575/1‑93 (MQ=255)
                                                          gCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGaa                                  >  2:2262141/1‑93 (MQ=255)
                                                          gCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGaa                                  >  1:3542244/1‑93 (MQ=255)
                                                           ccATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAAt                                 >  1:404393/1‑93 (MQ=255)
                                                            cATGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAAtt                                <  2:4765958/93‑1 (MQ=255)
                                                             aTGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTa                               >  1:2589301/1‑93 (MQ=255)
                                                             aTGTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTa                               <  1:3378971/93‑1 (MQ=255)
                                                              tGTACCCCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTAc                              <  2:3192129/93‑1 (MQ=255)
                                                               gTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACg                             >  2:6364676/1‑93 (MQ=255)
                                                               gTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACg                             <  2:6488899/93‑1 (MQ=255)
                                                               gTACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACg                             >  1:3473030/1‑93 (MQ=255)
                                                                tACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGc                            >  2:2679538/1‑93 (MQ=255)
                                                                tACCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGGGCCACGCCGCTGTTTACCACCGACGGATTACga                            >  1:6172381/1‑92 (MQ=255)
                                                                 aCCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCa                           <  2:3175003/93‑1 (MQ=255)
                                                                 aCCGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCa                           <  1:943546/93‑1 (MQ=255)
                                                                  ccGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAg                          <  2:5265289/93‑1 (MQ=255)
                                                                  ccGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAg                          <  2:5110582/93‑1 (MQ=255)
                                                                  ccGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAg                          <  2:324655/93‑1 (MQ=255)
                                                                  ccGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAg                          >  2:5569407/1‑93 (MQ=255)
                                                                   cGCCAGTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAgc                         <  2:6639125/93‑1 (MQ=255)
                                                                      cAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGcc                      >  2:3495291/1‑93 (MQ=255)
                                                                        gTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCa                    <  1:1634837/93‑1 (MQ=255)
                                                                        gTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCa                    <  2:3960242/93‑1 (MQ=255)
                                                                        gTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCa                    <  2:373578/93‑1 (MQ=255)
                                                                        gTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCa                    <  2:537569/93‑1 (MQ=255)
                                                                        gTTTGCCGATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCa                    <  2:2942124/93‑1 (MQ=255)
                                                                          ttGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATg                  >  2:4940791/1‑93 (MQ=255)
                                                                          ttGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATg                  <  1:3091173/93‑1 (MQ=255)
                                                                          ttGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATg                  <  1:743768/93‑1 (MQ=255)
                                                                              cGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGc              >  2:1663415/1‑93 (MQ=255)
                                                                               gATGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCa             >  1:3648805/1‑93 (MQ=255)
                                                                                 tGCGGTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAAt           <  1:2773501/93‑1 (MQ=255)
                                                                                    ggTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATgga        >  2:3065404/1‑92 (MQ=255)
                                                                                    ggTGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGgc        >  2:5025882/1‑93 (MQ=255)
                                                                                     gTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGgcg       >  1:6528688/1‑93 (MQ=255)
                                                                                      tGCAGGCGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGgcgc      <  2:6568438/93‑1 (MQ=255)
                                                                                            |                                                                                          
GCAGGCCTATGTTGAAGCGGGTGCCGAGATGTTGTTCCCGGAGGCGATTACCGAACTCGCCATGTACCGCCAGTTTGCCGATGCGGTGCAGGTGCCGATCCTCGCCAACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTA  >  NC_000913/348427‑348609

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