Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,485,264 G→A 16.6% W297* (TGG→TGA evgS → hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,485,2640GA16.6% 475.2 / 71.4 199W297* (TGG→TGAevgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base G (81/83);  new base A (21/12);  total (102/97)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 209 total aligned reads displayed.

AAATGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGA  >  NC_000913/2485165‑2485362
                                                                                                   |                                                                                                  
aaaTGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTgg                                                                                                    <  1:3149670/100‑1 (MQ=255)
aaaTGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTgg                                                                                                    >  1:3711041/1‑100 (MQ=255)
aaaTGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTgg                                                                                                    <  2:3123523/100‑1 (MQ=255)
 aaTGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGa                                                                                                   >  1:3584240/1‑100 (MQ=255)
  aTGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAAt                                                                                                  <  1:2452027/100‑1 (MQ=255)
    gAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTa                                                                                                <  1:3787603/100‑1 (MQ=255)
      aGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAag                                                                                              >  2:3526580/1‑100 (MQ=255)
      aGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAag                                                                                              <  2:3206607/100‑1 (MQ=255)
      aGTGCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGGAAAACAGTGGAtt                                                                                                 <  2:2579726/97‑1 (MQ=255)
       gTTCGGTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCAATGAACATGAAAAACAGTGGATTAagc                                                                                             >  2:3280418/1‑100 (MQ=255)
       gTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAagc                                                                                             <  1:656640/100‑1 (MQ=255)
        ttCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAagca                                                                                            <  1:3305125/100‑1 (MQ=255)
        ttCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAagca                                                                                            <  1:3606248/100‑1 (MQ=255)
         tCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAagcag                                                                                           >  1:2681546/1‑100 (MQ=255)
          cGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAagcagc                                                                                          <  1:2132802/100‑1 (MQ=255)
             tATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATc                                                                                       >  2:277537/1‑100 (MQ=255)
              aTGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATcc                                                                                      >  1:1577518/1‑100 (MQ=255)
               ttaaGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATccc                                                                                     >  1:2819978/3‑100 (MQ=255)
               tGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATccc                                                                                     <  1:87013/100‑1 (MQ=255)
               tGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATccc                                                                                     >  1:988689/1‑100 (MQ=255)
                gAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCa                                                                                    >  2:3071694/1‑100 (MQ=255)
                   gTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAAtt                                                                                 <  2:3496245/100‑1 (MQ=255)
                     aTCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTa                                                                               >  2:1361540/1‑100 (MQ=255)
                       cacaAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTaaa                                                                             >  2:1876560/1‑100 (MQ=255)
                       cacaAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTaaa                                                                             <  2:747308/100‑1 (MQ=255)
                        acaAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAg                                                                            >  2:3939589/1‑100 (MQ=255)
                        acaAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAg                                                                            >  1:2510778/1‑100 (MQ=255)
                         caAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAgg                                                                           <  2:3905699/100‑1 (MQ=255)
                         caAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAgg                                                                           <  1:2406586/100‑1 (MQ=255)
                          aaaaTTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGt                                                                          >  1:2373578/1‑100 (MQ=255)
                            aaTTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGc                                                                        <  2:3779127/100‑1 (MQ=255)
                                ggCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGa                                                                    <  2:1080161/100‑1 (MQ=255)
                                ggCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGa                                                                    <  1:77379/100‑1 (MQ=255)
                                 gCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGaa                                                                   >  1:1043789/1‑100 (MQ=255)
                                 gCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGaa                                                                   >  1:1292592/1‑100 (MQ=255)
                                  ccgatctcAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGaaa                                                                  <  1:3920828/93‑1 (MQ=255)
                                  cTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGaaa                                                                  >  2:429744/1‑100 (MQ=255)
                                   ttGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGaaaa                                                                 >  1:3035296/1‑100 (MQ=255)
                                   ttGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGaaaa                                                                 >  1:977745/1‑100 (MQ=255)
                                    tGATACAGGAAAGATGACCTTTCTGAACAAACCATTAGAACTCACTGAACATGCAATACAGTTGATTAAGCATCATTCAAATATATAGTTCCTGGtaaa                                                                 >  2:3443278/1‑99 (MQ=255)
                                     gaaaCAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATc                                                               >  2:2980244/4‑100 (MQ=255)
                                       tACAGGAAACCTGGCCTTTCTGAACAAACCATTGGAACTCACTGAAAATGAAAAACAGCGGATTAAGCAACATCCCAATTTACAGGTGCTGGAAACTcc                                                              >  1:3889721/1‑99 (MQ=255)
                                       tACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCt                                                             <  1:2490021/100‑1 (MQ=255)
                                       tACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCt                                                             >  1:3071808/1‑100 (MQ=255)
                                       tACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCt                                                             <  1:3512984/100‑1 (MQ=255)
                                        aCAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCtt                                                            <  1:1770553/100‑1 (MQ=255)
                                         cAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTa                                                           >  2:1723721/1‑100 (MQ=255)
                                         cAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTa                                                           >  1:2262635/1‑100 (MQ=255)
                                         cAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGaaaagatcgg                                                           >  2:3920798/1‑94 (MQ=255)
                                          aGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTAc                                                          <  2:3588807/100‑1 (MQ=255)
                                            gAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTc                                                        <  1:3394445/100‑1 (MQ=255)
                                            gAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTc                                                        <  1:2985479/100‑1 (MQ=255)
                                             aaaCCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCg                                                       <  1:3239155/100‑1 (MQ=255)
                                                ccTGGCCTTTCTGAACAAACCATTAGAACTCACTGATCATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGAGCTGGAAAATCCTTACTCGccc                                                    >  2:669265/1‑100 (MQ=255)
                                                   ggCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGcccccc                                                 <  2:3552419/100‑1 (MQ=255)
                                                   ggCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGcccccc                                                 <  1:3348615/100‑1 (MQ=255)
                                                      cTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTAt                                              <  2:3361660/100‑1 (MQ=255)
                                                      cTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCCCCCCCCTAt                                              >  2:2864403/1‑100 (MQ=255)
                                                       tttCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTAtt                                             >  1:4174106/1‑100 (MQ=255)
                                                          cTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGTTGTTGGAAAATCCTTACTCGCCCCCCTATTCTa                                          <  2:2744865/100‑1 (MQ=255)
                                                            gAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATg                                        <  1:1154398/100‑1 (MQ=255)
                                                            gAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATg                                        >  2:1800517/1‑100 (MQ=255)
                                                            gAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATg                                        >  1:3807469/1‑100 (MQ=255)
                                                              aCAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGAc                                      >  2:4098146/1‑100 (MQ=255)
                                                              aCAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGAc                                      >  2:1172736/1‑100 (MQ=255)
                                                               cAAAGCATTAGAAATCACTGAACATGAAAAACAGCGAATTAAGCAGCATCCCACTTTAAAGGTGCTGGAAAATCCTTACTCGCCCCGCTATTCTATGACg                                     <  2:1336589/100‑1 (MQ=255)
                                                               cAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACg                                     <  2:454416/100‑1 (MQ=255)
                                                               cAAACCATTAGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACg                                     <  1:2011116/100‑1 (MQ=255)
                                                               cAAACCATTAGAACGCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACg                                     <  2:1580534/100‑1 (MQ=255)
                                                                aaaCCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACgg                                    <  1:1671732/100‑1 (MQ=255)
                                                                 aaCCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGa                                   >  1:2774745/1‑100 (MQ=255)
                                                                   ccATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATg                                 >  1:482173/1‑100 (MQ=255)
                                                                   ccATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATg                                 >  2:4152050/1‑100 (MQ=255)
                                                                     aTTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGaa                               >  1:2928929/1‑100 (MQ=255)
                                                                     aTTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGaa                               >  1:593525/1‑100 (MQ=255)
                                                                      ttAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGaaa                              <  1:1128239/100‑1 (MQ=255)
                                                                        aGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAAt                            >  1:3213705/1‑100 (MQ=255)
                                                                        aGAACTCACTGAACATGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAAt                            >  1:326195/1‑100 (MQ=255)
                                                                         gAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATg                           >  1:909861/1‑100 (MQ=255)
                                                                           aCTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATgg                          >  2:2886118/1‑99 (MQ=255)
                                                                             tCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTc                       <  1:174171/100‑1 (MQ=255)
                                                                              cACTGAATATGAAACACAGTGGATTAAGCAGCATCCCAATTTAAAGATGCTGGAAAATCCTTCCTCGCCACCCTATTCTATGACGGATGAAAATGGCTCg                      >  2:1883231/1‑100 (MQ=255)
                                                                              cACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCg                      >  2:226866/1‑100 (MQ=255)
                                                                              cACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCg                      >  1:4070657/1‑100 (MQ=255)
                                                                                    aCATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCgg                >  1:2599812/1‑100 (MQ=255)
                                                                                     cgtGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCggg               >  1:3862860/3‑100 (MQ=255)
                                                                                     cATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTc                  >  2:1819703/1‑97 (MQ=255)
                                                                                      aTGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCgggg              <  1:4151367/100‑1 (MQ=255)
                                                                                       tGAAAAACAGTGGATTTAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGc             <  2:2431705/100‑1 (MQ=255)
                                                                                       tGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGc             >  2:2157366/1‑100 (MQ=255)
                                                                                       tGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGc             >  1:3397251/1‑100 (MQ=255)
                                                                                       tGAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGc             <  2:2404422/100‑1 (MQ=255)
                                                                                        gAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCg            <  1:3987972/100‑1 (MQ=255)
                                                                                        gAAAAACAGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCg            >  2:1464287/1‑100 (MQ=255)
                                                                                        gAAAAACAGTGAAGTAAGCAGCATCACAATTTAAAAGTGCTGGAAAATGCTTACTCGCGCCCCTAGTCTTTGACGGATGAATATGGCTCGGTTCGGGGCg            <  2:2548696/100‑1 (MQ=255)
                                                                                          aaaaCAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGATCGGTTCggg               >  2:4003017/1‑95 (MQ=255)
                                                                                           aaaCAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTa         <  1:2365389/100‑1 (MQ=255)
                                                                                               aGTGAATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATggg     >  1:569222/1‑100 (MQ=255)
                                                                                                gTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATgggg    <  1:1472317/100‑1 (MQ=255)
                                                                                                gTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATgggg    <  2:2338110/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
AAATGAAGTTCGCTATGAAGTATCACAAAATTGGCTTGATACAGGAAACCTGGCCTTTCTGAACAAACCATTAGAACTCACTGAACATGAAAAACAGTGGATTAAGCAGCATCCCAATTTAAAGGTGCTGGAAAATCCTTACTCGCCCCCCTATTCTATGACGGATGAAAATGGCTCGGTTCGGGGCGTTATGGGGGA  >  NC_000913/2485165‑2485362

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