Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 481,228 C→A 18.5% intergenic (‑520/+26) 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,2280CA18.5% 453.7 / 72.8 175intergenic (‑520/+26)tomB/acrBHha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (69/73);  new base A (19/14);  total (88/87)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

ATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACA  >  NC_000913/481129‑481326
                                                                                                   |                                                                                                  
aTTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAc                                                                                                    <  2:48129/100‑1 (MQ=255)
   ttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAcgcg                                                                                                 >  1:3730392/1‑100 (MQ=255)
   ttATCAGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTCCAAGTACTCCTTAATGCTCGTAGGTTATGCATAAAAAAGGCCGCTTAcgcg                                                                                                 <  1:2810973/100‑1 (MQ=255)
        cGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCtt                                                                                            >  2:336048/1‑100 (MQ=255)
         gTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCtt                                                                                            >  1:3696038/1‑99 (MQ=255)
         gTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttt                                                                                           <  1:1016569/100‑2 (MQ=255)
          tGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAg                                                                                          <  2:84815/100‑1 (MQ=255)
           ggggaataatGGTTTTCTTATGAGATCATGAATAAGAGGTTAAATTACTCCTTTTTGTTCGAAAAAAAAAAAAAAAAAAGGCCGCTTACGCGGCCTTAg                                                                                          <  1:2628417/95‑1 (MQ=255)
           ggTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttt                                                                                         <  2:354868/100‑4 (MQ=255)
            gTTAATACTGGTTTTCGTATGAGATCCTGATTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTg                                                                                        <  2:313512/100‑1 (MQ=255)
            gTTAATACTGGTATTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTg                                                                                        >  2:2505339/1‑100 (MQ=255)
            gTTAATACTAGTTTTCGTATGAGTTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGTTTATGCATAAAAAAGGCCGCGTAAGCGGCCtttttt                                                                                        <  1:3709502/100‑5 (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)
               aaTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttttatg                                                                                     <  2:2234409/100‑8 (MQ=38)
               aaTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttttatg                                                                                     <  1:2189853/100‑8 (MQ=38)
                 tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTac                                                                                   >  2:1521811/1‑100 (MQ=255)
                 tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTac                                                                                   <  2:1731168/100‑1 (MQ=255)
                 tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTac                                                                                   >  1:1832556/1‑100 (MQ=255)
                 tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttttatgca                                                                                   <  1:2257359/100‑10 (MQ=38)
                  aCTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTaca                                                                                  <  1:1960042/100‑1 (MQ=255)
                   cTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACCCGGCCTTAGTGATTacac                                                                                 >  1:3143173/1‑100 (MQ=255)
                    tGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACg                                                                                <  1:301118/100‑1 (MQ=255)
                     ggTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGt                                                                               <  1:424176/100‑1 (MQ=255)
                     ggTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGt                                                                               <  1:1940396/100‑1 (MQ=255)
                        tttCGTATGAGATCCTGAGTTGGTTGTTCAATTACTCCCTATTGTTCGTAGGTTATGCATAAAAAAGGTCGCGTAAGCGGCCttttttat                                                                                      <  1:3466757/90‑7 (MQ=255)
                        tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGt                                                                            <  1:1714982/100‑1 (MQ=255)
                        tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGt                                                                            >  2:341698/1‑100 (MQ=255)
                        tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGt                                                                            <  2:1611819/100‑1 (MQ=255)
                           cGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATc                                                                         <  2:2178181/100‑1 (MQ=255)
                            gTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATc                                                                         >  1:3347430/1‑99 (MQ=255)
                             tATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCaa                                                                       >  1:3259126/1‑100 (MQ=255)
                              aTGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAat                                                                      >  1:1614022/1‑100 (MQ=255)
                               tGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGATCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatg                                                                     >  1:2035898/1‑100 (MQ=255)
                                gagaTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTAGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatga                                                                    <  1:2069418/100‑1 (MQ=255)
                                 agaTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatgat                                                                   >  1:2732739/1‑100 (MQ=255)
                                    tCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATgatgat                                                                <  1:541466/100‑1 (MQ=255)
                                     ccTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATc                                                               <  2:3001628/100‑1 (MQ=255)
                                     ccTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATc                                                               >  1:2410169/1‑100 (MQ=255)
                                      cTGTGTTGGTAGTTCAATTAGTCCTTACTGTTCGTGGGTAATGCTTAAAAATGGCCGCTTGCGCGGGCTTAGTGATTACACGTTGTATCAATgatgattg                                                              <  1:2242775/100‑3 (MQ=255)
                                      cTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCAGAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCg                                                              >  2:1178525/1‑100 (MQ=255)
                                         aGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACa                                                           <  2:2410181/100‑1 (MQ=255)
                                            tggtggTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTa                                                        <  1:1781639/100‑1 (MQ=255)
                                              gtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATg                                                      >  2:3069172/1‑100 (MQ=255)
                                                ggtTCAATTACTCGTTAATGTTCGTAGGTTATGCATAAAAAAGGCCCCTTACGCGGCCTTAGTGATTACACAGTGTATCAATGATGAACTACAGTATGGc                                                    >  1:2112685/1‑100 (MQ=255)
                                                ggtTCAATTACTCCTTAATGTTCGTAGGTTATGCCTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGc                                                    <  1:3237273/100‑1 (MQ=255)
                                                ggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGc                                                    <  1:1242285/100‑1 (MQ=255)
                                                ggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGc                                                    <  2:3780538/100‑1 (MQ=255)
                                                  ttCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtg                                                  <  2:203443/100‑1 (MQ=255)
                                                   tCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgt                                                 <  1:2098582/100‑1 (MQ=255)
                                                    cAATTACTCCTTAATGTTCGTAGTTTATGCATAAAAAAGGCCGCTTACACGGCCTTAGTGATTAAACGTTGTATCAATAATGAGCGACAGTATGGCtgtg                                                >  1:1166308/1‑100 (MQ=255)
                                                    cAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgt                                                 >  2:3306610/1‑99 (MQ=255)
                                                        tACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCg                                            >  2:58311/1‑100 (MQ=255)
                                                          cTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGat                                          <  2:3127826/100‑1 (MQ=255)
                                                           tCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGata                                         <  1:2782392/100‑1 (MQ=255)
                                                           tCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATAGACAGTATGGCTGTGCGCGata                                         <  1:3581777/100‑1 (MQ=255)
                                                            ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGatat                                        >  2:964386/1‑100 (MQ=255)
                                                            ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGatat                                        >  1:2392863/1‑100 (MQ=255)
                                                            ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGata                                         >  1:1216558/1‑99 (MQ=255)
                                                             cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATc                                       >  1:1415604/1‑100 (MQ=255)
                                                               tAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCtt                                     <  2:3533988/100‑1 (MQ=255)
                                                               tAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCtt                                     <  1:341693/100‑1 (MQ=255)
                                                                aaTGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTc                                    >  1:3482539/1‑100 (MQ=255)
                                                                   gTTCGGAGGTTATGCATAAAAAAGGCCGCTTACTCGGCCTGATTGTTTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCAt                                  <  1:3612478/99‑1 (MQ=255)
                                                                         aGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                           <  2:983564/100‑1 (MQ=255)
                                                                         aGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                           <  2:3730416/100‑1 (MQ=255)
                                                                         aGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                           >  1:2298150/1‑100 (MQ=255)
                                                                           gTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGc                         <  1:1676791/100‑1 (MQ=255)
                                                                             tATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCt                        >  2:2815787/1‑99 (MQ=255)
                                                                                 cATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgc                   >  1:247163/1‑100 (MQ=255)
                                                                                 cATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgc                   >  2:2696458/1‑100 (MQ=255)
                                                                                 aatccctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGAGGATCGACAGTATGGGTGTGCAGGATATATTGATTCTTGCGGATa                       >  1:890440/8‑96 (MQ=255)
                                                                                  aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcg                  >  1:2275069/1‑100 (MQ=255)
                                                                                   tcactaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcgg                 >  1:1214683/6‑100 (MQ=255)
                                                                                   tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCCATGATGATCGACAGTATGGCTGTGCTCGATCTCTTCATTCTTGCGGCTAAAgcgg                 >  1:819315/1‑100 (MQ=255)
                                                                                    aaaaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggc                >  2:1687579/1‑100 (MQ=255)
                                                                                     actaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg               >  1:750474/4‑100 (MQ=255)
                                                                                     actaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg               >  1:657903/4‑100 (MQ=255)
                                                                                     actaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg               >  1:2933309/4‑100 (MQ=255)
                                                                                     aaaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg               >  1:115258/1‑100 (MQ=255)
                                                                                      ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg              >  1:2146105/3‑100 (MQ=255)
                                                                                      ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg              >  2:2316995/3‑100 (MQ=255)
                                                                                      aaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg              <  1:498860/100‑1 (MQ=255)
                                                                                       taaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggaag              >  1:199950/2‑96 (MQ=255)
                                                                                       aaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcggc             <  2:2948173/100‑1 (MQ=255)
                                                                                        aaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcggcg            >  1:3793372/1‑100 (MQ=255)
                                                                                        aaGGCCGCGTAAGCGGCCTAAGCGATTACACGTTGTATCAATGATCAGCCACAGTACGCCTGTGCGCGATATCTTCATTCATGCTGCTAAAg                    >  1:3605207/1‑92 (MQ=255)
                                                                                         aGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGa           <  1:3108069/100‑1 (MQ=255)
                                                                                          ggCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGaa          <  2:1565617/100‑1 (MQ=255)
                                                                                           gCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAac         <  1:1914924/100‑1 (MQ=255)
                                                                                            ccGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacc        >  2:1133003/1‑100 (MQ=255)
                                                                                             cccTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGGTCGATATCTTCATTCTTGCGGCTATAGCGGCGGCGAac         <  1:2359335/96‑1 (MQ=255)
                                                                                             cGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacca       <  1:806823/100‑1 (MQ=255)
                                                                                              gCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTGATTCTTGCGGCTAAAGCGGCGGCGAaccac      <  1:2505327/100‑1 (MQ=255)
                                                                                              gCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAaccac      <  1:1351069/100‑1 (MQ=255)
                                                                                                 tACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACcaccac   >  1:968207/1‑100 (MQ=255)
                                                                                                 tACGCGGCCTTAGGGATTACACGTTGTATCACTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCGAAAGCGGCGGCGAACcaccac   >  1:421492/1‑100 (MQ=255)
                                                                                                  aCGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca  <  2:630863/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
ATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACA  >  NC_000913/481129‑481326

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: