Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,680,121 C→T 13.9% intergenic (+173/+359) cbpD → / ← 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,1210CT13.9% 184.1 / 18.7 72intergenic (+173/+359)cbpD/USA300HOU_2537possible secretory antigen/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/35);  new base T (5/5);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

CCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATT  >  CP000730/2679981‑2680259
                                                                                                                                            |                                                                                                                                              
ccTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGc                                                                                                                                                >  2:165389/1‑141 (MQ=255)
 cTGTATGCGTGTGATAGAGCTGGTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCt                                                                                                                                               <  1:359546/141‑1 (MQ=255)
 cTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCt                                                                                                                                               <  1:255023/141‑1 (MQ=255)
 cTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCt                                                                                                                                               <  2:447868/141‑1 (MQ=255)
 cTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCt                                                                                                                                               >  1:118726/1‑141 (MQ=255)
 cTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCt                                                                                                                                               >  1:45477/1‑141 (MQ=255)
 cTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCt                                                                                                                                               >  1:191955/1‑141 (MQ=255)
  tGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCta                                                                                                                                              <  1:680571/141‑1 (MQ=255)
     aTGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCtatta                                                                                                                                           >  1:562332/1‑141 (MQ=255)
        cGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTAttt                                                                                                                                        <  2:422737/141‑1 (MQ=255)
           gtgATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTaaa                                                                                                                                     >  2:154507/1‑141 (MQ=255)
             gATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAaata                                                                                                                                   >  2:89440/1‑141 (MQ=255)
               tAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTTAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAaataat                                                                                                                                 <  2:101131/141‑1 (MQ=255)
               tAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAaataat                                                                                                                                 <  1:151485/141‑1 (MQ=255)
                  aaGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCg                                                                                                                              >  2:71232/1‑141 (MQ=255)
                    gTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCt                                                                                                                            <  2:335580/141‑1 (MQ=255)
                       tttattttatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCtata                                                                                                                         >  1:268760/1‑141 (MQ=255)
                           ttttatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACtt                                                                                                                     >  1:90242/1‑141 (MQ=255)
                           ttttatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACtt                                                                                                                     >  1:396059/1‑141 (MQ=255)
                              tatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAAc                                                                                                                  <  1:232305/141‑1 (MQ=255)
                              tatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAAc                                                                                                                  <  1:19903/141‑1 (MQ=255)
                              tatGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAAc                                                                                                                  >  1:194153/1‑141 (MQ=255)
                               atGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACa                                                                                                                 <  1:304368/141‑1 (MQ=255)
                                   aaTTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACg                                                                                                             >  2:114913/1‑141 (MQ=255)
                                           tattaCTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGat                                                                                                     >  1:81191/1‑141 (MQ=255)
                                              taCTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGAtata                                                                                                  <  1:681520/141‑1 (MQ=255)
                                               tcTTCTACAAATTTCAAATTGCCCTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATaa                                                                                                 <  2:643580/140‑1 (MQ=255)
                                                tttCTACAATTTTCAGATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATACTATTATTTAAATAATTCGCTATAACTTAACATGCGTTTTCGATATaaa                                                                                                <  2:317542/140‑1 (MQ=255)
                                                   cTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCt                                                                                             <  1:165389/141‑1 (MQ=255)
                                                   cTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCt                                                                                             >  2:386494/1‑141 (MQ=255)
                                                     aCAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTg                                                                                           <  2:665491/141‑1 (MQ=255)
                                                       aaaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGtt                                                                                         <  1:604382/141‑1 (MQ=255)
                                                        aaTTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAAAAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTc                                                                                        >  1:180995/1‑141 (MQ=255)
                                                              aaaTTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAAtc                                                                                  <  1:679995/141‑1 (MQ=255)
                                                                           tGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAAtttttt                                                                  <  2:194153/141‑1 (MQ=255)
                                                                                 tataTTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGtt                                                           <  2:191955/141‑1 (MQ=255)
                                                                                     tttCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCa                                                       <  2:51644/141‑1 (MQ=255)
                                                                                               aaCTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTAGATTTATACCTTGTTCTATCTCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAg                                             >  1:491093/1‑141 (MQ=255)
                                                                                                 cTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGtt                                           <  1:71232/141‑1 (MQ=255)
                                                                                                  caaattaccATCTTCGGCATCGTCTAAGTTTAAAATTTTATGCTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACTAACCTTTTTCTAAATCACT‑‑TACAATTTCTTGATTTTTTTCGTCGTCATCAGt                                            <  2:185798/133‑1 (MQ=255)
                                                                                                    ttattTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTaaa                                        >  1:243904/1‑141 (MQ=255)
                                                                                                    ttattTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTaaa                                        >  2:425276/1‑141 (MQ=255)
                                                                                                       tttCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTaaaaaa                                     <  1:662660/141‑1 (MQ=255)
                                                                                                       tttCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTaaaaaa                                     <  1:662637/141‑1 (MQ=255)
                                                                                                       tttCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTaaaaaa                                     <  1:662634/141‑1 (MQ=255)
                                                                                                           aTCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAaataa                                 >  1:474850/1‑141 (MQ=255)
                                                                                                                tAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTaa                            >  1:662603/1‑141 (MQ=255)
                                                                                                                 aGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAAc                           >  2:427991/1‑141 (MQ=255)
                                                                                                                 aGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAAc                           >  2:295574/1‑141 (MQ=255)
                                                                                                                  gCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACt                          <  1:401871/141‑1 (MQ=255)
                                                                                                                     tAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAATTGAg                       >  1:394452/1‑141 (MQ=255)
                                                                                                                     tAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAg                       >  1:525950/1‑141 (MQ=255)
                                                                                                                     tAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAg                       <  1:295574/141‑1 (MQ=255)
                                                                                                                     tAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAg                       <  1:583869/141‑1 (MQ=255)
                                                                                                                       aTCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGtt                     >  1:87450/1‑141 (MQ=255)
                                                                                                                         ctaagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAAcc                                                                                              >  2:6147/8‑70 (MQ=255)
                                                                                                                         ctaagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAAcc                                                                                              <  1:6147/63‑1 (MQ=255)
                                                                                                                          taagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGata                                                                                                    >  1:401215/7‑63 (MQ=255)
                                                                                                                          taagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGata                                                                                                    <  2:401215/57‑1 (MQ=255)
                                                                                                                           aagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGatac                                                                                                   >  1:408277/6‑62 (MQ=255)
                                                                                                                           aagtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGatac                                                                                                   <  2:408277/58‑2 (MQ=255)
                                                                                                                            agtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAAtcactt                                                                              >  1:523125/5‑79 (MQ=255)
                                                                                                                            agtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAAtcactt                                                                              <  2:523125/79‑5 (MQ=255)
                                                                                                                             gtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAAtcactta                                                                             <  1:669334/80‑6 (MQ=255)
                                                                                                                             gtttaAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAAtcactta                                                                             >  2:669334/4‑78 (MQ=255)
                                                                                                                               ataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGAc             <  1:126136/141‑1 (MQ=255)
                                                                                                                                  aaaTTTTATGCTATTATTTAAATAATTAGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTT‑‑‑AATAATTTTTTTCTGTTTTCATAGTAATAAGATAAAAAAATAATTTAACTGAGTTTTTGAGACCTg          <  2:383740/141‑1 (MQ=255)
                                                                                                                                   aaTTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTCAATAATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGAcc            <  1:154507/141‑1 (MQ=255)
                                                                                                                                      tttATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTa      <  2:708587/141‑1 (MQ=255)
                                                                                                                                       ttATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGAAAAAAAAAAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTaa     <  2:512610/141‑1 (MQ=255)
                                                                                                                                         aTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAAt   >  2:381549/1‑141 (MQ=255)
                                                                                                                                         aTGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAANTTAACTGAGTTTNTGAGACCTGTTTAAAt   <  1:394475/141‑1 (MQ=255)
                                                                                                                                          tGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAAtt  >  1:610512/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CCTGTATGCGTGTGATAGAAGTTTTTATTTTATGAAATTATATTATTACTTCTACAAATTTCAAATTGCCGTAATTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑‑‑AATAATTTTTTGCTGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATT  >  CP000730/2679981‑2680259

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