Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,847,265 A→G W167R (TGG→CGG)  b3668 ← sensor histidine protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,847,2650AG99.6% 796.2 / NA 253W167R (TGG→CGG) b3668sensor histidine protein kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (112/140);  total (113/140)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 272 total aligned reads displayed.

GACTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCAGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCACG  >  NC_000913/3847175‑3847354
                                                                                          |                                                                                         
gaCTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCggg                                                                                         >  2:5229895/1‑93 (MQ=255)
gaCTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCggg                                                                                         <  2:3037200/93‑1 (MQ=255)
    gATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGtt                                                                                     <  2:2981405/93‑1 (MQ=255)
      tAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACTGACCGAGCGGCAGCCGGGTGTTAt                                                                                   >  2:6180493/1‑93 (MQ=255)
      tAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTAt                                                                                   <  1:2214177/93‑1 (MQ=255)
       aacaaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAACCCAGTGACGGACCGAGCGGCAGCCGGGTGTTAtt                                                                                  <  2:3950026/93‑1 (MQ=255)
       aacaaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTAtt                                                                                  <  2:4802856/93‑1 (MQ=255)
        acaaacaGCAGCAAGTACCAGACCAGCTGTCGCCCGCGCTAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTg                                                                                 <  1:3361998/93‑1 (MQ=255)
        acaaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGGCGGGTGGTATTg                                                                                 >  2:3780067/1‑93 (MQ=255)
         caaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTgg                                                                                <  2:671964/93‑1 (MQ=255)
         caaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTgg                                                                                >  2:1747229/1‑93 (MQ=255)
          aaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTggg                                                                               >  2:650220/1‑92 (MQ=255)
            acaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCga                                                                             <  2:21060/93‑1 (MQ=255)
             cagcagCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCgag                                                                            >  1:1020749/1‑93 (MQ=255)
               gcagcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGAt                                                                          <  1:2685064/93‑1 (MQ=255)
                cagcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATa                                                                         <  1:6028555/93‑1 (MQ=255)
                 agcaAGTACCAGACCAGATGTCTCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAg                                                                        <  2:713151/93‑1 (MQ=255)
                 agcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAg                                                                        <  2:5665344/93‑1 (MQ=255)
                  gcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGt                                                                       >  2:1985485/1‑93 (MQ=255)
                    aaGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGCGCGGCAGCCGGGTGTTATTGGCGCGATAGTGc                                                                     >  1:1916940/1‑93 (MQ=255)
                      gTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCa                                                                   <  2:1725004/93‑1 (MQ=255)
                      gTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCa                                                                   >  1:800466/1‑93 (MQ=255)
                          cagaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAac                                                               <  2:1240493/93‑1 (MQ=255)
                              ccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGGCAGAACACCa                                                           >  2:5787984/1‑93 (MQ=255)
                              ccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCa                                                           >  1:4625117/1‑93 (MQ=255)
                               cagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAg                                                          >  1:3885212/1‑93 (MQ=255)
                                agaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGa                                                         >  1:1098317/1‑93 (MQ=255)
                                 gaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGAc                                                        >  2:1882492/1‑93 (MQ=255)
                                  aTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACa                                                       <  2:4513102/93‑1 (MQ=255)
                                   tGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGCGGTATTGGCGAGATAGTGCCAGAACACCAGACat                                                      >  1:1275590/1‑93 (MQ=255)
                                    gTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGGGTTATTGGCGAGATAGTGTCAGAACACCAGACata                                                     >  2:4527846/1‑93 (MQ=255)
                                     tCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACatat                                                    <  1:3259147/93‑1 (MQ=255)
                                     gcGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACatat                                                    <  1:540393/92‑1 (MQ=255)
                                      cGCCCGCGCCAGTTGATTGGCTGAGAAGCCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATc                                                   <  2:4086593/93‑1 (MQ=255)
                                       gCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCg                                                  <  2:4880262/93‑1 (MQ=255)
                                         ccGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGGCAGAACACCCGACATATCggg                                                >  1:4488008/1‑93 (MQ=255)
                                         ccGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCggg                                                >  1:5514722/1‑93 (MQ=255)
                                         ccGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGCGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCggg                                                >  2:2474065/1‑93 (MQ=255)
                                          cgcgCCAGTTGATTGTCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCgggg                                               <  1:775483/93‑1 (MQ=255)
                                          cgcgCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCgggg                                               >  2:931130/1‑93 (MQ=255)
                                           gcgcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGc                                              >  1:1866694/1‑93 (MQ=255)
                                           gcgcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGc                                              >  1:3030073/1‑93 (MQ=255)
                                            cgcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGcc                                             >  2:5001192/1‑93 (MQ=255)
                                             gcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCa                                            >  1:3315221/1‑93 (MQ=255)
                                             gcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCa                                            <  1:4668501/93‑1 (MQ=255)
                                             gcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCa                                            >  1:3479482/1‑93 (MQ=255)
                                             gcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCa                                            <  1:1252151/93‑1 (MQ=255)
                                               cAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGc                                          <  1:2863417/93‑1 (MQ=255)
                                               cAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGc                                          <  1:216117/93‑1 (MQ=255)
                                                aGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCCGACATATCGGGGCCAGCg                                         >  1:5278164/1‑93 (MQ=255)
                                                aGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCg                                         >  2:1163273/1‑93 (MQ=255)
                                                aGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCg                                         <  2:3430371/93‑1 (MQ=255)
                                                 gTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGt                                        <  1:3970044/93‑1 (MQ=255)
                                                 gTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGt                                        <  1:1156694/93‑1 (MQ=255)
                                                  ttGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTc                                       >  1:5606406/1‑93 (MQ=255)
                                                  ttGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTc                                       >  1:816993/1‑93 (MQ=255)
                                                  ttGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTc                                       >  1:3851041/1‑93 (MQ=255)
                                                  ttGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTc                                       <  2:3352381/93‑1 (MQ=255)
                                                   tGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCa                                      <  1:4075067/93‑1 (MQ=255)
                                                   tGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCa                                      <  2:1629905/93‑1 (MQ=255)
                                                   tGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCa                                      >  2:1497182/1‑93 (MQ=255)
                                                    gATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAg                                     <  2:45434/93‑1 (MQ=255)
                                                     aTTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAgg                                    <  1:4422645/93‑1 (MQ=255)
                                                     aTTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAgg                                    >  1:6149778/1‑93 (MQ=255)
                                                     aTTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAgg                                    <  2:4233506/93‑1 (MQ=255)
                                                        ggCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCGGCGTCAGgccg                                 <  2:3281885/93‑1 (MQ=255)
                                                        ggCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGgccg                                 >  1:1735050/1‑93 (MQ=255)
                                                         gCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGgccgc                                <  2:6249264/93‑1 (MQ=255)
                                                          cTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGgccgcc                               <  1:2537586/93‑1 (MQ=255)
                                                           tGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCa                              >  2:4715816/1‑93 (MQ=255)
                                                           tGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCa                              >  1:2065231/1‑93 (MQ=255)
                                                            gagaAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAg                             >  1:5111040/1‑93 (MQ=255)
                                                             agaAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGt                            >  2:1314405/1‑93 (MQ=255)
                                                              gaAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGtt                           <  2:380271/93‑1 (MQ=255)
                                                                  ccAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAg                       <  1:5217460/93‑1 (MQ=255)
                                                                   cAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGt                      >  2:5389583/1‑93 (MQ=255)
                                                                    aGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATTGGGGCCAGCGTCAGGCCGCCAGGTAAAGTc                     >  2:3835895/1‑93 (MQ=255)
                                                                      tGACGGACCGAGCGTCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAg                   <  2:5842406/93‑1 (MQ=255)
                                                                      tGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAg                   <  1:2371764/93‑1 (MQ=255)
                                                                        aCGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCa                 >  1:1654388/1‑93 (MQ=255)
                                                                        aCGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCa                 <  1:6047225/93‑1 (MQ=255)
                                                                         cGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCaa                <  2:2127649/93‑1 (MQ=255)
                                                                          ggACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAAc               >  1:3707356/1‑93 (MQ=255)
                                                                          ggACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAAc               <  1:67713/93‑1 (MQ=255)
                                                                           gACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACa              >  2:4603770/1‑93 (MQ=255)
                                                                            aCCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAg             <  2:5466064/93‑1 (MQ=255)
                                                                            aCCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAg             <  2:6058379/93‑1 (MQ=255)
                                                                            aCCGAGCGGCAGCCGGGGGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAg             >  1:589360/1‑93 (MQ=255)
                                                                               gAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAgcgc          <  2:5113101/93‑1 (MQ=255)
                                                                                 gcggcgGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCAt        <  2:1135944/93‑1 (MQ=255)
                                                                                 gCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCAt        <  2:2741821/93‑1 (MQ=255)
                                                                                 gCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCAt        >  1:4209332/1‑93 (MQ=255)
                                                                                 gCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCAt        <  2:1551066/93‑1 (MQ=255)
                                                                                    gCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTcc     <  2:5518724/93‑1 (MQ=255)
                                                                                     cAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCa    <  2:6117627/93‑1 (MQ=255)
                                                                                     cAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCa    <  1:993293/93‑1 (MQ=255)
                                                                                     cAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCa    >  2:840509/1‑93 (MQ=255)
                                                                                     aaGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCa    <  2:5557046/92‑1 (MQ=255)
                                                                                      aGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCAc   >  2:3542712/1‑93 (MQ=255)
                                                                                       gCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCACg  <  2:4113900/93‑1 (MQ=255)
                                                                                          |                                                                                         
GACTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCAGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCACG  >  NC_000913/3847175‑3847354

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