Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 481,225 T→G 18.5% intergenic (‑517/+29) tomB ← / ← acrB Hha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913481,2250TG18.5% 420.4 / 79.2 176intergenic (‑517/+29)tomB/acrBHha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (68/75);  new base G (19/14);  total (87/89)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

GAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACC  >  NC_000913/481127‑481323
                                                                                                  |                                                                                                  
gAATTTTATCCGTGGTTAATACTTGTTTTCGTATGATATCCTGAGTTTGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCtt                                                                                                   <  1:3410342/100‑1 (MQ=255)
 aaTTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTa                                                                                                  >  2:387750/1‑100 (MQ=255)
  aTTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAc                                                                                                 <  2:48129/100‑1 (MQ=255)
   ttttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAcg                                                                                                <  2:2473729/100‑1 (MQ=255)
     ttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAcgcg                                                                                              <  1:2053536/100‑1 (MQ=255)
     ttATCAGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTCCAAGTACTCCTTAATGCTCGTAGGTTATGCATAAAAAAGGCCGCTTAcgcg                                                                                              <  1:2810973/100‑1 (MQ=255)
         ccGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGcc                                                                                           <  2:2924837/99‑1 (MQ=255)
          cGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCtt                                                                                         <  2:3528510/100‑1 (MQ=255)
          cGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCtt                                                                                         <  1:2249829/100‑1 (MQ=255)
          cGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCtt                                                                                         >  1:1868812/1‑100 (MQ=255)
          cGTGGTTAATACTGGATTTCGTATGGGAGCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtt                                                                                         <  1:687561/100‑1 (MQ=255)
           gTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCtt                                                                                         >  1:3696038/1‑99 (MQ=255)
           gTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttt                                                                                        <  1:1016569/100‑2 (MQ=255)
           gTGGTTAATACTGGTTTTCGTATGAGAGCCTGAGTTGGTGGTTAAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGcc                                                                                           <  1:2641166/97‑1 (MQ=255)
            tGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAg                                                                                       <  2:84815/100‑1 (MQ=255)
             ggTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                      >  2:3690521/1‑100 (MQ=255)
              gTTAATACTGGTTTTCGTATGAGATCCTGATTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTg                                                                                     <  2:313512/100‑1 (MQ=255)
              gTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGTTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttt                                                                                      <  1:2073242/99‑4 (MQ=255)
              gTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                      <  1:3606751/99‑1 (MQ=255)
               ttAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtttttta                                                                                    <  2:223734/100‑6 (MQ=255)
                tAATATTGGTTTTCGTATGAGATCTTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAt                                                                                   <  1:872233/100‑1 (MQ=255)
                 aaTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAtt                                                                                  <  1:1190156/100‑1 (MQ=255)
                 aaTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttttatg                                                                                  <  2:101478/100‑8 (MQ=38)
                   tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTac                                                                                >  1:1832556/1‑100 (MQ=255)
                   tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTac                                                                                >  2:1521811/1‑100 (MQ=255)
                   tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTac                                                                                <  1:2158575/100‑1 (MQ=255)
                   tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTac                                                                                <  2:1731168/100‑1 (MQ=255)
                    aCTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTaca                                                                               <  1:1960042/100‑1 (MQ=255)
                     cTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACCCGGCCTTAGTGATTacac                                                                              >  1:3143173/1‑100 (MQ=255)
                      tGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACg                                                                             <  2:2317020/100‑1 (MQ=255)
                         ttttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCAGAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTg                                                                          <  2:3279950/100‑1 (MQ=255)
                          tttCGTATGAGATCCTGAGTTGGTTGTTCAATTACTCCCTATTGTTCGTAGGTTATGCATAAAAAAGGTCGCGTAAGCGGCCttttttat                                                                                   <  1:3466757/90‑7 (MQ=255)
                          tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGt                                                                         >  2:3383931/1‑100 (MQ=255)
                          tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGt                                                                         <  1:1714982/100‑1 (MQ=255)
                              gTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCa                                                                     >  1:1614136/1‑100 (MQ=255)
                              gTATGAGATCCTGAGATGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGCTTACACGTTGTATCa                                                                     >  1:3057143/1‑100 (MQ=255)
                               tATGAGATCCTGAGTTGGTTGTTCAATTACTCCTTAATGTTCGTAGTTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCaa                                                                    <  1:3139807/100‑1 (MQ=255)
                                aTGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAat                                                                   <  2:3041241/100‑1 (MQ=255)
                                aTGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAat                                                                   >  1:1614022/1‑100 (MQ=255)
                                 tGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGATCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatg                                                                  >  1:2035898/1‑100 (MQ=255)
                                      tCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATgatgat                                                             <  1:541466/100‑1 (MQ=255)
                                        cTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCg                                                           >  1:2977519/1‑100 (MQ=255)
                                              tggtggTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTa                                                     <  1:753975/100‑1 (MQ=255)
                                              tggtggTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTa                                                     <  1:1781639/100‑1 (MQ=255)
                                               ggtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTAt                                                    >  2:3745740/1‑100 (MQ=255)
                                               ggtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTAt                                                    >  2:1715644/1‑100 (MQ=255)
                                                  ggtTCAATTACTCGTTAATGTTCGTAGGTTATGCATAAAAAAGGCCCCTTACGCGGCCTTAGTGATTACACAGTGTATCAATGATGAACTACAGTATGGc                                                 >  1:2112685/1‑100 (MQ=255)
                                                  ggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGc                                                 <  1:1242285/100‑1 (MQ=255)
                                                  ggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGc                                                 >  1:2238855/1‑100 (MQ=255)
                                                    ttCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtg                                               <  2:203443/100‑1 (MQ=255)
                                                      cAATTACTCCTTAATGTTCGTAGTTTATGCATAAAAAAGGCCGCTTACACGGCCTTAGTGATTAAACGTTGTATCAATAATGAGCGACAGTATGGCtgtg                                             >  1:1166308/1‑100 (MQ=255)
                                                        aTTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCt                                           >  2:1622545/1‑100 (MQ=255)
                                                           aCTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGa                                        <  1:3757504/100‑1 (MQ=255)
                                                             tCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGata                                      <  1:2782392/100‑1 (MQ=255)
                                                              ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGatat                                     >  2:964386/1‑100 (MQ=255)
                                                              ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGatat                                     >  1:2824269/1‑100 (MQ=255)
                                                              ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTATGCTCGatat                                     >  1:1216431/1‑100 (MQ=255)
                                                               cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATc                                    >  1:1415604/1‑100 (MQ=255)
                                                                 tAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCtt                                  <  1:341693/100‑1 (MQ=255)
                                                                 tAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCtt                                  <  2:3533988/100‑1 (MQ=255)
                                                                     gTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCAtt                              >  2:2981538/1‑100 (MQ=255)
                                                                     gTTCGGAGGTTATGCATAAAAAAGGCCGCTTACTCGGCCTGATTGTTTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCAt                               <  1:3612478/99‑1 (MQ=255)
                                                                         gcaGGTTATGCATTAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTg                          >  2:2591269/3‑100 (MQ=255)
                                                                          tAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGc                         >  1:2403328/1‑100 (MQ=255)
                                                                           aGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                        <  2:983564/100‑1 (MQ=255)
                                                                           aGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                        >  1:2298150/1‑100 (MQ=255)
                                                                           aGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                        <  1:3383911/100‑1 (MQ=255)
                                                                             gTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGc                      >  2:510631/1‑100 (MQ=255)
                                                                               tATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCt                     >  2:2815787/1‑99 (MQ=255)
                                                                               tATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTa                    <  1:3558668/100‑1 (MQ=255)
                                                                                   cATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgc                >  1:253681/1‑100 (MQ=255)
                                                                                   cATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgc                >  1:3253209/1‑100 (MQ=255)
                                                                                   cATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgc                >  1:247163/1‑100 (MQ=255)
                                                                                   cATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgc                >  1:3726219/1‑100 (MQ=255)
                                                                                    aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcg               >  1:2275069/1‑100 (MQ=255)
                                                                                     tcactaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcgg              >  1:1214683/6‑100 (MQ=255)
                                                                                     tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCCATGATGATCGACAGTATGGCTGTGCTCGATCTCTTCATTCTTGCGGCTAAAgcgg              >  1:819315/1‑100 (MQ=255)
                                                                                     tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcgg              >  2:3082294/1‑100 (MQ=255)
                                                                                     tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcgg              >  2:3337197/1‑100 (MQ=255)
                                                                                      cactaaGGCGGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGCCAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggc             >  1:3604110/5‑100 (MQ=255)
                                                                                      cactaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggc             >  2:2180128/5‑100 (MQ=255)
                                                                                       actaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg            >  1:2314175/4‑100 (MQ=255)
                                                                                       actaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg            >  1:657903/4‑100 (MQ=255)
                                                                                       aaaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg            <  2:115260/100‑1 (MQ=255)
                                                                                       aaaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg            >  1:115258/1‑100 (MQ=255)
                                                                                       aaaaaGGCCGCTTACGCGGACTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTATGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg            >  1:874288/1‑100 (MQ=255)
                                                                                        ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg           >  1:2478875/3‑100 (MQ=255)
                                                                                        ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg           >  2:2316995/3‑100 (MQ=255)
                                                                                        ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg           >  1:2387621/3‑100 (MQ=255)
                                                                                        ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg           >  1:2146105/3‑100 (MQ=255)
                                                                                        ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGATGTCTCAATGATGATCGACAGTCTGGCTGTGCTCGATATATTCATTCTTGCCGCTAAAgc                >  1:3343145/3‑95 (MQ=255)
                                                                                        aaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg           <  1:498860/100‑1 (MQ=255)
                                                                                         taaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggaag           >  1:199950/2‑96 (MQ=255)
                                                                                          aaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcggcg         >  1:3793372/1‑100 (MQ=255)
                                                                                          aaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcggc          >  1:423246/1‑99 (MQ=255)
                                                                                            ggCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGaa       >  1:1919240/1‑100 (MQ=255)
                                                                                            ggCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGaa       <  2:1565617/100‑1 (MQ=255)
                                                                                                gCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAaccac   <  1:1351069/100‑1 (MQ=255)
                                                                                                gCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAaccac   >  2:1496575/1‑100 (MQ=255)
                                                                                                 cTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAaccacc  <  1:950995/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
GAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACC  >  NC_000913/481127‑481323

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