Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,792,026 A→G N138S (AAC→AGC)  rfaD → ADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,792,0260AG99.3% 1374.8 / ‑2.6 436N138S (AAC→AGC) rfaDADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (191/242);  total (193/243)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

Only 100 of 452 total aligned reads displayed.

AAATCCCGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATT  >  NC_000913/3791936‑3792115
                                                                                          |                                                                                         
   tCCCGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCt                                                                                      <  1:3951182/93‑1 (MQ=255)
     ccGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTAc                                                                                    >  2:7398371/1‑93 (MQ=255)
     ccGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTAc                                                                                    <  2:609496/93‑1 (MQ=255)
      cGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTTCg                                                                                   >  2:2452840/1‑93 (MQ=255)
      cGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACg                                                                                   <  1:1916372/93‑1 (MQ=255)
      cGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACg                                                                                   <  1:2081821/93‑1 (MQ=255)
       gTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACgg                                                                                  <  1:7886010/93‑1 (MQ=255)
        ttCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGt                                                                                 >  2:5702414/1‑93 (MQ=255)
        ttCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGAGTATACGGt                                                                                 >  1:6418955/1‑93 (MQ=255)
               aCGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGCCTACGGTTACTCaa                                                                          >  1:466317/1‑93 (MQ=255)
                cGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCaaa                                                                         >  1:2393894/1‑93 (MQ=255)
                 gCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCaaaa                                                                        >  1:6159088/1‑93 (MQ=255)
                    tcttcCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTc                                                                     >  1:4566574/1‑93 (MQ=255)
                      ttcCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCt                                                                   <  1:4956328/93‑1 (MQ=255)
                        ccGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGt                                                                 <  1:5315328/93‑1 (MQ=255)
                         cGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGtt                                                                <  1:1473213/93‑1 (MQ=255)
                         cGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGtt                                                                <  2:1008944/93‑1 (MQ=255)
                          gCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGttt                                                               <  2:6551864/93‑1 (MQ=255)
                          gCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGttt                                                               <  2:1914848/93‑1 (MQ=255)
                          ccAGCCACCTCGGGCGGACGCCCCTCCGACTTTATTGAATCCCGCGAGTACGAAAACCCGTTGAGCGTCTACGGTTACTCAAAATTCCTGttt                                                               <  1:6601985/92‑1 (MQ=255)
                           cAGCCACCTAGGGCGGGCGCCCCTCCGAGTTTATTGGAAGCCGCGAGTACGAAACACCGCTGAGCGACTACGGTTACTCAAAATTCCtgtgtg                                                              >  2:1198091/1‑93 (MQ=255)
                            aGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTtga                                                             <  2:1976560/93‑1 (MQ=255)
                            aGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGGTGAGCGTCTACGGTTACTCAAAATTCCTGCTtga                                                             >  1:3172783/1‑93 (MQ=255)
                             gCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTtgat                                                            <  1:2681342/93‑1 (MQ=255)
                              ccACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTtgatg                                                           >  1:7127984/1‑93 (MQ=255)
                               cACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTtgatga                                                          <  2:7616462/93‑1 (MQ=255)
                                cccTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGaa                                                         <  2:4692599/92‑1 (MQ=255)
                                 ccTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAATAACCGTTGAGCGTCCACGGTTACTCAAAATTCCTGTTTGATGAat                                                        >  1:2215201/1‑93 (MQ=255)
                                 ccTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAat                                                        <  2:3254534/93‑1 (MQ=255)
                                  cTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAata                                                       <  1:5620165/93‑1 (MQ=255)
                                  cTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAata                                                       <  1:6647000/93‑1 (MQ=255)
                                   tACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAatat                                                      <  2:7823990/93‑1 (MQ=255)
                                   tACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGATGAatat                                                      <  1:7211397/93‑1 (MQ=255)
                                    aCGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATg                                                     >  2:4181214/1‑93 (MQ=255)
                                    aCGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATg                                                     <  2:6944254/93‑1 (MQ=255)
                                      ggcggACGCGCCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCATCTACGGTTACTCAAAATTCCTGTTTGATGAATATGtt                                                   <  2:316074/93‑1 (MQ=255)
                                      ggcggACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGtt                                                   <  2:7472112/93‑1 (MQ=255)
                                      ggcggACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGtt                                                   <  1:3783890/93‑1 (MQ=255)
                                      ggcggACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGtt                                                   >  1:5891707/1‑93 (MQ=255)
                                       gcggACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTc                                                  >  2:5313220/1‑93 (MQ=255)
                                       gcggACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTc                                                  <  1:6039924/93‑1 (MQ=255)
                                        cggACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCg                                                 >  1:69394/1‑93 (MQ=255)
                                        cggACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCg                                                 <  2:4330155/93‑1 (MQ=255)
                                          gACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGTGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTAGTc                                               >  2:4598268/1‑93 (MQ=255)
                                           aCGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCa                                              <  2:1793429/93‑1 (MQ=255)
                                            cGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCaa                                             <  2:2721355/93‑1 (MQ=255)
                                             gCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCaaa                                            <  1:2861155/93‑1 (MQ=255)
                                             gCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCaaa                                            >  2:3803804/1‑93 (MQ=255)
                                               aCCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATc                                          >  1:3614838/1‑93 (MQ=255)
                                                ccTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAGATct                                         >  2:1585859/1‑92 (MQ=255)
                                                 cTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATcct                                        >  1:1973950/1‑93 (MQ=255)
                                                 cTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATcct                                        <  1:2621042/93‑1 (MQ=255)
                                                 cTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTCCGTCAAATcct                                        >  2:5620013/1‑93 (MQ=255)
                                                 cTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTCAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATcct                                        >  1:4430273/1‑93 (MQ=255)
                                                  tCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATcctg                                       >  2:473197/1‑93 (MQ=255)
                                                    cgccTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATcctgcc                                     <  1:2060966/90‑1 (MQ=255)
                                                    cGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATcctgcc                                     <  2:730376/93‑1 (MQ=255)
                                                     gACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATcctgcct                                    <  2:5058958/93‑1 (MQ=255)
                                                     aaCTTTATTGAGTCCCGCGAGTACGAAAAACCGTTGAGCGTCTTCGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATcctgcct                                    >  2:6246420/2‑93 (MQ=255)
                                                        tttATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGaa                                 >  2:4430817/1‑93 (MQ=255)
                                                        tttATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGaa                                 <  1:1206975/93‑1 (MQ=255)
                                                         ttATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAg                                >  2:3965073/1‑93 (MQ=255)
                                                         ttATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAg                                >  2:6727012/1‑93 (MQ=255)
                                                         ttATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAg                                <  2:827700/93‑1 (MQ=255)
                                                          tATTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGc                               >  2:6422647/1‑93 (MQ=255)
                                                              gAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAAc                           >  1:1742897/1‑93 (MQ=255)
                                                               aaTCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACt                          <  1:2904400/93‑1 (MQ=255)
                                                               aaTCCCGCCAGTACGAAAAACCGTTGAGCCGCGACGGTTACTCAAAATTCCTGTTTGATGAAGATGTTCGTCAAATCCTGCCTGAAGCGAACt                          <  1:3214351/93‑1 (MQ=255)
                                                                aTCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTc                         >  2:5918518/1‑93 (MQ=255)
                                                                 tCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCg                        >  1:4800017/1‑93 (MQ=255)
                                                                 tCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCg                        <  2:3050604/93‑1 (MQ=255)
                                                                 tCCCGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCg                        >  2:1290422/1‑93 (MQ=255)
                                                                   ccGCGAGTACGAAAAACCGTTGTGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCa                      <  2:1363809/93‑1 (MQ=255)
                                                                   ccGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCa                      >  2:2277797/1‑93 (MQ=255)
                                                                   ccGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCa                      <  1:803698/93‑1 (MQ=255)
                                                                   ccGCGAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCa                      >  1:779155/1‑93 (MQ=255)
                                                                    cgcgAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAg                     <  2:7758894/93‑1 (MQ=255)
                                                                     gcgAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGa                    <  2:2864536/93‑1 (MQ=255)
                                                                      cgAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGAt                   <  1:6347685/93‑1 (MQ=255)
                                                                       gAGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGAtt                  >  1:5170134/1‑93 (MQ=255)
                                                                        aGTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGAttg                 <  1:6755815/93‑1 (MQ=255)
                                                                         gTACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGAttgt                >  2:3421094/1‑93 (MQ=255)
                                                                          tACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGAttgtt               <  1:3367537/93‑1 (MQ=255)
                                                                          tACGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGAttgtt               <  2:3573891/93‑1 (MQ=255)
                                                                           aCGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGAttgttg              <  1:6084473/93‑1 (MQ=255)
                                                                           aCGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGAttgttg              >  1:1621008/1‑93 (MQ=255)
                                                                           aCGAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGAttgttg              <  2:2667507/93‑1 (MQ=255)
                                                                             gAAAAACCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGc            >  2:2585099/1‑93 (MQ=255)
                                                                              aaaaaCCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCt           <  2:6040551/93‑1 (MQ=255)
                                                                               aaaaCCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCtt          <  1:4076169/93‑1 (MQ=255)
                                                                               aaaaCCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCtt          <  2:2128085/93‑1 (MQ=255)
                                                                                 aaCCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTcc        <  1:3742807/93‑1 (MQ=255)
                                                                                 aaCCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTcc        <  1:842955/93‑1 (MQ=255)
                                                                                 aaCCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTcc        >  1:1001033/1‑93 (MQ=255)
                                                                                  aCCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCg       <  1:7315721/93‑1 (MQ=255)
                                                                                  aCCGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCg       <  2:1504379/93‑1 (MQ=255)
                                                                                   ccGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGt      <  1:5725961/93‑1 (MQ=255)
                                                                                   ccGTTGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGt      >  1:2728600/1‑93 (MQ=255)
                                                                                     gTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTa    >  1:7545966/1‑93 (MQ=255)
                                                                                      ttGAGCGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTAt   >  2:3319584/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAATCCCGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATT  >  NC_000913/3791936‑3792115

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