Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,810,132 C→T R229C (CGC→TGC)  emrA → EmrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,810,1320CT86.3% 710.6 / 68.8 263R229C (CGC→TGC) emrAEmrA
Reads supporting (aligned to +/- strand):  ref base C (16/20);  new base T (120/107);  total (136/127)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-02

Only 100 of 286 total aligned reads displayed.

CCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAA  >  NC_000913/2810040‑2810223
                                                                                             |                                                                                           
  gTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCcgctgc                                                                                           >  1:7452444/1‑93 (MQ=255)
  gTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCcgctgc                                                                                           <  1:5259442/93‑1 (MQ=255)
  gTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCcgccgc                                                                                           >  2:1475564/1‑93 (MQ=255)
  gTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCcgccgc                                                                                           >  1:831890/1‑93 (MQ=255)
  gTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCcgccgc                                                                                           >  2:6001145/1‑93 (MQ=255)
    gCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCtgcgc                                                                                         <  1:6103043/93‑1 (MQ=255)
    gCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCtgcgc                                                                                         >  2:3116664/1‑93 (MQ=255)
      aaCAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCgg                                                                                       >  1:2737118/1‑93 (MQ=255)
      aaCAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCgg                                                                                       >  1:4073525/1‑93 (MQ=255)
      aaCAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCgg                                                                                       <  1:5702136/93‑1 (MQ=255)
      aaCAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCgg                                                                                       >  1:189788/1‑93 (MQ=255)
       aCAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGt                                                                                      >  1:7066289/1‑93 (MQ=255)
       aCAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGt                                                                                      >  2:4413250/1‑93 (MQ=255)
        cAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTa                                                                                     >  1:2538349/1‑93 (MQ=255)
         aGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTAc                                                                                    <  2:3793932/93‑1 (MQ=255)
          ggCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACa                                                                                   >  1:4894277/1‑93 (MQ=255)
          ggCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACa                                                                                   <  2:5275679/93‑1 (MQ=255)
           gCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAg                                                                                  <  1:5023039/93‑1 (MQ=255)
           gCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAg                                                                                  >  2:6923401/1‑93 (MQ=255)
             tGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTTCTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGcc                                                                                <  1:6144918/93‑1 (MQ=255)
             tGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGcc                                                                                >  1:584135/1‑93 (MQ=255)
               ccACCGAAGTACGTAACGCCTGGCTGGCGCTTGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCCCTGCGCGGTACAGCCTg                                                                              >  1:7420545/1‑93 (MQ=255)
               ccACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTg                                                                              <  2:1600173/93‑1 (MQ=255)
               ccACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTg                                                                              <  2:6976685/93‑1 (MQ=255)
                 aCCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTggg                                                                            <  1:4481055/93‑1 (MQ=255)
                 aCCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTggg                                                                            <  1:4716212/93‑1 (MQ=255)
                 aCCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTggg                                                                            >  1:2560167/1‑93 (MQ=255)
                 aCCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGCACAGCCTggg                                                                            >  1:4362046/1‑93 (MQ=255)
                   cGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGgc                                                                          <  2:2794902/93‑1 (MQ=255)
                     aaGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGgcgc                                                                        >  1:2837329/1‑93 (MQ=255)
                      aGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCa                                                                       >  2:4763456/1‑93 (MQ=255)
                        tACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCaaa                                                                     <  2:1908441/93‑1 (MQ=255)
                           gTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTa                                                                  >  2:784946/1‑93 (MQ=255)
                            tAACGCCTGGGTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAg                                                                 <  1:1018390/93‑1 (MQ=255)
                            tAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAg                                                                 >  2:1068370/1‑93 (MQ=255)
                            tAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAg                                                                 >  2:5369545/1‑93 (MQ=255)
                             aaCGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGc                                                                >  2:1193457/1‑93 (MQ=255)
                              aCGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGcc                                                               >  2:2199982/1‑93 (MQ=255)
                               cGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGccc                                                              >  1:7410152/1‑93 (MQ=255)
                                gCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCa                                                             >  1:2172367/1‑93 (MQ=255)
                                gCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCa                                                             <  2:5664809/93‑1 (MQ=255)
                                 ccTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCaa                                                            <  2:7458347/93‑1 (MQ=255)
                                 ccTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCaa                                                            >  2:5107688/1‑93 (MQ=255)
                                  ctggctggCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAac                                                           <  2:2318982/93‑1 (MQ=255)
                                  ctggctggCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAac                                                           <  1:5384175/93‑1 (MQ=255)
                                  ctggctggCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAac                                                           >  2:294626/1‑93 (MQ=255)
                                   tggctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAacg                                                          >  1:205313/1‑93 (MQ=255)
                                   tggctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAacg                                                          >  1:6497577/1‑93 (MQ=255)
                                    ggctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAacga                                                         <  2:1775655/93‑1 (MQ=255)
                                      ctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAacgacg                                                       >  1:4253365/1‑93 (MQ=255)
                                         gcgcTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGAcgccg                                                    >  2:350958/1‑93 (MQ=255)
                                            cTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTg                                                 >  1:7204479/1‑93 (MQ=255)
                                             tGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGa                                                >  1:3982116/1‑93 (MQ=255)
                                               gAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATg                                              >  1:3332321/1‑93 (MQ=255)
                                               gAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATg                                              >  2:573576/1‑93 (MQ=255)
                                                aGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATgg                                             <  2:3313983/93‑1 (MQ=255)
                                                 gCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGc                                            <  1:6592500/93‑1 (MQ=255)
                                                 gCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGc                                            >  2:5870394/1‑93 (MQ=255)
                                                 gCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGc                                            >  2:7156602/1‑93 (MQ=255)
                                                  cGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCg                                           <  1:5066165/93‑1 (MQ=255)
                                                   gTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCgg                                          >  2:3152206/1‑93 (MQ=255)
                                                   gTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCgg                                          <  2:6510351/93‑1 (MQ=255)
                                                    tACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGt                                         <  1:4884273/93‑1 (MQ=255)
                                                    tACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGt                                         >  1:684789/1‑93 (MQ=255)
                                                    tACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGt                                         <  2:6290495/93‑1 (MQ=255)
                                                     aCTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTc                                        >  2:2061227/1‑93 (MQ=255)
                                                          tattatCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTcc                                   <  2:1409478/93‑1 (MQ=255)
                                                           attatCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCa                                  >  2:2043875/1‑93 (MQ=255)
                                                           attatCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCa                                  >  2:859307/1‑93 (MQ=255)
                                                           attatCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCa                                  <  1:5284944/93‑1 (MQ=255)
                                                            ttatCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAg                                 >  1:496007/1‑93 (MQ=255)
                                                             tatCAGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGc                                >  2:6605910/1‑93 (MQ=255)
                                                                 aGTCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGccacc                            <  2:18431/93‑1 (MQ=255)
                                                                   tCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCaa                          <  1:2627121/93‑1 (MQ=255)
                                                                   tCCGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCaa                          <  2:6830850/93‑1 (MQ=255)
                                                                    ccGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAat                         <  1:6410499/93‑1 (MQ=255)
                                                                     cGATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAata                        >  1:6927154/1‑93 (MQ=255)
                                                                      gATGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAatat                       >  1:2509656/1‑93 (MQ=255)
                                                                        tGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATAtgt                     <  2:414611/93‑1 (MQ=255)
                                                                        tGACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATAtgt                     >  1:170838/1‑93 (MQ=255)
                                                                         gACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATAtgtg                    <  1:2908218/93‑1 (MQ=255)
                                                                         gACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATAtgtg                    <  1:3724943/93‑1 (MQ=255)
                                                                         gACCGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATAtgtg                    <  2:4536469/93‑1 (MQ=255)
                                                                          aCCGGTTATGTCTCGCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATAtgtg                    <  1:1189802/93‑1 (MQ=255)
                                                                           ccGGTTATGTCTCCCGCTGCGCGGGACAGCCTGGGGAGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATgtggg                  >  1:5185524/1‑93 (MQ=255)
                                                                           ccGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGGCGCCGCTGATGGCGGTCGTTCCAGCCACCAATATgtggg                  <  1:6214533/93‑1 (MQ=255)
                                                                            cGGTTATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATgtgggt                 >  1:3520871/1‑93 (MQ=255)
                                                                             ggTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATgtgggtg                >  1:3132716/1‑93 (MQ=255)
                                                                              gTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATgtgggtgg               >  2:985645/1‑93 (MQ=255)
                                                                               ttATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGa              >  2:3001919/1‑93 (MQ=255)
                                                                                tATGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGAt             <  1:3594710/93‑1 (MQ=255)
                                                                                  tGTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGc           >  2:7371264/1‑93 (MQ=255)
                                                                                  tGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGc           >  1:1372796/1‑93 (MQ=255)
                                                                                   gTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGcc          >  2:4946377/1‑93 (MQ=255)
                                                                                   gTCTCCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGcc          <  1:5885892/93‑1 (MQ=255)
                                                                                      tcCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAAc       >  2:6421740/1‑93 (MQ=255)
                                                                                      tcCCGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAAc       >  1:3112868/1‑93 (MQ=255)
                                                                                       cccGCTGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACt      <  1:5508806/93‑1 (MQ=255)
                                                                                          cgctgcGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACttt    <  2:5549535/93‑1 (MQ=255)
                                                                                          cgctgcGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACttt    >  1:412271/1‑93 (MQ=255)
                                                                                             |                                                                                           
CCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAA  >  NC_000913/2810040‑2810223

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