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,9600TC99.6% 805.4 / ‑4.8 249D216G (GAC→GGC) rfbXputative O‑antigen transporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (136/112);  total (136/113)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

Only 100 of 253 total aligned reads displayed.

GACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCA  >  NC_000913/2106870‑2107049
                                                                                          |                                                                                         
gACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCt                                                                                         <  2:1924225/93‑1 (MQ=255)
gACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCt                                                                                         >  2:302483/1‑93 (MQ=255)
gACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCt                                                                                         <  2:3154275/93‑1 (MQ=255)
  cGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGc                                                                                       <  2:3363171/93‑1 (MQ=255)
  cGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGc                                                                                       <  2:4350957/93‑1 (MQ=255)
   gTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAGCCGCCTGCt                                                                                      <  1:438096/93‑1 (MQ=255)
    ttGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTa                                                                                     >  2:601176/1‑93 (MQ=255)
      ggTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaaa                                                                                   <  2:3939727/93‑1 (MQ=255)
      ggTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaaa                                                                                   >  2:6200055/1‑93 (MQ=255)
      ggTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaaa                                                                                   <  1:2231529/93‑1 (MQ=255)
      ggTCAAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaaa                                                                                   <  1:3470174/93‑1 (MQ=255)
       gTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAg                                                                                  >  2:6299686/1‑93 (MQ=255)
        tCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGa                                                                                 >  2:1768294/1‑93 (MQ=255)
         ccAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGAt                                                                                <  2:2857090/93‑1 (MQ=255)
         ccAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGAt                                                                                <  1:3295548/93‑1 (MQ=255)
            gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                             <  1:5920256/93‑1 (MQ=255)
            gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                             <  2:2418846/93‑1 (MQ=255)
            gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                             <  1:5796861/93‑1 (MQ=255)
             atatAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGAc                                                                            >  2:1336222/1‑93 (MQ=255)
              tataAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACg                                                                           >  2:2855579/1‑93 (MQ=255)
              tataAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACg                                                                           >  1:4210718/1‑93 (MQ=255)
               ataAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGc                                                                          <  1:4152417/93‑1 (MQ=255)
                 aaaTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcac                                                                        >  1:4249843/1‑93 (MQ=255)
                  aaTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcaca                                                                       <  1:1957905/93‑1 (MQ=255)
                      ccAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTAt                                                                   >  2:5566128/1‑93 (MQ=255)
                       cAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGGACATTATg                                                                  >  1:2554661/1‑93 (MQ=255)
                        aGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGt                                                                 >  1:1014905/1‑93 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTa                                                                >  2:488991/1‑93 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTa                                                                >  2:4038632/1‑93 (MQ=255)
                          gATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTaa                                                               >  1:2015235/1‑93 (MQ=255)
                          gATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTaa                                                               <  1:1082599/93‑1 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGa                                                            >  1:3844444/1‑93 (MQ=255)
                              ataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGAt                                                           >  2:1524790/1‑93 (MQ=255)
                               taaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATa                                                          <  1:6635666/93‑1 (MQ=255)
                                aaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAg                                                         <  1:5417034/93‑1 (MQ=255)
                                 aCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGc                                                        >  1:6380050/1‑93 (MQ=255)
                                 aCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGATAAAGATCGACGCACATTATGTAATGATAGc                                                        <  1:3131642/93‑1 (MQ=255)
                                   tATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGt                                                      >  2:749964/1‑93 (MQ=255)
                                   tATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGt                                                      >  2:4990898/1‑93 (MQ=255)
                                    aTTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTa                                                     >  2:5733040/1‑93 (MQ=255)
                                    aTTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTa                                                     <  1:4597312/93‑1 (MQ=255)
                                      tCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAAc                                                   <  2:3711711/93‑1 (MQ=255)
                                      tCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAAc                                                   <  2:6127429/93‑1 (MQ=255)
                                        ccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACtt                                                 <  1:744720/93‑1 (MQ=255)
                                             gAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCa                                            >  1:1527834/1‑93 (MQ=255)
                                              aaTATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCaa                                           >  1:1364895/1‑93 (MQ=255)
                                               atatAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAAt                                          >  2:2972910/1‑93 (MQ=255)
                                               atatAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAAt                                          <  2:1038162/93‑1 (MQ=255)
                                                tataAACTAATGGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATc                                         >  1:2358942/1‑93 (MQ=255)
                                                tataAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATc                                         >  1:4898607/1‑93 (MQ=255)
                                                    aaCTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAAc                                     >  1:1428441/1‑93 (MQ=255)
                                                    aaCTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAAc                                     <  1:1284216/93‑1 (MQ=255)
                                                     aCTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAAcc                                    <  1:3104084/93‑1 (MQ=255)
                                                      cTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCt                                   >  1:2899029/1‑93 (MQ=255)
                                                        aaTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTc                                 <  2:4515032/93‑1 (MQ=255)
                                                         aTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCa                                >  2:3417671/1‑93 (MQ=255)
                                                         aTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCa                                <  1:6526866/93‑1 (MQ=255)
                                                         aTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCa                                >  2:267176/1‑93 (MQ=255)
                                                           aGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATg                              >  1:1883318/1‑93 (MQ=255)
                                                            gCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGa                             >  2:2007528/1‑93 (MQ=255)
                                                            acAGATGTGGAGATAAAAACATGATAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTGCCAATCCAACCTTCATGa                             <  1:6025042/92‑1 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                           <  1:832888/93‑1 (MQ=255)
                                                               gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                          >  2:4801150/1‑93 (MQ=255)
                                                               gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                          >  1:4985789/1‑93 (MQ=255)
                                                                 tgtgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaacaac                        >  2:4217566/1‑93 (MQ=255)
                                                                 tgtgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaacaac                        >  2:5527680/1‑93 (MQ=255)
                                                                    ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                     <  1:5467950/93‑1 (MQ=255)
                                                                      aaaTAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAg                   <  2:2182403/93‑1 (MQ=255)
                                                                        aTAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCt                 >  1:94090/1‑93 (MQ=255)
                                                                         tAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTa                <  2:5835874/93‑1 (MQ=255)
                                                                         aaaaaaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTa                <  2:4641802/92‑1 (MQ=255)
                                                                           aaaaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGt              >  1:3654996/1‑93 (MQ=255)
                                                                           aaaaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGt              >  1:5550458/1‑93 (MQ=255)
                                                                            aaaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTg             >  1:5187324/1‑93 (MQ=255)
                                                                              aCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCa           >  1:2057866/1‑93 (MQ=255)
                                                                               cATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaa          >  2:3831711/1‑93 (MQ=255)
                                                                               cATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaa          <  1:6104335/93‑1 (MQ=255)
                                                                               cATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGCGCaa          >  1:2976541/1‑93 (MQ=255)
                                                                                aTGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaat         >  2:5785438/1‑93 (MQ=255)
                                                                                aTGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaat         <  2:1546072/93‑1 (MQ=255)
                                                                                aTGAAAACCGCCTGCGAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaat         <  2:1301690/93‑1 (MQ=255)
                                                                                 tGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaata        >  2:6574904/1‑93 (MQ=255)
                                                                                 agAAAACCGCCGGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaata        <  1:1193721/92‑1 (MQ=255)
                                                                                  gAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataa       <  2:861289/93‑1 (MQ=255)
                                                                                  gAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataa       <  1:2457023/93‑1 (MQ=255)
                                                                                  gAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataa       >  1:4005715/1‑93 (MQ=255)
                                                                                  gAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataa       <  1:5249902/93‑1 (MQ=255)
                                                                                   aaaaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataat      <  1:4299236/93‑1 (MQ=255)
                                                                                   aaaaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATCAACAACGATAGCTAGTGCaataat      >  1:6589308/1‑93 (MQ=255)
                                                                                   aaaaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACAAATCCAACCTTCATGAACAACGATAGCTAGTGCaataat      <  1:4061861/93‑1 (MQ=255)
                                                                                    aaaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAAtt     >  2:3753958/1‑93 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc    >  2:5538918/1‑93 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc    <  2:3644327/93‑1 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc    >  1:4768770/1‑93 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc    >  1:124979/1‑93 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc    <  1:983388/93‑1 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc    >  2:6697278/1‑93 (MQ=255)
                                                                                      aCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTcc   >  2:5091776/1‑93 (MQ=255)
                                                                                       ccGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCa  <  1:179293/93‑1 (MQ=255)
                                                                                       ccGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCa  <  1:5644310/93‑1 (MQ=255)
                                                                                          |                                                                                         
GACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCA  >  NC_000913/2106870‑2107049

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