Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,106,960 T→C D216G (GAC→GGC)  rfbX ← putative O‑antigen transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,106,9600TC98.5% 407.8 / ‑2.5 133D216G (GAC→GGC) rfbXputative O‑antigen transporter
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (80/51);  total (81/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

Only 100 of 135 total aligned reads displayed.

GACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCC  >  NC_000913/2106870‑2107048
                                                                                          |                                                                                        
  cGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGc                                                                                      <  2:332441/93‑1 (MQ=255)
   gTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCt                                                                                     >  2:4068246/1‑93 (MQ=255)
   gTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCt                                                                                     >  1:4026254/1‑93 (MQ=255)
    ttGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTa                                                                                    <  2:2249749/93‑1 (MQ=255)
    ttGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTa                                                                                    >  1:1786588/1‑93 (MQ=255)
     tGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaa                                                                                   <  2:4387400/93‑1 (MQ=255)
      ggTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATGTAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaaa                                                                                  <  2:4619787/93‑1 (MQ=255)
      ggTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaaa                                                                                  <  1:3089221/93‑1 (MQ=255)
      ggTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACACGAAAACCGCCTGCTaaa                                                                                  <  1:2543315/93‑1 (MQ=255)
       gTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAg                                                                                 >  1:4072032/1‑93 (MQ=255)
        tCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGa                                                                                <  2:597136/93‑1 (MQ=255)
         ccAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGAt                                                                               <  2:3231185/93‑1 (MQ=255)
          cAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATc                                                                              <  2:1495920/93‑1 (MQ=255)
           aGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCg                                                                             >  1:2707941/1‑93 (MQ=255)
            gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                            >  1:2921758/1‑93 (MQ=255)
            gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                            >  1:3095685/1‑93 (MQ=255)
            gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGGGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                            <  1:2468241/93‑1 (MQ=255)
              tataAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACg                                                                          >  2:2970725/1‑93 (MQ=255)
               ataAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGc                                                                         >  1:962598/1‑93 (MQ=255)
                taAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcg                                                                        >  1:898730/1‑92 (MQ=255)
                taAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcg                                                                        >  2:4027523/1‑92 (MQ=255)
                 aaaTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcgc                                                                       >  1:4000095/1‑91 (MQ=255)
                 aaaTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcgc                                                                       <  1:434801/93‑3 (MQ=255)
                  aaTCCCAGGATAATACCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcgca                                                                      <  1:2033720/93‑1 (MQ=255)
                      ccAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTAt                                                                  >  2:1751500/1‑93 (MQ=255)
                       cAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATg                                                                 >  2:3275547/1‑93 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTa                                                               <  2:1568440/93‑1 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTa                                                               >  1:4471384/1‑93 (MQ=255)
                           ataataACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAAt                                                             >  1:1958819/1‑93 (MQ=255)
                            taataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATg                                                            >  2:130509/1‑93 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           <  2:4348244/93‑1 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           >  2:2001059/1‑93 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           <  1:710490/93‑1 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           >  1:4500021/1‑93 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           <  1:2870861/93‑1 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGa                                                           <  2:1565436/93‑1 (MQ=255)
                              ataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGAt                                                          >  1:2635852/1‑93 (MQ=255)
                               taaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAATATCGACGCGCATTATGTAATGATa                                                         >  2:4539698/1‑93 (MQ=255)
                               taaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATa                                                         <  1:185119/93‑1 (MQ=255)
                               taaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATa                                                         >  2:3552136/1‑93 (MQ=255)
                                aaCTATTCCCGTGGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAg                                                        >  1:528614/1‑93 (MQ=255)
                                  cTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCg                                                      >  2:2680368/1‑93 (MQ=255)
                                  cTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCg                                                      >  1:1655777/1‑93 (MQ=255)
                                  cTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCg                                                      >  2:2808981/1‑93 (MQ=255)
                                    aTTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTa                                                    >  2:2624619/1‑93 (MQ=255)
                                       cccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACt                                                 >  1:3906053/1‑93 (MQ=255)
                                       cccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACt                                                 <  2:2317930/93‑1 (MQ=255)
                                        ccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACtt                                                >  1:588377/1‑93 (MQ=255)
                                           tAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTAc                                             >  2:3671648/1‑93 (MQ=255)
                                              aaTATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCaa                                          >  2:2826448/1‑93 (MQ=255)
                                              aaTATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCaa                                          >  2:854349/1‑93 (MQ=255)
                                               atatAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAAt                                         <  1:3480153/93‑1 (MQ=255)
                                                tataAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATc                                        >  1:2038745/1‑93 (MQ=255)
                                                 ataAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATcc                                       <  2:2637341/93‑1 (MQ=255)
                                                 ataAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATcc                                       >  2:2320193/1‑93 (MQ=255)
                                                  taAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCa                                      >  1:1901619/1‑93 (MQ=255)
                                                  taAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCa                                      <  1:1644025/93‑1 (MQ=255)
                                                    aaCTAATAGCAGATGTGGCAATAAAAACATGAAAACCGCCTGCTCAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAAc                                    >  2:707954/1‑93 (MQ=255)
                                                      cTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCt                                  <  2:4310940/93‑1 (MQ=255)
                                                      cTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCt                                  <  1:1843078/93‑1 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                 >  1:4309260/1‑93 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                 >  2:1209759/1‑93 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                 >  1:2251113/1‑93 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                 >  2:634151/1‑93 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCACCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                 >  2:4042853/1‑93 (MQ=255)
                                                         aTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCa                               <  1:2988887/93‑1 (MQ=255)
                                                         aTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCa                               >  2:2478997/1‑93 (MQ=255)
                                                           aGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATg                             >  2:3599834/1‑93 (MQ=255)
                                                            gCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGa                            >  1:3444889/1‑93 (MQ=255)
                                                             cAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaa                           <  2:3212479/93‑1 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                          >  2:4591200/1‑93 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                          <  2:954730/93‑1 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                          >  2:1998150/1‑93 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                          >  2:1580096/1‑93 (MQ=255)
                                                               gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                         >  1:1315150/1‑93 (MQ=255)
                                                               gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                         >  1:2660029/1‑93 (MQ=255)
                                                               gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                         >  1:4304209/1‑93 (MQ=255)
                                                                aTGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaacaa                        <  1:4551270/93‑1 (MQ=255)
                                                                aTGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaacaa                        <  1:835393/93‑1 (MQ=255)
                                                                 tgtgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAATCTTCATGaacaac                       >  2:3028899/1‑93 (MQ=255)
                                                                    ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                    <  2:1829571/93‑1 (MQ=255)
                                                                    ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                    <  1:2330033/93‑1 (MQ=255)
                                                                    ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                    <  1:2647455/93‑1 (MQ=255)
                                                                     gAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATa                   <  1:2477070/93‑1 (MQ=255)
                                                                       aaTAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGc                 >  2:3672294/1‑93 (MQ=255)
                                                                       aaTAAAAACATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGc                 >  1:308977/1‑93 (MQ=255)
                                                                         aaaaaaCATGAAACCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTa               <  2:289811/92‑1 (MQ=255)
                                                                          aaaaaCATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAg              <  1:3223861/93‑1 (MQ=255)
                                                                          aaaaaCATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAg              >  1:2350715/1‑93 (MQ=255)
                                                                           aaaaCATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGt             >  2:3631466/1‑93 (MQ=255)
                                                                           aaaaCATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGt             >  1:2751081/1‑93 (MQ=255)
                                                                              aCATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCa          >  1:4290687/1‑93 (MQ=255)
                                                                               cATGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaa         >  1:1418110/1‑93 (MQ=255)
                                                                                aTGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaat        >  1:1268232/1‑93 (MQ=255)
                                                                                 tGAAAACCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaata       >  2:2435888/1‑93 (MQ=255)
                                                                                   aaaaCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataat     <  2:120582/93‑1 (MQ=255)
                                                                                   aaaaCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataat     <  2:3378531/93‑1 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc   <  2:3072561/93‑1 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc   >  2:2416528/1‑93 (MQ=255)
                                                                                      aCCGCCTGCTAAAGATCGACGCGCATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTcc  >  2:3908749/1‑93 (MQ=255)
                                                                                          |                                                                                        
GACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCC  >  NC_000913/2106870‑2107048

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