Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,680,105 T→A 6.5% intergenic (+157/+375) USA300HOU_2536 → / ← USA300HOU_2537 possible secretory antigen/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,680,1050TA 6.5% 177.1 / 2.4 61intergenic (+157/+375)USA300HOU_2536/USA300HOU_2537possible secretory antigen/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/26);  new base A (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GAACGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGA  >  CP000730/2679967‑2680245
                                                                                                                                          |                                                                                                                                               
gAACGTTTCAATTTCCGTTATGCGTGTGATAGAATTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCtatat                                                                                                                                               <  1:347123/141‑1 (MQ=255)
  aCGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCGTAATCtatat‑at                                                                                                                                            >  2:624711/1‑141 (MQ=255)
  aCGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCtatat‑at                                                                                                                                            <  2:329884/141‑1 (MQ=255)
  aCGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCtatat‑at                                                                                                                                            <  2:712870/141‑1 (MQ=255)
   cGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTAGAAATTTCAAATTGCCGTAATTGAACGGATATTTCTTCTTCAACTATTATTTCATCTTTAGCGTAGTCtatat‑att                                                                                                                                           >  1:83034/1‑140 (MQ=255)
    gCTTCAATTTCCTGTATGCGTGTGATAGAACTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATTTTTAGCATTATCTATAT‑ATaa                                                                                                                                          <  1:399064/141‑1 (MQ=255)
     cTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTGCAAATTTCAAATTGCGGTAATTGAACGTATGTTTCTTCTTCAACTATTATTTCATATTTCGCATAATCTATAT‑ATaaa                                                                                                                                         >  2:249027/1‑141 (MQ=39)
     cTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATaaa                                                                                                                                         >  1:746177/1‑141 (MQ=255)
     cTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATaaa                                                                                                                                         >  1:232470/1‑141 (MQ=255)
     cTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATaaa                                                                                                                                         >  2:470687/1‑141 (MQ=255)
             tCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATg                                                                                                                                 >  2:99653/1‑141 (MQ=255)
             tCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATg                                                                                                                                 <  1:513852/141‑1 (MQ=255)
             tCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ACAAAATTTTATg                                                                                                                                 >  1:211850/1‑141 (MQ=255)
               cTGTATGCGTGTGAAAGAAGTTTTTTTTTTTTTTTATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAAtttttttct                                                                                                                               >  2:521445/1‑136 (MQ=255)
                   aTGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTTATCTTTAGCATAATCTATAT‑ATACAATTTTATGCtatta                                                                                                                           <  2:291506/141‑1 (MQ=255)
                      cGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTAttt                                                                                                                        <  1:68143/141‑1 (MQ=255)
                              aGAAGGTTTTATTTTATGACATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAAtt                                                                                                                <  2:379506/141‑1 (MQ=255)
                                  gTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCt                                                                                                            <  2:465879/141‑1 (MQ=255)
                                   tttttATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCta                                                                                                           >  1:687252/1‑141 (MQ=255)
                                             atGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACa                                                                                                 <  2:121930/141‑1 (MQ=255)
                                              tGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACAt                                                                                                >  1:447300/1‑141 (MQ=255)
                                               gAAATTATATTATTACTTCAACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATa                                                                                               <  2:691720/141‑1 (MQ=255)
                                                      tattatTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTc                                                                                        <  1:694148/141‑1 (MQ=255)
                                                          attaCTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGata                                                                                    >  2:516597/1‑141 (MQ=255)
                                                          actaCTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGata                                                                                    <  1:227235/139‑1 (MQ=255)
                                                            taCTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGAtata                                                                                  <  1:695709/141‑1 (MQ=255)
                                                               ttCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAAc                                                                               >  1:642193/1‑141 (MQ=255)
                                                                 cTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ACAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCt                                                                             <  2:211850/141‑1 (MQ=255)
                                                                   aCAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATCAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATCAACCtgg                                                                           >  2:22013/1‑139 (MQ=255)
                                                                     aaaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTCTTTGAACTTTAGCACAAGCTATAT‑ATTAACTTTTATGCTACTATTTAATTGATTAGCTTTAACTTAACATAGGTTTTCGATATAACACTTGtt                                                                         >  2:278381/1‑141 (MQ=255)
                                                                     aaaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATAGTTTTTCGATATAAACATTGtt                                                                         >  1:303610/1‑141 (MQ=255)
                                                                      aaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACTTTTTCGATATAAACCTTTTTc                                                                        >  2:33779/1‑141 (MQ=255)
                                                                      aaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTc                                                                        <  2:676199/141‑1 (MQ=255)
                                                                       aTTTCAAATTGCCGTAATTGAACGTATATTTCCTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCt                                                                       <  1:39866/141‑1 (MQ=255)
                                                                        tttCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTTTTCTa                                                                      >  1:702189/1‑141 (MQ=255)
                                                                        tttCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTa                                                                      >  1:304389/1‑141 (MQ=255)
                                                                        tttCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTa                                                                      >  1:502967/1‑141 (MQ=255)
                                                                             aaTTGCCGTAATTGAACGTATATTTCTTCTACAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATATGTTTTCGATATAAAACTTTTTTTAAAtct                                                                 >  1:183472/1‑141 (MQ=255)
                                                                              aTTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAATTTAACATACGTTTTCGATATAAACCTTGTTCTAActctt                                                                >  2:457899/1‑140 (MQ=255)
                                                                                tGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑aa                                                             >  1:443802/1‑141 (MQ=255)
                                                                                    gTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACATTGTTCTAAATCTC‑aataat                                                         >  2:489994/1‑141 (MQ=255)
                                                                                            aCGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑Ct                                                <  2:232470/141‑1 (MQ=255)
                                                                                            aCGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑Ct                                                <  2:746177/141‑1 (MQ=255)
                                                                                             cGTATATTTCTTCTTCAACTATTATTTCATATTTAGTATAGTATCTTTGTTATCACGTTCATC‑AGTGTGGAAATAATTCGCTATCCTCTGTCATACTTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTg                                               <  2:350155/141‑1 (MQ=255)
                                                                                              gTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGt                                              <  1:470687/141‑1 (MQ=255)
                                                                                              gTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGt                                              <  1:461613/141‑1 (MQ=255)
                                                                                                      cttcttcAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTCTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCg                                      <  1:213309/141‑1 (MQ=255)
                                                                                                      cttcttcAACTATTAGTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTCTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCg                                      <  1:213139/141‑1 (MQ=255)
                                                                                                        tcttcAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑ACTAATTTTTTG‑CTGTTTTCATCGTc                                    >  2:760860/1‑141 (MQ=255)
                                                                                                         cttcAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCa                                   <  2:671814/141‑1 (MQ=255)
                                                                                                         cttcAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCa                                   <  2:131742/141‑1 (MQ=255)
                                                                                                          ttcAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAGTAATTCGCTATAACTTAAAATACGTTTTCGATATAAACCTTGTTCTAAATATC‑AATAATTTTTTG‑CTGTTTTGATTGTCAt                                  >  1:254246/1‑141 (MQ=255)
                                                                                                            cAACTATTATTTCATCTTTAGCATAATCTATAT‑GTAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGttatta                                <  1:521445/141‑1 (MQ=255)
                                                                                                              aCTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCATTAGt                              >  1:714437/1‑141 (MQ=255)
                                                                                                                  ttattTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCTTCGTTATTAGTTaaa                          >  2:115200/1‑141 (MQ=255)
                                                                                                                        catcttcggCATAGTCTAAGT‑GTAAAATTTTATGCTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAATCACTTACAATTTCTTGATTTTTTTCGTCGTCATCAGTTaaaaagt                      >  2:277158/9‑139 (MQ=255)
                                                                                                                        catcttcgacATAGTCTAAGT‑ATAAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACGAACCTTTTTCTAAATCACTTACAATTTCTTGATTTTTTTCATCATCATCCGTTaaaaaat                      >  1:397679/10‑140 (MQ=255)
                                                                                                                        cATCTTCGGCATCGTCTAAGT‑TTAAAATTTTATGCTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAATCACTTACAATTTCTTGATTTTTTTCGTCGTCATCAGTTaaaaagt                      <  1:468003/141‑3 (MQ=255)
                                                                                                                         atcgtcggCATAGTCTAAGT‑TTAAAATTTTATGCTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAATCACTTACAATTTTTTGATTTTTTTCGTCGTCATCAGTTaaaaagtc                     <  1:219701/134‑4 (MQ=255)
                                                                                                                              tAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTTTAAATCTC‑AATAATTTTTTG‑CTGTTTTTATCGTCATTAGTTAAAAAAATAATTTat              >  2:564485/1‑140 (MQ=255)
                                                                                                                                gCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACt            >  2:546254/1‑141 (MQ=255)
                                                                                                                                gCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTAGATATAAACCTTGTTATAAATCTC‑ACTAATTTTTTG‑CTGTTTTCATCGGCATTAGTTAAAAAAATAATTTCACt            >  2:149684/1‑141 (MQ=255)
                                                                                                                                  aTAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATACGTTCACATGCGTTTTTGATATAAACCGTGTTCTAAATCTC‑AAAAATTTTTTT‑CTGTTTTCATCGTCCTTAGTTATAAAAATAATTTAATTgg          >  1:368957/1‑140 (MQ=255)
                                                                                                                                  aTAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGa          <  1:188527/141‑1 (MQ=255)
                                                                                                                                     aTCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCACTATAATTTAAACTCCATTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATCTTTTG‑CTGTTTTTATCTTCATTAGTTAAAAAAATCAtttcaatgtctt       >  2:446534/1‑132 (MQ=255)
                                                                                                                                      tCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTT‑AATAACTATTTT‑CTGGTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGttt      <  2:304389/141‑1 (MQ=255)
                                                                                                                                         ttat‑ataAATTTTTATGCTATTATTTAATTAATTCGCTATAACTTAACATACTTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTg   <  2:59289/140‑1 (MQ=255)
                                                                                                                                         atat‑ataAAATTTTATGGTGTTATTTGAGTAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTg   >  1:185125/1‑141 (MQ=255)
                                                                                                                                         atat‑ataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATCCGTTTTCGATATAAACCTTGTTCTAAATCTC‑CATCAACTTTTG‑CTGCGTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTg   <  2:735278/141‑1 (MQ=255)
                                                                                                                                         atat‑ataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTTGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGttttta   >  2:141951/1‑140 (MQ=255)
                                                                                                                                          tat‑ataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTga  <  1:624711/141‑1 (MQ=255)
                                                                                                                                          tat‑ataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTga  <  2:687252/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                               
GAACGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATAT‑ATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑CTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGA  >  CP000730/2679967‑2680245

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