Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 481,225 T→G 9.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 9.9% 335.8 / 29.0 111intergenic (‑517/+29)tomB/acrBHha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (46/54);  new base G (6/5);  total (52/59)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

Only 100 of 113 total aligned reads displayed.

ATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCAT  >  NC_000913/481077‑481373
                                                                                                                                                    |                                                                                                                                                    
aTCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCt                                                                                                                                                      >  3:113713/1‑149 (MQ=255)
   tATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAc                                                                                                                                                   <  3:686725/149‑1 (MQ=255)
       cgcgcTACTTAATATAAATTAACTAACACTCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCCCTTACGCgg                                                                                                                                               <  4:931033/148‑1 (MQ=255)
       aGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCgg                                                                                                                                               >  2:309835/1‑149 (MQ=255)
         cgcTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGcc                                                                                                                                             <  2:625382/149‑1 (MQ=255)
          gcTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCt                                                                                                                                            <  4:1064937/149‑1 (MQ=255)
             aCTTAATATAAATTAACTATCAATCGGTATTATTTATGTTTTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtttt                                                                                                                                         <  1:52659/149‑3 (MQ=255)
              cTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                                                                        <  4:714809/149‑1 (MQ=255)
              cTTAATATAAATTAACTAACAATCAGTATTTTTTATGAACTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                                                                        <  3:663425/149‑1 (MQ=255)
              cTTAATATAAATTAACTAACAATCAGTATTTTTAATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                                                                        <  4:714286/149‑1 (MQ=255)
                tAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGa                                                                                                                                      <  1:54971/149‑1 (MQ=255)
                 aaTATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAt                                                                                                                                     <  4:652497/149‑1 (MQ=255)
                 aaTATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAt                                                                                                                                     <  1:601081/149‑1 (MQ=255)
                 aaTATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAt                                                                                                                                     <  2:400826/149‑1 (MQ=255)
                   tataAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTTATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTa                                                                                                                                   >  2:160969/1‑149 (MQ=255)
                   tataAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTa                                                                                                                                   >  4:1080560/1‑149 (MQ=255)
                     taAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTaca                                                                                                                                 >  2:225345/1‑149 (MQ=255)
                      aaaTTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTacac                                                                                                                                <  4:934946/149‑1 (MQ=255)
                             ctaacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTAt                                                                                                                         <  4:1071132/149‑1 (MQ=255)
                              taacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATCCTGGTTTTCGTATGATATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATc                                                                                                                        >  2:907075/1‑149 (MQ=255)
                              taacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATc                                                                                                                        >  4:175014/1‑149 (MQ=255)
                              taacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATc                                                                                                                        >  1:206331/1‑149 (MQ=255)
                               aacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCa                                                                                                                       >  2:21486/1‑149 (MQ=255)
                                 cAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAat                                                                                                                     <  1:457080/149‑1 (MQ=255)
                                   aTCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatga                                                                                                                   >  4:651691/1‑149 (MQ=255)
                                    tCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatgat                                                                                                                  >  1:774374/1‑149 (MQ=255)
                                    tCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatgat                                                                                                                  <  3:368364/149‑1 (MQ=255)
                                       gTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATgatgat                                                                                                               <  4:864987/149‑1 (MQ=255)
                                        tATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATc                                                                                                              >  3:906704/1‑149 (MQ=255)
                                          ttttttATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGACAGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGa                                                                                                            >  4:357397/1‑149 (MQ=255)
                                           tttatATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGAc                                                                                                           <  3:1065342/145‑1 (MQ=255)
                                            ttttATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACa                                                                                                          <  3:324106/149‑1 (MQ=255)
                                               tATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTa                                                                                                       <  3:102923/149‑1 (MQ=255)
                                                 tGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATg                                                                                                     >  3:863693/1‑149 (MQ=255)
                                                  gAATTTTATCCGTTGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCTGCCTTAGTATTTACACGTTGTATCAATGATGATCGACAGTATgg                                                                                                    <  4:922503/149‑1 (MQ=255)
                                                      tttATCCGTTGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgt                                                                                                <  1:390848/149‑1 (MQ=255)
                                                      tttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgt                                                                                                <  1:21486/149‑1 (MQ=255)
                                                       ttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgtg                                                                                               >  1:232008/1‑149 (MQ=255)
                                                        tATCCGTGGTTAATACTGTTTTTCGTATGAGATCCAGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTGACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGc                                                                                              <  1:318532/149‑1 (MQ=255)
                                                          tCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtttttt                                                                                                                                       >  3:44558/1‑102 (MQ=255)
                                                          tCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtttttt                                                                                                                                       <  4:44558/106‑5 (MQ=255)
                                                            cGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGa                                                                                          <  2:921320/149‑1 (MQ=255)
                                                             gTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGat                                                                                         <  4:113713/149‑1 (MQ=255)
                                                             gTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtttt                                                                                                                                         >  4:284346/1‑99 (MQ=255)
                                                             gTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtttt                                                                                                                                         <  3:284346/101‑3 (MQ=255)
                                                              tGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGata                                                                                        <  2:863489/149‑1 (MQ=255)
                                                                     tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCAt                                                                                 >  2:799415/1‑149 (MQ=255)
                                                                     tACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCttttt                                                                                                                                        >  1:917010/1‑91 (MQ=255)
                                                                            tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGt                                                                                                                              <  2:425190/97‑1 (MQ=255)
                                                                            tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGt                                                                                                                              >  1:425190/1‑97 (MQ=255)
                                                                            tttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                                                                          <  3:651691/149‑1 (MQ=255)
                                                                                  aTGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTaaa                                                                    <  2:892737/149‑1 (MQ=255)
                                                                                  aTGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTaaa                                                                    >  2:152084/1‑149 (MQ=255)
                                                                                   tGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAg                                                                   >  2:87816/1‑149 (MQ=255)
                                                                                       aTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggc                                                               >  4:1045358/1‑149 (MQ=255)
                                                                                        tCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg                                                              <  4:906704/149‑1 (MQ=255)
                                                                                         ccTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg                                                             <  2:504558/149‑1 (MQ=255)
                                                                                          cTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCaa                                                                                                                      >  1:86549/1‑91 (MQ=255)
                                                                                          cTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCaa                                                                                                                      <  2:86549/91‑1 (MQ=255)
                                                                                          cTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcggc                                                            >  3:324670/1‑149 (MQ=255)
                                                                                               ttGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacc                                                       <  3:1045358/149‑1 (MQ=255)
                                                                                                  gtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAaccacc                                                    >  4:1043059/1‑149 (MQ=255)
                                                                                                     gtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCt                                                                                                  >  1:888780/1‑100 (MQ=255)
                                                                                                     gtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCt                                                                                                  <  2:888780/100‑1 (MQ=255)
                                                                                                     gtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                 >  2:823755/1‑149 (MQ=255)
                                                                                                         aaTTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGa                                             >  4:192860/1‑149 (MQ=255)
                                                                                                         aaTTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGa                                             <  2:174781/149‑1 (MQ=255)
                                                                                                         aaTTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGa                                             >  2:753799/1‑149 (MQ=255)
                                                                                                          aTTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTc                                                                               <  1:500091/114‑1 (MQ=255)
                                                                                                          aTTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATAGCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGaa                                            <  4:758053/149‑1 (MQ=255)
                                                                                                                 cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGa                                     >  3:1077981/1‑149 (MQ=255)
                                                                                                                 cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGa                                     >  3:769270/1‑149 (MQ=255)
                                                                                                                 cTTAATGTTCGTAGGCTATGCATAAAAACGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCTTTCTTGCGGCTCAAGCGGCGGCGAACCACCACAAAGAATACCGGa                                     >  4:327691/1‑149 (MQ=255)
                                                                                                                  ttAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGc                                                                                                   >  4:684110/1‑86 (MQ=255)
                                                                                                                  ttAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGc                                                                                                   <  3:684110/86‑1 (MQ=255)
                                                                                                                  ttAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGaa                                    <  3:182866/149‑1 (MQ=255)
                                                                                                                  ttAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGaa                                    >  2:376807/1‑149 (MQ=255)
                                                                                                                   tAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAAc                                   <  4:863693/149‑1 (MQ=255)
                                                                                                                   tAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAAc                                   >  1:693342/1‑149 (MQ=255)
                                                                                                                      tGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGata                                                                                        <  1:720318/93‑1 (MQ=255)
                                                                                                                      tGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGata                                                                                        >  2:720318/1‑93 (MQ=255)
                                                                                                                      tGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAAGACCGGAACgaa                                >  2:593150/1‑149 (MQ=255)
                                                                                                                       gTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaag                               >  1:681906/1‑149 (MQ=255)
                                                                                                                         tCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaagaa                             <  1:236366/149‑1 (MQ=255)
                                                                                                                           gTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAagaaga                           >  2:905049/1‑149 (MQ=255)
                                                                                                                             aGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGAtt                         <  2:810273/149‑1 (MQ=255)
                                                                                                                               gTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGc                       >  3:597728/1‑149 (MQ=255)
                                                                                                                                  aTGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAg                    <  4:597728/149‑1 (MQ=255)
                                                                                                                                     cATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGat                                                                                         >  2:474562/1‑77 (MQ=255)
                                                                                                                                     aatcactaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGTTATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTAc                 >  1:932190/8‑149 (MQ=255)
                                                                                                                                      aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                                                                          <  4:716614/91‑1 (MQ=255)
                                                                                                                                      aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                                                                          >  3:716614/1‑91 (MQ=255)
                                                                                                                                      aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTAcc                >  3:431860/1‑149 (MQ=255)
                                                                                                                                       tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCg               >  4:164888/1‑149 (MQ=255)
                                                                                                                                         aaaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGtt             >  2:800076/1‑149 (MQ=255)
                                                                                                                                         aaaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGtt             <  3:1043059/149‑1 (MQ=255)
                                                                                                                                          ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAg                    >  4:200187/3‑141 (MQ=255)
                                                                                                                                          ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAg                    <  3:200187/139‑1 (MQ=255)
                                                                                                                                          ctaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTg            >  4:944560/3‑149 (MQ=255)
                                                                                                                                                    ttACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTAGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCAt  >  1:945467/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCAT  >  NC_000913/481077‑481373

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