Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 481,228 C→A 15.2% 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,2280CA15.2% 332.2 / 40.8 123intergenic (‑520/+26)tomB/acrBHha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (60/43);  new base A (12/7);  total (72/51)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

Only 100 of 130 total aligned reads displayed.

TTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAA  >  NC_000913/481130‑481327
                                                                                                  |                                                                                                   
ttttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAcg                                                                                                    >  2:1244724/1‑100 (MQ=255)
 tttATCCGTGGTTAATTCTGGTTTTCGTATGAGTTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAcgc                                                                                                   <  1:415861/100‑1 (MQ=255)
   tATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCgg                                                                                                 >  1:885554/1‑100 (MQ=255)
     tCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGcc                                                                                               <  1:2857642/100‑1 (MQ=255)
        gTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTa                                                                                            >  1:2560004/1‑100 (MQ=255)
          ggTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                          <  1:1423285/100‑1 (MQ=255)
          ggTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                          >  1:22861/1‑100 (MQ=255)
           gTTAATACTGGTTTTCGTATGAGATCTTGAGTTGGTGGTTTAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                          <  1:1646605/99‑1 (MQ=255)
            ttAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                          <  1:522524/98‑1 (MQ=255)
            ttAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGa                                                                                        >  2:1379553/1‑100 (MQ=255)
            ttAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtttttta                                                                                        <  2:397229/100‑6 (MQ=255)
            ttAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtttttta                                                                                        <  1:2743964/100‑6 (MQ=255)
              aaTACTGGTTTTCGTATGAGTTTCTGAGTTGTTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttttatg                                                                                      <  1:335754/100‑8 (MQ=37)
              aaTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttttatg                                                                                      <  2:1065490/100‑8 (MQ=38)
              aaTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttttatg                                                                                      <  1:2612851/100‑8 (MQ=38)
               aTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTa                                                                                     <  1:396080/100‑1 (MQ=255)
               aTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTa                                                                                     >  2:2468690/1‑100 (MQ=255)
                tACTGGTTTTCGTATGAGATCTTGAGTTTGTGGTTTAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttttatgca                                                                                    <  1:1895185/100‑10 (MQ=37)
                tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTac                                                                                    >  2:2857703/1‑100 (MQ=255)
                   tGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACg                                                                                 >  1:2973853/1‑100 (MQ=255)
                    ggTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACg                                                                                 >  1:1453418/1‑99 (MQ=255)
                     gTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGtt                                                                               >  1:1036880/1‑100 (MQ=255)
                      ttttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTg                                                                              <  2:2880437/100‑1 (MQ=255)
                       tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGt                                                                             >  2:1655656/1‑100 (MQ=255)
                        ttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTATGTTATGCATAAAAAAAGCCACATAAGCGGCCTTAGTGATTACACGATg                                                                              >  1:2578893/1‑98 (MQ=255)
                           gTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCa                                                                         >  2:2990200/1‑100 (MQ=255)
                             aTGAGATCCTGAGTTGGTGGTTCAATTATTTCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAat                                                                       <  1:318157/100‑1 (MQ=255)
                             aTGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTCCGCGGCCTTAGTGATTATACGTTGTATCact                                                                       <  1:2043483/100‑3 (MQ=255)
                              tGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatg                                                                      >  1:634968/1‑100 (MQ=255)
                              tGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatg                                                                      <  1:1721205/100‑1 (MQ=255)
                              tGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatg                                                                      >  1:2596952/1‑100 (MQ=255)
                                  aTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCagatcgga                                                                  >  1:481141/1‑93 (MQ=255)
                                  aTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAAtgatga                                                                  >  1:1958080/1‑100 (MQ=255)
                                    ccTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATc                                                                <  1:2350785/100‑1 (MQ=255)
                                      tGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGa                                                              >  1:202304/1‑100 (MQ=255)
                                        aGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACa                                                            >  2:1308324/1‑100 (MQ=255)
                                         gTTGGTGGTTCAATTACTCCTTAATGTTCGTAGTTTATTCATAAAAAAGGCAGCTTACGCGGCATTAGTGATTACACGTTGTATCAATGATGATCGACAg                                                           <  1:30171/100‑1 (MQ=255)
                                          ttGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGt                                                          <  1:2629600/100‑1 (MQ=255)
                                            ggtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTAt                                                        >  1:772476/1‑100 (MQ=255)
                                            ggtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTAt                                                        >  1:412088/1‑100 (MQ=255)
                                            ggtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCAACAGTAt                                                        >  1:1763067/1‑100 (MQ=255)
                                             gtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATg                                                       >  2:415665/1‑100 (MQ=255)
                                              tggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATAAATGATGATCGACAGTATgg                                                      >  1:578484/1‑100 (MQ=255)
                                                 ttCAATTACTCCTTTATGATCGTAGTTTATGCATAAAAAATGCCGCTTACGCGTCTTTAGTGATTATACGTTGTATCAATTATGATCGACAGTATGGCtg                                                   <  2:1426617/100‑1 (MQ=255)
                                                  tCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgt                                                  <  1:1507941/100‑1 (MQ=255)
                                                   cAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgtg                                                 >  2:788775/1‑100 (MQ=255)
                                                   cAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgtg                                                 >  1:1770295/1‑100 (MQ=255)
                                                   cAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgtg                                                 <  2:578481/100‑1 (MQ=255)
                                                    aaTTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGc                                                <  2:2559973/100‑1 (MQ=255)
                                                      ttACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTc                                              >  1:217002/1‑100 (MQ=255)
                                                          tCCTTACTGTTCGTATGTTATGCATAAAAAATGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAAGGATGATCGACAGTATGGCTGTGCTCGata                                          <  1:824326/100‑1 (MQ=255)
                                                          tCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGata                                          <  2:2834839/100‑1 (MQ=255)
                                                           ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGatat                                         >  2:2301709/1‑100 (MQ=255)
                                                           ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGatat                                         >  1:1084004/1‑100 (MQ=255)
                                                            cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATc                                        <  2:950942/100‑1 (MQ=255)
                                                                aTGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCa                                    >  2:2750251/1‑100 (MQ=255)
                                                                aTGTTCGTAGGTTATGCATAAAAAAAGCCGCTTACGCGGCCTTCGTGATTACACGTTCTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCa                                    >  2:262129/1‑100 (MQ=255)
                                                                aTGTTCGTAGGTTATACATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTATGCTCGATATCTTCa                                    >  2:2790327/1‑100 (MQ=255)
                                                                        aGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                            <  2:2596920/100‑1 (MQ=255)
                                                                        aGGTTATACATAAAAAAGGCAGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCAATCTTGc                             >  2:3017546/1‑99 (MQ=255)
                                                                         ggTTATGCATAAAAAAGGCCGTTTATGCGGTCTTAGTGATTAAACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                            <  1:364762/99‑1 (MQ=255)
                                                                          gTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACTTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                            <  2:502262/98‑1 (MQ=255)
                                                                          gTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGc                          <  1:718496/100‑1 (MQ=255)
                                                                                cctAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGCTCGACAGTATGGCTGAGCTCGATATCTTCATTCTTGCAGCTa                        >  2:2409261/3‑96 (MQ=255)
                                                                                cATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgc                    <  1:2629631/100‑1 (MQ=255)
                                                                                cATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGACTAAAgc                    >  1:451840/1‑100 (MQ=255)
                                                                                 atcactaaGGACGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcg                   >  2:2592543/7‑100 (MQ=255)
                                                                                 aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGTCTAAAgcg                   >  2:43537/1‑100 (MQ=255)
                                                                                 aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcg                   >  1:2133205/1‑100 (MQ=255)
                                                                                 aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcg                   >  2:1741237/1‑100 (MQ=255)
                                                                                 aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTCTGCTCGATATCTTCATTCTTGCGGCTaaacc                    >  2:1269338/1‑97 (MQ=255)
                                                                                  tcactaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCAGCTAAAgcgg                  >  1:581137/6‑100 (MQ=255)
                                                                                  tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTCTCAATGATGATCTACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcgg                  >  2:578980/1‑100 (MQ=255)
                                                                                  tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcgg                  >  1:958024/1‑100 (MQ=255)
                                                                                   cactaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggc                 >  1:1246341/5‑100 (MQ=255)
                                                                                   cactaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggc                 >  1:2098225/5‑100 (MQ=255)
                                                                                    actaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg                >  1:2500610/4‑100 (MQ=255)
                                                                                    actaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTACCAATGATGATCGACAGTATAACTGTGCACGATAACTTCATACTTGCAGCTAAAgcggc                 >  1:443024/4‑99 (MQ=255)
                                                                                    actaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGAATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg                >  2:1870753/4‑100 (MQ=255)
                                                                                    aaaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggc                 >  1:1632870/1‑99 (MQ=255)
                                                                                     ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg               >  2:405257/3‑100 (MQ=255)
                                                                                     aaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg               >  1:1268090/1‑100 (MQ=255)
                                                                                     aaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg               <  1:1809685/100‑1 (MQ=255)
                                                                                      taaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcggc              >  2:439108/2‑100 (MQ=255)
                                                                                      aaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcggc              >  1:2440033/1‑100 (MQ=255)
                                                                                      aaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcggc              <  2:2795321/100‑1 (MQ=255)
                                                                                       aaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGACATCTTCATTCTTGCGGCTAAAgcggcggcg             >  1:747547/1‑100 (MQ=255)
                                                                                       aaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTAGATATCTTCATTCTTGCGGCTAAAgcggcggcg             >  2:1240489/1‑100 (MQ=255)
                                                                                         ggCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGaa           <  2:1036873/100‑1 (MQ=255)
                                                                                           ccGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacc         <  2:186349/100‑1 (MQ=255)
                                                                                            cGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacca        >  1:1596599/1‑100 (MQ=255)
                                                                                              cTTACGCGGCCTTAGTGATTACACTTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAaccacc      <  1:945834/100‑1 (MQ=255)
                                                                                              cTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAaccacc      <  2:1268078/100‑1 (MQ=255)
                                                                                              cTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAaccacc      <  2:885548/100‑1 (MQ=255)
                                                                                               ttACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAAccacca     <  2:1770275/100‑1 (MQ=255)
                                                                                               ttACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAAccacca     <  2:1860018/100‑1 (MQ=255)
                                                                                                 aCGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca   >  1:2243939/1‑100 (MQ=255)
                                                                                                  cgcgGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACaa  <  1:1909222/100‑1 (MQ=255)
                                                                                                  cgcgGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACaa  >  2:2066644/1‑100 (MQ=255)
                                                                                                  cgcgGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTCTGCTCGATAGCTTCATTCTTGCGGCAAACGCGGCGGCGAACCACCCCaa  >  2:1007877/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAA  >  NC_000913/481130‑481327

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: