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.4% 543.5 / ‑4.7 170W167R (TGG→CGG) b3668sensor histidine protein kinase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (73/96);  total (73/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 179 total aligned reads displayed.

GACTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCAGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCACGA  >  NC_000913/3847175‑3847355
                                                                                          |                                                                                          
 aCTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCgggg                                                                                         >  1:1225605/1‑92 (MQ=255)
 aCTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgt                                                                                         <  1:2084656/93‑1 (MQ=255)
 aCTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgt                                                                                         <  1:2357184/93‑1 (MQ=255)
  cTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgtg                                                                                        >  1:174462/1‑93 (MQ=255)
  cTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGgtg                                                                                        <  1:4360787/93‑1 (MQ=255)
    gATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGGGtt                                                                                      >  1:3213040/1‑93 (MQ=255)
     aTAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTa                                                                                     <  1:2894493/93‑1 (MQ=255)
     aTAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCAGGTGTTa                                                                                     >  2:4406044/1‑93 (MQ=255)
      tAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTAt                                                                                    <  2:1338151/93‑1 (MQ=255)
        acaaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTg                                                                                  >  1:409053/1‑93 (MQ=255)
         caaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTgg                                                                                 >  1:723487/1‑93 (MQ=255)
         caaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTgg                                                                                 <  2:2733667/93‑1 (MQ=255)
         caaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTgg                                                                                 <  2:1946083/93‑1 (MQ=255)
          aaacaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGc                                                                                >  2:25384/1‑93 (MQ=255)
            acaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCga                                                                              <  2:3132198/93‑1 (MQ=255)
            acaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCga                                                                              <  2:2825558/93‑1 (MQ=255)
            acaGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCga                                                                              <  1:4034529/93‑1 (MQ=255)
                 agcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTTTTATTGGCGAAATAg                                                                         <  2:1275271/93‑1 (MQ=255)
                 agcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAg                                                                         >  1:2223676/1‑93 (MQ=255)
                 agcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAg                                                                         >  2:1254393/1‑93 (MQ=255)
                  gcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGt                                                                        >  1:146527/1‑93 (MQ=255)
                  gcaAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGCGTTATTGGCGAGATAGt                                                                        >  1:2502530/1‑93 (MQ=255)
                     aGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCAGGTGTTATTGGCGAGATAGTGcc                                                                     <  1:745433/93‑1 (MQ=255)
                      gTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGGGTTATTGGCGAGATAGTGCCa                                                                    >  2:3400945/1‑93 (MQ=255)
                        accagaccagATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGa                                                                  >  1:4543712/1‑93 (MQ=255)
                         ccagaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGaa                                                                 >  1:2173584/1‑93 (MQ=255)
                          cagaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAac                                                                <  2:3274032/93‑1 (MQ=255)
                          cagaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAac                                                                >  1:103934/1‑93 (MQ=255)
                           agaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAaca                                                               <  2:838115/93‑1 (MQ=255)
                           agaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAaca                                                               >  2:2580217/1‑93 (MQ=255)
                            gaccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGGGTTATTGGCGAGATAGTGCCAGAacac                                                              >  1:3607721/1‑93 (MQ=255)
                            gaccagaTGTCGCCCGCGCCAGGCGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGGGAGTGGCGAGATAGGGCCAGAacac                                                              >  1:3617868/1‑93 (MQ=255)
                              ccagaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCa                                                            >  1:2018545/1‑93 (MQ=255)
                                agaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGa                                                          >  1:1901107/1‑93 (MQ=255)
                                agaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGa                                                          <  2:380760/93‑1 (MQ=255)
                                agaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGCGTTATTGGCGAGATAGTGGCAGAACACCAGa                                                          >  1:2080316/1‑93 (MQ=255)
                                 gaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGAc                                                         >  1:2311292/1‑93 (MQ=255)
                                 gaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGAc                                                         >  1:3097808/1‑93 (MQ=255)
                                 gaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGAc                                                         >  1:1932517/1‑93 (MQ=255)
                                 gaTGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGAc                                                         <  2:1331112/93‑1 (MQ=255)
                                   tGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACat                                                       <  2:2038955/93‑1 (MQ=255)
                                   tGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACat                                                       <  1:4739489/93‑1 (MQ=255)
                                    gTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACata                                                      <  1:1569479/93‑1 (MQ=255)
                                    gTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACata                                                      <  1:1665328/93‑1 (MQ=255)
                                     tagCCCACGCCTGTAGATTGACTGAGAAGCCAGTGACGGACCGAGCGGCAGATGGGTGTTATTGGCGAGATAGTGCCAGAAAACCAGACctat                                                     <  2:1469952/91‑1 (MQ=255)
                                      cGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATc                                                    <  2:904805/93‑1 (MQ=255)
                                      cGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATc                                                    <  2:160289/93‑1 (MQ=255)
                                       gCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCg                                                   >  2:2819662/1‑93 (MQ=255)
                                       gCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCg                                                   >  2:827052/1‑93 (MQ=255)
                                       gCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCg                                                   >  2:827051/1‑93 (MQ=255)
                                         ccGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCggg                                                 <  1:1080509/93‑1 (MQ=255)
                                         ccGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACACATAggg                                                 >  1:1522268/1‑93 (MQ=255)
                                          cgcgCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCgggg                                                >  2:700660/1‑93 (MQ=255)
                                          cgcgCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCgggg                                                >  1:4151374/1‑93 (MQ=255)
                                           gcgcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGc                                               >  2:2174472/1‑93 (MQ=255)
                                            cgcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACAAGACATATCGGGGcc                                              <  2:1192272/93‑1 (MQ=255)
                                             gcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCa                                             >  1:4459964/1‑93 (MQ=255)
                                             gcCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCa                                             >  1:3149827/1‑93 (MQ=255)
                                                aGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCg                                          <  2:4589002/93‑1 (MQ=255)
                                                aGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCg                                          >  2:3947380/1‑93 (MQ=255)
                                                 gTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGt                                         >  1:2701056/1‑93 (MQ=255)
                                                  ttGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTc                                        >  2:1861989/1‑93 (MQ=255)
                                                   tGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCa                                       <  2:1071590/93‑1 (MQ=255)
                                                   tGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCa                                       <  2:2305791/93‑1 (MQ=255)
                                                   tGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCa                                       <  2:1690370/93‑1 (MQ=255)
                                                    gATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAg                                      <  1:2549405/93‑1 (MQ=255)
                                                    gATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAg                                      <  2:661659/93‑1 (MQ=255)
                                                    gATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAg                                      >  1:292246/1‑93 (MQ=255)
                                                         gCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGgccgc                                 >  2:2319602/1‑93 (MQ=255)
                                                         gCTGAGAAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGgccgc                                 <  2:1220520/93‑1 (MQ=255)
                                                            gagaAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAg                              <  2:3221063/93‑1 (MQ=255)
                                                             agaAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGt                             <  1:3004186/93‑1 (MQ=255)
                                                             agaAACCAGTGACGGACCGAGCGCCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGt                             <  2:2633952/93‑1 (MQ=255)
                                                              gaAACCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGGGATAGTGCCAGAACACCAGACATATTGGGGCCATCGTTAGgccgcccgat                            >  2:2855634/1‑89 (MQ=255)
                                                               aaaCCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTa                           <  1:849580/93‑1 (MQ=255)
                                                                 aCCAGTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTaaa                         >  1:4092390/1‑93 (MQ=255)
                                                                     gTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCa                     >  2:3955600/1‑93 (MQ=255)
                                                                     gTGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCa                     >  2:2914101/1‑93 (MQ=255)
                                                                      tGACGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAg                    <  1:3358620/93‑1 (MQ=255)
                                                                         cGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCaa                 <  2:100902/93‑1 (MQ=255)
                                                                         cGGACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCaa                 <  2:1817066/93‑1 (MQ=255)
                                                                          ggACCGAGCGGCAGCCGGGGGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGGCAGCGTCAGGCCGCCAGTTCAAGTCAGCAAc                >  2:4490133/1‑93 (MQ=255)
                                                                           gACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACa               >  2:4345003/1‑93 (MQ=255)
                                                                           gACCGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACa               <  1:3027343/93‑1 (MQ=255)
                                                                             ccGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAgc             <  1:2737798/93‑1 (MQ=255)
                                                                             ccGAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAgc             <  2:2645609/93‑1 (MQ=255)
                                                                               gAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAgcgc           >  1:3925718/1‑93 (MQ=255)
                                                                               gAGCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAgcgc           <  2:785956/93‑1 (MQ=255)
                                                                                 gCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGCCATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCAt         <  1:3509135/93‑1 (MQ=255)
                                                                                 gCGGCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCAt         <  2:2590246/93‑1 (MQ=255)
                                                                                   ggCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGCCAGCAACAGCGCATTc       <  2:872505/93‑1 (MQ=255)
                                                                                    gCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTcc      >  2:40871/1‑93 (MQ=255)
                                                                                    gCAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTcc      <  1:4549822/93‑1 (MQ=255)
                                                                                    gCAGCCAGGTGTTATTAGCGAGATAGTGCCAGAACACCAGACATATCGGGACCAGCGTCAGGCCGCCAATTAAAGTCAGCAACAGCGCATTcc      <  1:3769671/93‑1 (MQ=255)
                                                                                     cAGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGAAACAGCGCATTCCa     >  2:3659567/1‑93 (MQ=255)
                                                                                      aGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCAc    >  2:250444/1‑93 (MQ=255)
                                                                                      aGCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCAc    >  1:3581442/1‑93 (MQ=255)
                                                                                       gCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCACg   >  1:2141303/1‑93 (MQ=255)
                                                                                       gCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCACg   >  1:4310894/1‑93 (MQ=255)
                                                                                       gCCGGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCACGa  >  1:1770438/1‑93 (MQ=255)
                                                                                          |                                                                                          
GACTGATAACAAACAGCAGCAAGTACCAGACCAGATGTCGCCCGCGCCAGTTGATTGGCTGAGAAACCAGTGACGGACCGAGCGGCAGCCAGGTGTTATTGGCGAGATAGTGCCAGAACACCAGACATATCGGGGCCAGCGTCAGGCCGCCAGTTAAAGTCAGCAACAGCGCATTCCACGA  >  NC_000913/3847175‑3847355

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