Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,680,132 A→G 8.1% intergenic (+184/+348) 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,1320AG 8.1% 188.2 / 3.7 70intergenic (+184/+348)USA300HOU_2536/USA300HOU_2537possible secretory antigen/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/31);  new base G (4/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCA  >  CP000730/2679997‑2680268
                                                                                                                                        |                                                                                                                                          
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:502967/1‑141 (MQ=255)
                                          tttCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTa                                                                                             >  1:304389/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)
                                                                                                           atatataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTTGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGttttta                          >  2:141951/1‑140 (MQ=255)
                                                                                                            tatataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTga                         <  2:687252/141‑1 (MQ=255)
                                                                                                            tatataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTga                         <  1:624711/141‑1 (MQ=255)
                                                                                                             atataAAATTTTTGGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTgag                        <  2:291856/141‑1 (MQ=255)
                                                                                                                ataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGAc                      <  1:647447/141‑1 (MQ=255)
                                                                                                                 taAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGTGAcc                     >  1:294142/1‑141 (MQ=255)
                                                                                                                       tttATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTa               <  1:251338/141‑1 (MQ=255)
                                                                                                                          aTGCTATTATTTAAATAATTCGCTATAACTTAACCTAGGCCGATGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAAt            >  1:181540/1‑141 (MQ=255)
                                                                                                                          aTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAATTGATTTTTTgtttactgttttaat            >  2:306007/1‑127 (MQ=255)
                                                                                                                          aTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAAt            >  2:251906/1‑141 (MQ=255)
                                                                                                                          aTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAAt            >  2:131039/1‑141 (MQ=255)
                                                                                                                              tattatTTAAATAATTAGCTACAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTGAATtgat        >  1:395127/1‑141 (MQ=255)
                                                                                                                              tattagTTAAATAATTGGCTCTAACTTAACAAGCGTTTGCGATAGAAGCCTGATTCTAAATCTAAATAATTGTTGGTTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTTAGTTTTTGAGACCTGTTTAAATtgat        <  2:169744/141‑1 (MQ=255)
                                                                                                                                ttattTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAACAtgatga      >  1:109170/1‑141 (MQ=255)
                                                                                                                                  attTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACt    <  1:99653/141‑1 (MQ=255)
                                                                                                                                   tttAAACAATTCGTTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTc   <  2:147816/141‑1 (MQ=255)
                                                                                                                                    ttAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAAACTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCTGTATTAATTAAAAAAATAATTTAACTGTGTTTTTAAGACCGTTTTAAATTCATACCTCa  >  2:153173/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                          
AGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCA  >  CP000730/2679997‑2680268

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: