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.5% 659.8 / ‑4.9 214W167R (TGG→CGG) b3668sensor histidine protein kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (82/131);  total (83/131)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 239 total aligned reads displayed.

GACTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCAGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCACG  >  NC_000913/3847175‑3847354
                                                                                          |                                                                                         
 aCTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgt                                                                                        <  1:676988/93‑1 (MQ=255)
  cTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCgggat                                                                                       >  1:4038020/1‑91 (MQ=255)
  cTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgtg                                                                                       >  1:4227803/1‑93 (MQ=255)
   tGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgtgt                                                                                      >  2:3231272/1‑93 (MQ=255)
   tGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgtgt                                                                                      <  2:2510818/93‑1 (MQ=255)
   tGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgtgt                                                                                      <  1:2184968/93‑1 (MQ=255)
   tGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgtgt                                                                                      >  2:2147126/1‑93 (MQ=255)
   tGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgtgt                                                                                      >  1:1818402/1‑93 (MQ=255)
   tGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgtgt                                                                                      <  1:1762794/93‑1 (MQ=255)
   tGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGACAAACCAGTGACGGACCGAGCGGCAGCCGGgggt                                                                                      >  1:2807180/1‑93 (MQ=255)
     aTAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGTCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTa                                                                                    <  1:812149/93‑1 (MQ=255)
     aTAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTa                                                                                    <  2:3137593/93‑1 (MQ=255)
     aTAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTa                                                                                    <  1:1188555/93‑1 (MQ=255)
      tatcaaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTAt                                                                                   >  1:3780903/4‑93 (MQ=255)
       aacaaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACTGAGCGGCAGCCGGGCGTTAtt                                                                                  >  1:4590602/1‑93 (MQ=255)
       aacaaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTAtt                                                                                  <  1:3741054/93‑1 (MQ=255)
       aacaaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTAtt                                                                                  <  2:1311338/93‑1 (MQ=255)
        acaaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTg                                                                                 <  2:3681452/93‑1 (MQ=255)
        acaaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTg                                                                                 <  1:2999533/93‑1 (MQ=255)
         caaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTgg                                                                                <  2:2755709/93‑1 (MQ=255)
           acgtgcagcaAGTGCCAGACCAGATGCCGCGCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCg                                                                              <  1:3338549/89‑1 (MQ=255)
           aacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCg                                                                              >  2:2373538/1‑93 (MQ=255)
           aacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCg                                                                              >  1:1210556/1‑93 (MQ=255)
           aacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGACGGGTGTTATTGGCg                                                                              >  1:10819/1‑93 (MQ=255)
           aacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGAAGCCGGGTGTTATTGTCg                                                                              <  2:1460097/93‑1 (MQ=255)
              agcagcAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCgaga                                                                           <  2:592637/93‑1 (MQ=255)
               gcagcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGAt                                                                          <  2:2745837/93‑1 (MQ=255)
               gcagcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGAt                                                                          >  1:1655214/1‑93 (MQ=255)
                 agcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAg                                                                        <  1:3271568/93‑1 (MQ=255)
                 agcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAg                                                                        >  2:1050702/1‑93 (MQ=255)
                 agcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAg                                                                        <  2:3738825/93‑1 (MQ=255)
                  gcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAgg                                                                       >  1:2146019/1‑92 (MQ=255)
                  gcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGt                                                                       <  2:1825597/93‑1 (MQ=255)
                  gcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGt                                                                       <  2:1973719/93‑1 (MQ=255)
                   caAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTg                                                                      >  2:449388/1‑93 (MQ=255)
                    aaGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGGGTTATTGGCGAGATAGGGc                                                                     >  1:3976505/1‑93 (MQ=255)
                     aGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGcc                                                                    >  2:1502875/1‑93 (MQ=255)
                      gTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCa                                                                   >  1:2538677/1‑93 (MQ=255)
                       tACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAg                                                                  >  1:3197394/1‑93 (MQ=255)
                        accagaccagATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGa                                                                 >  2:392979/1‑93 (MQ=255)
                          cagaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAac                                                               <  1:870521/93‑1 (MQ=255)
                           agaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAaca                                                              >  2:2234919/1‑93 (MQ=255)
                           agaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAaca                                                              <  1:1854288/93‑1 (MQ=255)
                           agaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAaca                                                              <  1:4384282/93‑1 (MQ=255)
                           agaccagaTGTCGCCCGCGCCAGCTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGGGTTATTGGCGATGAAGTGCCCGAaca                                                              >  2:3582748/1‑93 (MQ=255)
                            gaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGGCAGAacac                                                             >  1:849060/1‑93 (MQ=255)
                            gaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAacac                                                             >  1:535794/1‑93 (MQ=255)
                             accagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACAcc                                                            >  1:1951725/1‑93 (MQ=255)
                               cagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAg                                                          <  1:2578930/93‑1 (MQ=255)
                               cagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGGGACGGACCGAGGGGCAGCCGGGTGTTATTTGAGAGATAGTTCCAGAACACCaa                                                          >  2:20648/1‑92 (MQ=255)
                                cgacgTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGa                                                         <  1:14959/89‑1 (MQ=255)
                                 gaTGTCGCCCGCGCCAGTTGATTGGCTGGGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGAc                                                        <  2:3789755/93‑1 (MQ=255)
                                 gaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGAc                                                        >  1:735057/1‑93 (MQ=255)
                                 gaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGAc                                                        <  1:4351635/93‑1 (MQ=255)
                                 gaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGAc                                                        >  2:3330352/1‑93 (MQ=255)
                                  aTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGCGTTATTGGCGAGATAGTGCCAGAACACCAGACa                                                       >  2:3628965/1‑93 (MQ=255)
                                  aTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGAGTTATTGGCGAGATAGTGCCAGAACACCAGACa                                                       >  2:767640/1‑93 (MQ=255)
                                    gTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACata                                                     <  1:3723873/93‑1 (MQ=255)
                                       gCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGCGCGGCAGCCGGGGGGTATTTGCGCGATAGTGCCAGAACACCAGACATATCg                                                  >  2:2312437/1‑93 (MQ=255)
                                       gCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCg                                                  <  1:3085769/93‑1 (MQ=255)
                                       gCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCg                                                  <  2:4193754/93‑1 (MQ=255)
                                         ccGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCggg                                                <  2:4553727/93‑1 (MQ=255)
                                         ccGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCggg                                                >  1:2684431/1‑93 (MQ=255)
                                         ccGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCggg                                                <  1:3588757/93‑1 (MQ=255)
                                         ccGCGCAAGTTGATTGGCTGAGAAACCAGTGACGTACCGAGCGTCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACGCCAGACATATCggg                                                <  2:2994987/93‑1 (MQ=255)
                                             gcCAGTTGATTGGCTGAGAAGCCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCa                                            <  2:42921/93‑1 (MQ=255)
                                              ccAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAg                                           <  2:2282171/93‑1 (MQ=255)
                                              ccAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAg                                           >  1:2886541/1‑93 (MQ=255)
                                                aGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACAAATCGGGGCCGGCg                                         >  2:2597393/1‑93 (MQ=255)
                                                  ttGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGTCAGCGTc                                       >  1:2042625/1‑93 (MQ=255)
                                                    gATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAg                                     >  1:2455342/1‑93 (MQ=255)
                                                    gATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAg                                     <  2:4462025/93‑1 (MQ=255)
                                                    gATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGGACCCAAACCCCAGACATATCGGGGCCAGCGTCAg                                     >  1:3522724/1‑93 (MQ=255)
                                                     aTTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAgg                                    >  1:2883249/1‑93 (MQ=255)
                                                      ttGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGgc                                   >  1:3846126/1‑93 (MQ=255)
                                                      ttGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGgc                                   >  2:2869624/1‑93 (MQ=255)
                                                       tGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGGGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCATgcc                                  >  1:4613922/1‑93 (MQ=255)
                                                          cTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGgccgcc                               <  2:2339403/93‑1 (MQ=255)
                                                           tGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACTCCAGACATATCGGGGCCAGCGTCAGGCCGCCa                              >  2:4362070/1‑93 (MQ=255)
                                                            gagaAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAg                             >  1:1061140/1‑93 (MQ=255)
                                                             agaAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGt                            <  2:2153954/93‑1 (MQ=255)
                                                               aaaCCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTa                          <  1:2697110/93‑1 (MQ=255)
                                                                 aCCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTaaa                        <  1:804505/93‑1 (MQ=255)
                                                                  ccAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAg                       <  2:1665525/93‑1 (MQ=255)
                                                                  ccAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAg                       <  2:1872772/93‑1 (MQ=255)
                                                                  ccAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAg                       >  2:1368319/1‑93 (MQ=255)
                                                                   cAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGt                      <  1:1294281/93‑1 (MQ=255)
                                                                    aGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTc                     <  2:3559090/93‑1 (MQ=255)
                                                                    aGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTc                     <  1:1029588/93‑1 (MQ=255)
                                                                     gTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCa                    >  1:4106511/1‑93 (MQ=255)
                                                                     gTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCa                    <  1:999232/93‑1 (MQ=255)
                                                                           gACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACa              <  1:3892144/93‑1 (MQ=255)
                                                                            aCCGAGCGCCACCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAg             <  2:2565501/93‑1 (MQ=255)
                                                                              cGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCCGACATATCGGGGTCAGCGTCAGGCCGCCCGTTAAAGCAAGCAACAgcg           >  2:3796998/1‑93 (MQ=255)
                                                                               gcgCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAgcgc          >  2:2835988/3‑93 (MQ=255)
                                                                               gAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAgcgc          >  1:2324363/1‑93 (MQ=255)
                                                                                 gCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCAt        >  2:2655332/1‑93 (MQ=255)
                                                                                    gCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTcc     <  2:2591155/93‑1 (MQ=255)
                                                                                    gCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTcc     <  2:4352444/93‑1 (MQ=255)
                                                                                      aGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCAc   <  2:161180/93‑1 (MQ=255)
                                                                                          |                                                                                         
GACTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCAGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCACG  >  NC_000913/3847175‑3847354

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