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.5% 681.6 / ‑5.1 211D216G (GAC→GGC) rfbXputative O‑antigen transporter
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (134/76);  total (135/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 213 total aligned reads displayed.

GACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCA  >  NC_000913/2106870‑2107049
                                                                                          |                                                                                         
gACGTTGGTCCAGATATAAATCCCAGTAGAATCACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCt                                                                                         <  1:3163431/93‑1 (MQ=255)
gACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCt                                                                                         <  2:2700500/93‑1 (MQ=255)
gACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCt                                                                                         >  2:2395578/1‑93 (MQ=255)
  cGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGc                                                                                       <  1:1750430/93‑1 (MQ=255)
   gTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCt                                                                                      <  2:1466754/93‑1 (MQ=255)
   gTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCt                                                                                      >  1:3225025/1‑93 (MQ=255)
      ggTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTaaa                                                                                   >  1:1556184/1‑93 (MQ=255)
       gTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAg                                                                                  >  1:2999546/1‑93 (MQ=255)
        tCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGa                                                                                 <  2:4473394/93‑1 (MQ=255)
        tCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGa                                                                                 <  1:1381089/93‑1 (MQ=255)
         ccAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGAt                                                                                <  2:1151385/93‑1 (MQ=255)
          cAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATc                                                                               >  1:2063261/1‑93 (MQ=255)
           aGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCg                                                                              >  2:2474305/1‑93 (MQ=255)
           aGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCg                                                                              >  1:3721009/1‑93 (MQ=255)
             atatAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGAc                                                                            <  1:1631216/93‑1 (MQ=255)
              tataAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACg                                                                           >  2:2601040/1‑93 (MQ=255)
              tataAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACg                                                                           >  1:3310036/1‑93 (MQ=255)
                taAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGca                                                                         >  2:513318/1‑93 (MQ=255)
                  aaTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcaca                                                                       <  1:4639994/93‑1 (MQ=255)
                    tCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACAtt                                                                     >  2:4184869/1‑93 (MQ=255)
                    tCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACAtt                                                                     <  2:3809020/93‑1 (MQ=255)
                     cccAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTa                                                                    >  1:2852592/1‑93 (MQ=255)
                       cAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATg                                                                  >  1:3274276/1‑93 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTa                                                                >  1:2694016/1‑93 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTa                                                                <  1:2654288/93‑1 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTa                                                                <  2:666143/93‑1 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTa                                                                >  1:3391889/1‑93 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTa                                                                >  2:4276007/1‑93 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTa                                                                <  1:1845852/93‑1 (MQ=255)
                         ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTa                                                                >  2:2017003/1‑93 (MQ=255)
                          gATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTaa                                                               <  2:3371958/93‑1 (MQ=255)
                          gATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTaa                                                               >  2:710189/1‑93 (MQ=255)
                           ataataACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAAt                                                              <  2:3575183/93‑1 (MQ=255)
                            taataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATg                                                             >  2:1691418/1‑93 (MQ=255)
                            taataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATg                                                             >  2:4610221/1‑93 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGa                                                            >  2:851205/1‑93 (MQ=255)
                             aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGa                                                            >  2:3118539/1‑93 (MQ=255)
                              ataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGAt                                                           >  2:4145875/1‑93 (MQ=255)
                              ataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGAt                                                           >  2:3744500/1‑93 (MQ=255)
                              ataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGAt                                                           >  2:1512130/1‑93 (MQ=255)
                               taaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACAGCCTGCTAAAGATCGACGCACATTATGTAATGATa                                                          >  2:2361489/1‑93 (MQ=255)
                                aaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGGAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAg                                                         >  1:3057314/1‑93 (MQ=255)
                                aaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAg                                                         <  2:4021026/93‑1 (MQ=255)
                                 aCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGc                                                        <  2:2912432/93‑1 (MQ=255)
                                 aCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGc                                                        >  1:4028619/1‑93 (MQ=255)
                                 aCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGc                                                        <  2:2017598/93‑1 (MQ=255)
                                 aCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGc                                                        >  2:1984801/1‑93 (MQ=255)
                                  cTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCg                                                       >  1:1911412/1‑93 (MQ=255)
                                   tATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGt                                                      >  2:3904706/1‑93 (MQ=255)
                                     ttCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTaa                                                    >  2:4125537/1‑93 (MQ=255)
                                     ttCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTaa                                                    >  2:398631/1‑93 (MQ=255)
                                     ttCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTaa                                                    >  2:155023/1‑93 (MQ=255)
                                     ttCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTaa                                                    >  2:155143/1‑93 (MQ=255)
                                      tCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAAc                                                   <  2:3185052/93‑1 (MQ=255)
                                       cccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACt                                                  <  2:3520455/93‑1 (MQ=255)
                                         cGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACttt                                                >  1:2327224/1‑93 (MQ=255)
                                         cGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACttt                                                <  2:4410793/93‑1 (MQ=255)
                                          gTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTa                                               >  2:3265496/1‑93 (MQ=255)
                                          gTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTa                                               >  1:4544363/1‑93 (MQ=255)
                                           tAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTAc                                              <  2:2409249/93‑1 (MQ=255)
                                           tAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTAc                                              <  1:4345224/93‑1 (MQ=255)
                                           tAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTAc                                              <  2:4520947/93‑1 (MQ=255)
                                            aGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTAcc                                             >  2:1845211/1‑93 (MQ=255)
                                            aGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTAcc                                             <  2:1878452/93‑1 (MQ=255)
                                             gAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCa                                            >  2:2326568/1‑93 (MQ=255)
                                             gAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCa                                            >  2:1336540/1‑93 (MQ=255)
                                              aaTATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCaa                                           >  1:1002012/1‑93 (MQ=255)
                                               atatAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAAt                                          >  1:1484802/1‑93 (MQ=255)
                                                 ataAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATcc                                        <  1:859846/93‑1 (MQ=255)
                                                  taAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCa                                       >  1:1042646/1‑93 (MQ=255)
                                                   aaaCTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCac                                      >  1:3336256/1‑92 (MQ=255)
                                                    aaCTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAAc                                     <  1:2478728/93‑1 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                  >  1:3550933/1‑93 (MQ=255)
                                                       tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATACAACCtt                                  >  1:1667281/1‑93 (MQ=255)
                                                         aTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCa                                <  1:197673/93‑1 (MQ=255)
                                                         aTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCa                                >  1:2244271/1‑93 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACTAATCCAACCTTCATGaaa                           >  2:2231495/1‑92 (MQ=255)
                                                              aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                           >  2:4614592/1‑93 (MQ=255)
                                                               gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                          >  1:2029760/1‑93 (MQ=255)
                                                                 tgtgGAAATAAACACATGAAAAACGCCTGCTAAAGATCGTCGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCACGaacaac                        <  1:677329/93‑1 (MQ=255)
                                                                 tgtgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaacaac                        >  1:3747994/1‑93 (MQ=255)
                                                                 tgtgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaacaac                        >  2:1283050/1‑93 (MQ=255)
                                                                 tgaggAAATAAAAACATGAAAACCGCCTGCTAAAGATAGACGCACATTATGTAATGATAGCGAAAGTTTACCAATCCAACCTTCATGaacaac                        >  2:3866321/4‑93 (MQ=255)
                                                                   tgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGa                      >  2:2659851/1‑93 (MQ=255)
                                                                    ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                     >  2:4372218/1‑93 (MQ=255)
                                                                      aaaTAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAg                   >  2:2577356/1‑93 (MQ=255)
                                                                         tAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTa                <  2:3802534/93‑1 (MQ=255)
                                                                          aaaaaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAg               <  1:744349/93‑1 (MQ=255)
                                                                           aaaaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGt              >  2:2098722/1‑93 (MQ=255)
                                                                             aaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGc            >  1:550952/1‑93 (MQ=255)
                                                                              aCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTCACTTTACCACTCCAACCTTCATGAACAACGATAGCTAGTGCa           >  1:3095155/1‑93 (MQ=255)
                                                                                aTGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaat         >  2:1260497/1‑93 (MQ=255)
                                                                                aTGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaat         <  1:2963007/93‑1 (MQ=255)
                                                                                 tGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaata        >  2:1022690/1‑93 (MQ=255)
                                                                                   aaaaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataat      <  2:1323614/93‑1 (MQ=255)
                                                                                   aaaaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaataat      >  1:2649487/1‑93 (MQ=255)
                                                                                    aaaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAAtt     <  2:566068/93‑1 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc    >  2:1356166/1‑93 (MQ=255)
                                                                                     aaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc    <  2:3653583/93‑1 (MQ=255)
                                                                                       ccGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCa  <  1:121509/93‑1 (MQ=255)
                                                                                          |                                                                                         
GACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCA  >  NC_000913/2106870‑2107049

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