Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,581 C→T 42.2% intergenic (‑30/‑493) USA300HOU_0026 ← / → USA300HOU_0028 IS431mec transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,5810CT42.2% 58.4 / 81.8 79intergenic (‑30/‑493)USA300HOU_0026/USA300HOU_0028IS431mec transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/24);  new base T (12/21);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

AATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTATTTAG  >  CP000730/36448‑36720
                                                                                                                                     |                                                                                                                                           
aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAAt                                                                                                                                      <  2:120899/141‑1 (MQ=21)
   atTTCAGTTATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                                                                                                                                   <  2:341181/141‑1 (MQ=25)
   aattCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                                                                                                                                   <  2:313371/139‑1 (MQ=25)
             tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTACAAc                                                                                                                         >  2:390579/1‑141 (MQ=25)
             aatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                                                                                                                         <  1:334007/140‑1 (MQ=25)
               ttacaatatCTGAATGCATATCTTAGATAATAGCCAACTGCTACAGTTATAACATCCTTGTTAAATTTTTTATATCTGAAATAGTTCATACATAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAg                                                                                                                       <  1:132207/136‑1 (MQ=1)
                     tagcccaaggcacatCTTAGATAGTAGCCAACGGCTACAGTTATAACATCCTTGTTAAATTTTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                                                                                                 <  2:255695/128‑1 (MQ=11)
                     tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACATTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                                                                                                 >  2:231374/4‑141 (MQ=9)
                     tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                                                                                                 <  2:266808/138‑1 (MQ=9)
                     tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                                                                                                 <  1:114297/138‑1 (MQ=9)
                     tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                                                                                                 <  1:341964/138‑1 (MQ=12)
                     tAACTCAATGTATATCTTAGATAGTAGCCAACGTCTACAGTCATAACATACTATTTAAATTGTTTACATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt                                                                                                                 <  2:295754/141‑1 (MQ=11)
                     tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt                                                                                                                 <  2:312671/141‑1 (MQ=37)
                     tAACTCAATGTATATCTTAGATAGTAGACAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAGACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt                                                                                                                 <  2:391034/141‑1 (MQ=255)
                     atactcgaggcatatCTTAGATAGTAGACAGCAGCTACAGTGCTAACATACTTGTTAAATTGTTTATATCTGAAATAGTTCACACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCACCAGAAcctc                                                                                                                 <  2:335425/130‑3 (MQ=11)
                     aaaCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt                                                                                                                 <  2:245004/140‑1 (MQ=37)
                        gtCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                                                                                                              >  2:102875/2‑141 (MQ=16)
                        cgcAATGCGTATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                                                                                                              <  2:187865/139‑1 (MQ=255)
                        cTCAATGCATATCTTAGATAGTCGCCAACGGCTACAGTGATAACATCCTTGTTAAGTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACAGtatt                                                                                                              >  2:288262/1‑141 (MQ=11)
                        cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGCAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                                                                                                              >  1:150085/1‑141 (MQ=12)
                        cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAACTTATAGTATAAATTCAACTTTGCAACATAACCGtgat                                                                                                              >  1:138063/1‑138 (MQ=11)
                         tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                                                                                                             >  1:26978/1‑140 (MQ=12)
                          cAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtattat                                                                                                            >  1:25919/1‑141 (MQ=40)
                          cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt                                                                                                            >  1:32293/1‑139 (MQ=12)
                          cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTCAAAATTAGATTATAAATTCAACTTTGCAACAGAACCGtatttt                                                                                                            <  1:21100/141‑3 (MQ=11)
                          cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTAAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTTCAACAGAACCGtatttt                                                                                                            >  1:261939/1‑139 (MQ=11)
                           aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATTCTTGTTAAATTGTTTATAGCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta                                                                                                           <  1:380664/141‑4 (MQ=11)
                           aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatttta                                                                                                           >  2:321788/1‑138 (MQ=21)
                           aaTGCATATCGTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTACATTGTTTGTAACTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGATCCGtatttta                                                                                                           <  1:91112/141‑4 (MQ=11)
                            atgcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat                                                                                                          >  1:151237/5‑137 (MQ=255)
                             tgcatatCTTAGATAGTAGCCAACTGCTACAGTGATAACATCCTTGTTAATTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatg                                                                                                         <  2:258476/138‑6 (MQ=11)
                              tcatatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTTAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactat                                                                                                        <  2:227239/139‑12 (MQ=11)
                              gcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatga                                                                                                        <  2:103847/139‑7 (MQ=255)
                              gcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatga                                                                                                        <  1:84353/139‑7 (MQ=255)
                               catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtat                                                                                                       <  1:132482/140‑13 (MQ=11)
                               catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgag                                                                                                       >  1:100155/2‑134 (MQ=255)
                               catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgag                                                                                                       <  1:38526/140‑8 (MQ=255)
                               catatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcctttcatata                                                                                                       >  1:346257/2‑129 (MQ=11)
                               catatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactata                                                                                                       <  1:298743/140‑13 (MQ=11)
                               catatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGATAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactata                                                                                                       >  2:216056/2‑129 (MQ=11)
                                atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagt                                                                                                      >  1:181274/1‑133 (MQ=255)
                                atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAGTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatt                                                                                                      >  1:329376/1‑128 (MQ=11)
                                 tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatta                                                                                                     >  1:120056/1‑127 (MQ=11)
                                 tatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                                                                                                     >  1:335379/1‑127 (MQ=11)
                                 tatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                                                                                                     <  2:198132/141‑15 (MQ=11)
                                  attttAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATaga                                                                                                    <  2:335219/138‑1 (MQ=255)
                                  atCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAATTAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACGTTGCAACGGAACCGtattttatgagtac                                                                                                    <  1:256203/141‑11 (MQ=11)
                                   tCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtaca                                                                                                   <  2:172975/141‑12 (MQ=255)
                                   tCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtaca                                                                                                   <  2:29210/141‑12 (MQ=11)
                                   tCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGCTTATATCGGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCTtattttatgagtaca                                                                                                   <  1:291933/141‑12 (MQ=11)
                                   tCTTAGATAGTAGCCAACCGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtaca                                                                                                   <  1:284944/141‑12 (MQ=11)
                                     ttAGATAGTAGCCAGCGGCTACAGTGATAACATCTTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGAt                                                                                                 <  2:224247/141‑1 (MQ=255)
                                     ttAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact                                                                                                 <  2:157853/141‑14 (MQ=11)
                                            gtagCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTa                                                                                          >  1:361194/1‑141 (MQ=255)
                                              agCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATGCCACGTTGCAACAGAACCGAATTATGGAATAGACATGTTGGTAAc                                                                                        >  1:246424/1‑141 (MQ=17)
                                               gCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACa                                                                                       <  1:60889/141‑1 (MQ=255)
                                                   accgCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACCTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTta                                                                                   <  2:156645/138‑1 (MQ=255)
                                                           aGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGAt                                                                           >  2:292694/1‑141 (MQ=255)
                                                                 aaCATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATtat                                                                     >  1:79161/1‑141 (MQ=255)
                                                                      ccTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTATAATTTTACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTa                                                                >  2:181657/1‑141 (MQ=255)
                                                                      ccTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTa                                                                >  2:354561/1‑141 (MQ=255)
                                                                        ttGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAg                                                              <  2:393003/141‑1 (MQ=255)
                                                                           ttAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGttt                                                           <  2:251940/141‑1 (MQ=255)
                                                                              aaTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTTTTTTTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAGCATTTATACAGGATCATTATACTTAAGTTTAAt                                                        >  2:309765/1‑141 (MQ=25)
                                                                                                  gTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACcaca                                    <  1:149589/141‑1 (MQ=255)
                                                                                                  gTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGGTTGATTTAGTTATTTCGGAACcaga                                    >  2:366364/1‑139 (MQ=37)
                                                                                                   ttCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACcacac                                   >  1:272460/1‑141 (MQ=255)
                                                                                                   ttCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACcacac                                   >  1:134226/1‑141 (MQ=255)
                                                                                                   ttCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTCATTACAGAACcacac                                   >  2:77776/1‑141 (MQ=255)
                                                                                                            aagaCTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAAcc                          >  1:384148/1‑141 (MQ=255)
                                                                                                            aagaCTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATCACAGAACCACACATTCCAAcc                          <  1:204671/141‑1 (MQ=255)
                                                                                                            aagaCCCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTGTGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTAAAGAACAAAACATTCCAAcc                          >  2:134166/1‑141 (MQ=38)
                                                                                                              gccTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCag                        <  1:224855/139‑1 (MQ=255)
                                                                                                                       ttGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAg               >  2:393610/1‑141 (MQ=255)
                                                                                                                        tGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGCTGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGt              >  2:304334/1‑141 (MQ=255)
                                                                                                                                acataCTAAAAATTCAACTTTGCAACAGAACCGTATTACGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTAt      <  1:140241/139‑1 (MQ=255)
                                                                                                                                   taCTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTATTTa   <  1:354561/141‑1 (MQ=255)
                                                                                                                                    aCTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTATTTAg  <  1:314020/141‑1 (MQ=255)
                                                                                                                                    aCTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTATTTAg  <  1:181657/141‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                           
AATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTATTTAG  >  CP000730/36448‑36720

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