Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 481,225 T→G 8.9% 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,2250TG 8.9% 441.4 / 30.9 134intergenic (‑517/+29)tomB/acrBHha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (71/51);  new base G (6/6);  total (77/57)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-02

Only 100 of 135 total aligned reads displayed.

GCCAGATGAAATTGGTGACACTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAGGTCATCAGGATCGGACGT  >  NC_000913/480955‑481488
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                       
gCCAGATGAAATTGGTGACACTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAcg                                                                                                                                                                                                                                                                     <  4:651154/275‑1 (MQ=255)
                                             gACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAAtgatga                                                                                                                                                                                                                                   >  2:6014/1‑264 (MQ=255)
                                                 taataaCCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGa                                                                                                                                                                                                                               <  4:6014/264‑1 (MQ=255)
                                                                                                                              aTTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAcg                                                                                                                                                                                                                                                                     >  6:605330/1‑149 (MQ=255)
                                                                                                                               ttAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAcgc                                                                                                                                                                                                                                                                    >  5:492527/1‑149 (MQ=255)
                                                                                                                                  gcgcTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGc                                                                                                                                                                                                                                                                 <  6:636729/149‑1 (MQ=255)
                                                                                                                                     cTACTTAATATAAATTAACTAACATTCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAGTGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCt                                                                                                                                                                                                                                                               <  7:877857/148‑1 (MQ=255)
                                                                                                                                       aCTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTa                                                                                                                                                                                                                                                             <  8:472332/148‑1 (MQ=255)
                                                                                                                                       aCTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAg                                                                                                                                                                                                                                                            <  6:9039/149‑1 (MQ=255)
                                                                                                                                         ttAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                                                                                                                                                                                           >  7:924960/1‑148 (MQ=255)
                                                                                                                                         ttAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGAATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTg                                                                                                                                                                                                                                                          >  5:726721/1‑149 (MQ=255)
                                                                                                                                           aaTATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAt                                                                                                                                                                                                                                                        <  5:559058/149‑1 (MQ=255)
                                                                                                                                             tataAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTa                                                                                                                                                                                                                                                      <  6:248936/149‑1 (MQ=255)
                                                                                                                                              ataAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTac                                                                                                                                                                                                                                                     >  6:441816/1‑149 (MQ=255)
                                                                                                                                                 aaTTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTacac                                                                                                                                                                                                                                                   >  8:851885/1‑148 (MQ=255)
                                                                                                                                                  aTTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATg                                                                                                                                                                                                                        >  2:481754/1‑174 (MQ=255)
                                                                                                                                                    taactaacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTg                                                                                                                                                                                                                                               <  5:83641/149‑1 (MQ=255)
                                                                                                                                                      actaacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtg                                                                                                                                                                                                                    <  4:481754/174‑1 (MQ=255)
                                                                                                                                                       ctaacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGt                                                                                                                                                                                                                                              <  8:940992/147‑1 (MQ=255)
                                                                                                                                                           cAATCAGTATTTTTTATGAATTTTATCGGTGGCTAATACTGGTTTTCGTAGGAGATCCTGAGTTGGTGGCGCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAATGGCTGCTTACGCGGGCTTAGTGCTTACACGGTGTATCaa                                                                                                                                                                                                                                         >  7:343350/1‑148 (MQ=255)
                                                                                                                                                                aGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTaca                                                                                                                                                                                                                                                    <  1:239380/132‑1 (MQ=255)
                                                                                                                                                                aGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTaca                                                                                                                                                                                                                                                    >  3:239380/1‑132 (MQ=255)
                                                                                                                                                                aGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAAtgatga                                                                                                                                                                                                                                   >  6:92484/1‑149 (MQ=255)
                                                                                                                                                                  tATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtt                                                                                                                                                                                                                                                              >  5:716315/1‑120 (MQ=255)
                                                                                                                                                                  tATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtt                                                                                                                                                                                                                                                              <  6:716315/120‑1 (MQ=255)
                                                                                                                                                                  tATTTTTTATGAATTTTATCCGCGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTAGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGt                                                                                                   >  2:70276/1‑275 (MQ=255)
                                                                                                                                                                     tttttATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGAc                                                                                                                                                                                                                              <  6:181197/149‑1 (MQ=255)
                                                                                                                                                                      ttttATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAAtgatga                                                                                                                                                                                                                                   >  2:737050/1‑143 (MQ=255)
                                                                                                                                                                      ttttATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGAc                                                                                                                                                                                                                              <  7:971018/148‑1 (MQ=255)
                                                                                                                                                                         tATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTa                                                                                                                                                                                                                          <  6:378553/149‑1 (MQ=255)
                                                                                                                                                                          aTGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGa                                                                                                                                                                                                                               <  4:737050/143‑1 (MQ=255)
                                                                                                                                                                               ttttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATgg                                                                                                                                                                                                                       >  6:104966/1‑146 (MQ=255)
                                                                                                                                                                               ttttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATgg                                                                                                                                                                                                                       <  5:104966/146‑1 (MQ=255)
                                                                                                                                                                                        tGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGat                                                                                                                                                                                                            >  8:266597/1‑148 (MQ=255)
                                                                                                                                                                                             aaTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTc                                                                                                                                                                                                      <  6:575059/149‑1 (MQ=255)
                                                                                                                                                                                             aaTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTACGTAGGTTATGCATAAAAACGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATACCGCCACTTAACCCGGTACCTACTGCGTTCTGCTCCCCGTAAcc                                                                        >  4:632651/1‑275 (MQ=255)
                                                                                                                                                                                               tACTGGTTTTCTTATGAGATCCTGAGTTGGTGGTTAAATTCCGCCTTAATGTTCGTCGGTTATGCATAAAAAAGGACGCTTACGGGGCCTTAGCGAATACTCGTTGTATCAATGATGATCGACAGAATGCCTATGCTCGATATCTTCa                                                                                                                                                                                                     <  8:343350/148‑1 (MQ=255)
                                                                                                                                                                                                aCTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTaa                                                                                                                                                                                                                                                                      <  8:355830/82‑2 (MQ=255)
                                                                                                                                                                                                  tGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCt                                                                                                                                                                                                 >  5:417669/1‑149 (MQ=255)
                                                                                                                                                                                                    gTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAACGAATACCGGCACGACGAAGATTGCCAGTACCGTTGCGGTCCCCCTCCCCCCCATTACACCGGTACCCACTACGTTCTGCGCCACGGCACCCGcccc                                                                  >  4:968/1‑274 (MQ=255)
                                                                                                                                                                                                    gTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCCCCAGTa                                                              >  1:59947/1‑278 (MQ=255)
                                                                                                                                                                                                    gTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCAcc                                                                  >  2:226120/1‑274 (MQ=255)
                                                                                                                                                                                                     ttttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTg                                                                                                                                                                                                                                                          >  8:248152/1‑89 (MQ=255)
                                                                                                                                                                                                     ttttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTg                                                                                                                                                                                                                                                          <  7:248152/89‑1 (MQ=255)
                                                                                                                                                                                                     ttttCGTATGAGAGCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAAcc                                                                        <  1:142565/267‑1 (MQ=255)
                                                                                                                                                                                                      tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCtt                                                                                                                                                                                                       >  7:670894/1‑139 (MQ=255)
                                                                                                                                                                                                      tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCtt                                                                                                                                                                                                       <  8:670894/139‑1 (MQ=255)
                                                                                                                                                                                                        tCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATc                                                                                                                                                                                                                                 <  2:67003/111‑1 (MQ=255)
                                                                                                                                                                                                           tATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCgg                                                                                                                                                                                            <  7:12362/145‑1 (MQ=255)
                                                                                                                                                                                                           tATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCgg                                                                                                                                                                                            >  8:12362/1‑145 (MQ=255)
                                                                                                                                                                                                            aTGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTaa                                                                                                                                                                                        <  8:859904/148‑1 (MQ=255)
                                                                                                                                                                                                                         ttGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacc                                                                                                                                                                          >  6:959086/1‑149 (MQ=255)
                                                                                                                                                                                                                         ttGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaa                                                                                                                                                   >  2:738312/1‑172 (MQ=255)
                                                                                                                                                                                                                           ggtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGATCTTCGCCCCGGAACCAGCACCAGTACTGCAAACCAGCGGCATaa                                           >  4:64566/1‑274 (MQ=255)
                                                                                                                                                                                                                            gtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATc                                                                                                                                                                                                                                 <  3:202648/91‑1 (MQ=255)
                                                                                                                                                                                                                            gtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATc                                                                                                                                                                                                                                 >  1:202648/1‑91 (MQ=255)
                                                                                                                                                                                                                             tggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGaagaag                                                                                                                                               <  4:738312/172‑1 (MQ=255)
                                                                                                                                                                                                                                 tCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatgat                                                                                                                                                                                                                                     >  4:284839/1‑82 (MQ=255)
                                                                                                                                                                                                                                   aaTTACTCCTTAATGTTCGTAGGCTATGCATAAAAAAGGCCGCTTACGCGGCCTTCGTGATTACACGTTGTATCAATGATGATTGACAGTATGGTTGTGCTTGATATCTTCATTCTTGCGTCTAAATCGGCCGCGAACCACCACAAAGAATACCGGAACGAAGCAGACTCCCATTACCGTCTCGGTCACCATCCCACCCATTAAACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTTCTGATAACCAGCGGCATAACGcccgagg                                  >  4:96286/1‑270 (MQ=255)
                                                                                                                                                                                                                                     ttACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATc                                                                                                                                                                                                                                 <  2:284839/82‑1 (MQ=255)
                                                                                                                                                                                                                                     ttACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAAt                                                                                                                                                              >  6:543420/1‑149 (MQ=255)
                                                                                                                                                                                                                                     ttACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAAt                                                                                                                                                              <  5:605330/149‑1 (MQ=255)
                                                                                                                                                                                                                                     ttACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGAtt                                                                                                                                            <  4:744284/167‑1 (MQ=255)
                                                                                                                                                                                                                                     ttACTCCTTAATGTACGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTATATCAGTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGACTCAAGCGGCGGCGAACTACCACAAAGaa                                                                                                                                                               >  8:655212/1‑148 (MQ=255)
                                                                                                                                                                                                                                        cTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAcc                                                                                                                                                           <  6:1038477/149‑1 (MQ=255)
                                                                                                                                                                                                                                          ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCg                                                                                                                                                          >  7:826447/1‑148 (MQ=255)
                                                                                                                                                                                                                                          ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTAc                                                                                                 <  3:240881/205‑1 (MQ=255)
                                                                                                                                                                                                                                          ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTAc                                                                                                 >  1:240881/1‑205 (MQ=255)
                                                                                                                                                                                                                                             tAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTaca                                                                                                                                                                                                                                                    >  1:230433/1‑55 (MQ=255)
                                                                                                                                                                                                                                             tAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTaca                                                                                                                                                                                                                                                    <  3:230433/55‑1 (MQ=255)
                                                                                                                                                                                                                                             tAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGaa                                                                                                                                                       >  8:120242/1‑148 (MQ=255)
                                                                                                                                                                                                                                                 gTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAg                    <  2:96286/275‑1 (MQ=255)
                                                                                                                                                                                                                                                    cGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTTCACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCCCCAGTACTGATAACCAGCGGCATAACCCCCAGGATAAACCCCAGCGAGGt                  >  4:108061/1‑274 (MQ=255)
                                                                                                                                                                                                                                                    cGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTACGTTCTGCTCGCCGGAACCAGCACCAGTACTGATAACAAGCGGCATAACGCCGAGGATAAACGCCCGAGAGGtcatca             >  1:125835/1‑279 (MQ=255)
                                                                                                                                                                                                                                                     gTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGat                                                                                                                                                                                                            >  1:211315/1‑87 (MQ=255)
                                                                                                                                                                                                                                                     gTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGat                                                                                                                                                                                                            <  3:211315/87‑1 (MQ=255)
                                                                                                                                                                                                                                                     gTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTTTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAGGtca                >  2:585716/1‑275 (MQ=255)
                                                                                                                                                                                                                                                         gTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGa                                                                                                                                                                                                             <  2:775694/82‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGc                                                                                                                                          >  6:732733/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCg                                                                                                                                  >  5:83752/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAAcc                                                                        >  1:128994/1‑207 (MQ=255)
                                                                                                                                                                                                                                                                 tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAAcc                                                                        <  3:128994/207‑1 (MQ=255)
                                                                                                                                                                                                                                                                   actaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaa                                                                                                                                                   <  6:757696/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aaaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTTCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCTCCGGAACCAGCCCCAGTACTGATAACCAGCGGCATAACGCCGATGATAAACGCCCACGAGGTCATCTGCATCGGACg   >  2:4091/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                    ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGtt                                                                                                                                >  7:49304/3‑148 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGaagaag                                                                                                                                               >  7:323945/1‑133 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGaagaag                                                                                                                                               <  8:323945/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTg                                                                                                                               >  6:779448/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTTCGCGCCGGACCCAGCACCAGTACTGATACCCAGCGGCCTAACGCCGAGGACAAACGCCAGCGACGTCATCAGGATCGGACGt  >  2:658666/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                     taaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAc                                                                                                                                                            >  5:944139/2‑119 (MQ=255)
                                                                                                                                                                                                                                                                     taaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAc                                                                                                                                                            <  6:944139/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                     taaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGc                                                                                                                              >  6:812740/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                     aaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTg                                                                                                                               <  7:851885/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aaGNCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATg                                                                                                                                                                                                                        >  7:931038/1‑58 (MQ=255)
                                                                                                                                                                                                                                                                      aaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATg                                                                                                                                                                                                                        <  8:931038/58‑1 (MQ=38)
                                                                                                                                                                                                                                                                          ccGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAGAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTc                                                                                                                          >  8:864229/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                            gCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaa                                                                                                                                                   >  5:211315/1‑121 (MQ=255)
                                                                                                                                                                                                                                                                            gCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCAc                                                                                                                        >  8:45983/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                            gCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCAc                                                                                                                        >  7:624857/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                              ttACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCAt                                                                                                                     >  5:214723/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                       
GCCAGATGAAATTGGTGACACTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAGGTCATCAGGATCGGACGT  >  NC_000913/480955‑481488

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