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

GAACGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGA  >  CP000730/2679967‑2680245
                                                                                                                                          |                                                                                                                                               
gAACGTTTCAATTTCCGTTATGCGTGTGATAGAATTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCtatat                                                                                                                                               <  1:347123/141‑1 (MQ=255)
  aCGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCGTAATCtatatat                                                                                                                                            >  2:624711/1‑141 (MQ=255)
  aCGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCtatatat                                                                                                                                            <  2:329884/141‑1 (MQ=255)
  aCGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCtatatat                                                                                                                                            <  2:712870/141‑1 (MQ=255)
   cGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTAGAAATTTCAAATTGCCGTAATTGAACGGATATTTCTTCTTCAACTATTATTTCATCTTTAGCGTAGTCtatatatt                                                                                                                                           >  1:83034/1‑140 (MQ=255)
    gCTTCAATTTCCTGTATGCGTGTGATAGAACTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATTTTTAGCATTATCTATATATaa                                                                                                                                          <  1:399064/141‑1 (MQ=255)
     cTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTGCAAATTTCAAATTGCGGTAATTGAACGTATGTTTCTTCTTCAACTATTATTTCATATTTCGCATAATCTATATATaaa                                                                                                                                         >  2:249027/1‑141 (MQ=39)
     cTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATaaa                                                                                                                                         >  1:746177/1‑141 (MQ=255)
     cTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATaaa                                                                                                                                         >  1:232470/1‑141 (MQ=255)
     cTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATaaa                                                                                                                                         >  2:470687/1‑141 (MQ=255)
             tCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATg                                                                                                                                 >  2:99653/1‑141 (MQ=255)
             tCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATg                                                                                                                                 <  1:513852/141‑1 (MQ=255)
             tCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATACAAAATTTTATg                                                                                                                                 >  1:211850/1‑141 (MQ=255)
               cTGTATGCGTGTGAAAGAAGTTTTTTTTTTTTTTTATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAAtttttttct                                                                                                                               >  2:521445/1‑136 (MQ=255)
                   aTGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTTATCTTTAGCATAATCTATATATACAATTTTATGCtatta                                                                                                                           <  2:291506/141‑1 (MQ=255)
                      cGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTAttt                                                                                                                        <  1:68143/141‑1 (MQ=255)
                              aGAAGGTTTTATTTTATGACATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAAtt                                                                                                                <  2:379506/141‑1 (MQ=255)
                                  gTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCt                                                                                                            <  2:465879/141‑1 (MQ=255)
                                   tttttATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCta                                                                                                           >  1:687252/1‑141 (MQ=255)
                                             atGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACa                                                                                                 <  2:121930/141‑1 (MQ=255)
                                              tGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACAt                                                                                                >  1:447300/1‑141 (MQ=255)
                                               gAAATTATATTATTACTTCAACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATa                                                                                               <  2:691720/141‑1 (MQ=255)
                                                      tattatTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTc                                                                                        <  1:694148/141‑1 (MQ=255)
                                                          attaCTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGata                                                                                    >  2:516597/1‑141 (MQ=255)
                                                          actaCTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGata                                                                                    <  1:227235/139‑1 (MQ=255)
                                                            taCTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGAtata                                                                                  <  1:695709/141‑1 (MQ=255)
                                                               ttCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAAc                                                                               >  1:642193/1‑141 (MQ=255)
                                                                 cTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATACAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCt                                                                             <  2:211850/141‑1 (MQ=255)
                                                                   aCAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATCAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATCAACCtgg                                                                           >  2:22013/1‑139 (MQ=255)
                                                                     aaaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTCTTTGAACTTTAGCACAAGCTATATATTAACTTTTATGCTACTATTTAATTGATTAGCTTTAACTTAACATAGGTTTTCGATATAACACTTGtt                                                                         >  2:278381/1‑141 (MQ=255)
                                                                     aaaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATAGTTTTTCGATATAAACATTGtt                                                                         >  1:303610/1‑141 (MQ=255)
                                                                      aaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACTTTTTCGATATAAACCTTTTTc                                                                        >  2:33779/1‑141 (MQ=255)
                                                                      aaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTc                                                                        <  2:676199/141‑1 (MQ=255)
                                                                       aTTTCAAATTGCCGTAATTGAACGTATATTTCCTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCt                                                                       <  1:39866/141‑1 (MQ=255)
                                                                        tttCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTTTTCTa                                                                      >  1:702189/1‑141 (MQ=255)
                                                                        tttCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTa                                                                      >  1:304389/1‑141 (MQ=255)
                                                                        tttCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTa                                                                      >  1:502967/1‑141 (MQ=255)
                                                                             aaTTGCCGTAATTGAACGTATATTTCTTCTACAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATATGTTTTCGATATAAAACTTTTTTTAAAtct                                                                 >  1:183472/1‑141 (MQ=255)
                                                                              aTTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAATTTAACATACGTTTTCGATATAAACCTTGTTCTAActctt                                                                >  2:457899/1‑140 (MQ=255)
                                                                                tGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCaa                                                             >  1:443802/1‑141 (MQ=255)
                                                                                    gTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACATTGTTCTAAATCTCaataat                                                         >  2:489994/1‑141 (MQ=255)
                                                                                            aCGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCt                                                <  2:232470/141‑1 (MQ=255)
                                                                                            aCGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCt                                                <  2:746177/141‑1 (MQ=255)
                                                                                             cGTATATTTCTTCTTCAACTATTATTTCATATTTAGTATAGTATCTTTGTTATCACGTTCATCAGTGTGGAAATAATTCGCTATCCTCTGTCATACTTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTg                                               <  2:350155/141‑1 (MQ=255)
                                                                                              gTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGt                                              <  1:470687/141‑1 (MQ=255)
                                                                                              gTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGt                                              <  1:461613/141‑1 (MQ=255)
                                                                                                      cttcttcAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTCTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCg                                      <  1:213309/141‑1 (MQ=255)
                                                                                                      cttcttcAACTATTAGTTCATCTTTAGCATAATCTATATATAAAATTTTCTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCg                                      <  1:213139/141‑1 (MQ=255)
                                                                                                        tcttcAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCACTAATTTTTTGCTGTTTTCATCGTc                                    >  2:760860/1‑141 (MQ=255)
                                                                                                         cttcAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCa                                   <  2:671814/141‑1 (MQ=255)
                                                                                                         cttcAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCa                                   <  2:131742/141‑1 (MQ=255)
                                                                                                          ttcAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAGTAATTCGCTATAACTTAAAATACGTTTTCGATATAAACCTTGTTCTAAATATCAATAATTTTTTGCTGTTTTGATTGTCAt                                  >  1:254246/1‑141 (MQ=255)
                                                                                                            cAACTATTATTTCATCTTTAGCATAATCTATATGTAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGttatta                                <  1:521445/141‑1 (MQ=255)
                                                                                                              aCTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGt                              >  1:714437/1‑141 (MQ=255)
                                                                                                                  ttattTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCTTCGTTATTAGTTaaa                          >  2:115200/1‑141 (MQ=255)
                                                                                                                        catcttcggCATAGTCTAAGTGTAAAATTTTATGCTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAATCACTTACAATTTCTTGATTTTTTTCGTCGTCATCAGTTaaaaagt                      >  2:277158/9‑139 (MQ=255)
                                                                                                                        catcttcgacATAGTCTAAGTATAAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACGAACCTTTTTCTAAATCACTTACAATTTCTTGATTTTTTTCATCATCATCCGTTaaaaaat                      >  1:397679/10‑140 (MQ=255)
                                                                                                                        cATCTTCGGCATCGTCTAAGTTTAAAATTTTATGCTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAATCACTTACAATTTCTTGATTTTTTTCGTCGTCATCAGTTaaaaagt                      <  1:468003/141‑3 (MQ=255)
                                                                                                                         atcgtcggCATAGTCTAAGTTTAAAATTTTATGCTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAATCACTTACAATTTTTTGATTTTTTTCGTCGTCATCAGTTaaaaagtc                     <  1:219701/134‑4 (MQ=255)
                                                                                                                              tAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTTTAAATCTCAATAATTTTTTGCTGTTTTTATCGTCATTAGTTAAAAAAATAATTTat              >  2:564485/1‑140 (MQ=255)
                                                                                                                                gCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACt            >  2:546254/1‑141 (MQ=255)
                                                                                                                                gCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTAGATATAAACCTTGTTATAAATCTCACTAATTTTTTGCTGTTTTCATCGGCATTAGTTAAAAAAATAATTTCACt            >  2:149684/1‑141 (MQ=255)
                                                                                                                                  aTAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATACGTTCACATGCGTTTTTGATATAAACCGTGTTCTAAATCTCAAAAATTTTTTTCTGTTTTCATCGTCCTTAGTTATAAAAATAATTTAATTgg          >  1:368957/1‑140 (MQ=255)
                                                                                                                                  aTAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGa          <  1:188527/141‑1 (MQ=255)
                                                                                                                                     aTCTATATATAAAATTTTATGCTATTATTTAAATAATTCACTATAATTTAAACTCCATTTTCGATATAAACCTTGTTCTAAATCTCAATAATCTTTTGCTGTTTTTATCTTCATTAGTTAAAAAAATCAtttcaatgtctt       >  2:446534/1‑132 (MQ=255)
                                                                                                                                      tCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTTAATAACTATTTTCTGGTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGttt      <  2:304389/141‑1 (MQ=255)
                                                                                                                                         ttatataAATTTTTATGCTATTATTTAATTAATTCGCTATAACTTAACATACTTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTg   <  2:59289/140‑1 (MQ=255)
                                                                                                                                         atatataAAATTTTATGGTGTTATTTGAGTAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTg   >  1:185125/1‑141 (MQ=255)
                                                                                                                                         atatataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATCCGTTTTCGATATAAACCTTGTTCTAAATCTCCATCAACTTTTGCTGCGTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTg   <  2:735278/141‑1 (MQ=255)
                                                                                                                                         atatataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTTGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGttttt  >  2:141951/1‑140 (MQ=255)
                                                                                                                                          tatataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTga  <  1:624711/141‑1 (MQ=255)
                                                                                                                                          tatataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTga  <  2:687252/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                               
GAACGCTTCAATTTCCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGA  >  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: