Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,303 A→T 10.2% E717D (GAA→GAT sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273557,3030AT10.2% 856.5 / 64.3 262E717D (GAA→GATsfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base A (122/111);  new base T (13/14);  total (135/127)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

GGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCA  >  CP009273/557154‑557452
                                                                                                                                                     |                                                                                                                                                     
gggTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCGGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaa                                                                                                                                                       <  1:2777751/150‑1 (MQ=255)
 ggTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaa                                                                                                                                                      <  2:4085982/150‑1 (MQ=255)
    ataACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATc                                                                                                                                                   >  1:3285641/1‑150 (MQ=255)
       aCTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCaaa                                                                                                                                                >  1:3727205/1‑150 (MQ=255)
           cGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTgg                                                                                                                                            <  2:3471236/150‑1 (MQ=255)
            gACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGc                                                                                                                                           >  2:1816057/1‑150 (MQ=255)
               cgcgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAtt                                                                                                                                        <  1:1196076/150‑1 (MQ=255)
                gcgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCAttt                                                                                                                                       <  1:3808027/150‑1 (MQ=255)
                gcgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAttt                                                                                                                                       <  1:3439527/150‑1 (MQ=255)
                 cgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAtttt                                                                                                                                      <  2:736599/150‑1 (MQ=255)
                 cgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAtttt                                                                                                                                      <  2:372228/150‑1 (MQ=255)
                 cgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAtttt                                                                                                                                      <  1:967715/150‑1 (MQ=255)
                  gATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTa                                                                                                                                     >  1:4080123/1‑150 (MQ=255)
                   aTCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTaa                                                                                                                                    >  2:4060170/1‑150 (MQ=255)
                   aTCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTaa                                                                                                                                    >  1:3482553/1‑150 (MQ=255)
                     cAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAAcc                                                                                                                                  >  1:178233/1‑150 (MQ=255)
                       aCATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGa                                                                                                                                >  1:804696/1‑150 (MQ=255)
                         aTGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACt                                                                                                                              <  1:135346/150‑1 (MQ=255)
                          tGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTg                                                                                                                             <  2:2946194/150‑1 (MQ=255)
                          tGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTg                                                                                                                             <  2:3735703/150‑1 (MQ=255)
                             cGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGgcg                                                                                                                          >  2:1216835/1‑150 (MQ=255)
                               ttAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGcgcg                                                                                                                        <  1:2470630/150‑1 (MQ=255)
                                tAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCgg                                                                                                                       >  1:2321587/1‑150 (MQ=255)
                                   cTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTa                                                                                                                    <  2:466737/150‑1 (MQ=255)
                                   cTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTa                                                                                                                    <  2:929972/150‑1 (MQ=255)
                                   cTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTa                                                                                                                    <  1:2883064/150‑1 (MQ=255)
                                    tGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTAt                                                                                                                   >  1:2596661/1‑150 (MQ=255)
                                     ggCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATg                                                                                                                  >  1:742308/1‑150 (MQ=255)
                                     ggCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATg                                                                                                                  >  1:1176273/1‑150 (MQ=255)
                                       cAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATACCGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCg                                                                                                                <  1:1074168/150‑1 (MQ=255)
                                         gCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGt                                                                                                              <  1:1817393/150‑1 (MQ=255)
                                            gTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGAt                                                                                                           <  2:2982498/150‑1 (MQ=255)
                                            gTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGAt                                                                                                           <  2:3120887/150‑1 (MQ=255)
                                            gTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTAAGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGAt                                                                                                           <  2:932834/150‑1 (MQ=255)
                                             tCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATg                                                                                                          <  2:3303091/150‑1 (MQ=255)
                                             tCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATg                                                                                                          <  2:4112705/150‑1 (MQ=255)
                                               tGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCt                                                                                                        >  1:1045011/1‑150 (MQ=255)
                                                    gtgtGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATg                                                                                                   >  2:2789663/1‑150 (MQ=255)
                                                    gtgtGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATg                                                                                                   >  2:616641/1‑150 (MQ=255)
                                                      gtgGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGcc                                                                                                 >  1:1589186/1‑150 (MQ=255)
                                                       tgGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCa                                                                                                >  2:3902533/1‑150 (MQ=255)
                                                        ggTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCAc                                                                                               >  2:2462818/1‑150 (MQ=255)
                                                             ggCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGttt                                                                                          <  2:3286510/150‑1 (MQ=255)
                                                              gCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTa                                                                                         >  2:1007520/1‑150 (MQ=255)
                                                                aTGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATc                                                                                       >  2:3225272/1‑150 (MQ=255)
                                                                aTGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTATACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATc                                                                                       >  2:2842112/1‑150 (MQ=255)
                                                                 tGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCg                                                                                      <  2:65678/150‑1 (MQ=255)
                                                                  gAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCgg                                                                                     <  1:2417758/150‑1 (MQ=255)
                                                                  gAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCgg                                                                                     >  1:1690535/1‑150 (MQ=255)
                                                                   aaaaTGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGt                                                                                    <  2:3832502/150‑1 (MQ=255)
                                                                    aaaTGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGta                                                                                   <  1:2852876/150‑1 (MQ=255)
                                                                    aaaTGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGta                                                                                   >  1:4028042/1‑150 (MQ=255)
                                                                    aaaTGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGta                                                                                   <  2:81725/150‑1 (MQ=255)
                                                                          cATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAAACAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCg                                                                             >  2:185575/1‑150 (MQ=255)
                                                                          cATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCACCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCg                                                                             <  1:414663/150‑1 (MQ=255)
                                                                            tAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTCTAACCGTa                                                                           >  2:4333218/1‑150 (MQ=255)
                                                                            tAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTa                                                                           >  1:3659978/1‑150 (MQ=255)
                                                                             aaCGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTAt                                                                          >  1:1292742/1‑150 (MQ=255)
                                                                             aaCGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTAt                                                                          >  1:352253/1‑150 (MQ=255)
                                                                                 cTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATcgcg                                                                      >  2:2036533/1‑150 (MQ=255)
                                                                                 cTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATcgcg                                                                      >  1:1329900/1‑150 (MQ=255)
                                                                                    agccagccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCtt                                                                   <  2:712411/150‑1 (MQ=255)
                                                                                    agccagccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCtt                                                                   >  1:1935836/1‑150 (MQ=255)
                                                                                          ccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATAcc                                                             >  2:1912471/1‑150 (MQ=255)
                                                                                               aGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATAc                                                        <  1:4265395/150‑1 (MQ=255)
                                                                                               aGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATAc                                                        >  2:1471782/1‑150 (MQ=255)
                                                                                                  ggATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGag                                                     <  1:4327457/150‑2 (MQ=255)
                                                                                                  ggATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGAt                                                     >  1:918809/1‑150 (MQ=255)
                                                                                                        cAATGTTTGGATGAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGaa                                               <  2:4299979/150‑1 (MQ=255)
                                                                                                             ttttGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCa                                          <  1:3845227/150‑1 (MQ=255)
                                                                                                               ttGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATc                                        <  2:1589186/150‑1 (MQ=255)
                                                                                                               ttGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATc                                        <  1:2122180/150‑1 (MQ=255)
                                                                                                                tGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCg                                       <  1:1102822/150‑1 (MQ=255)
                                                                                                                 gATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGa                                      >  2:368347/1‑150 (MQ=255)
                                                                                                                 gATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGa                                      >  1:453611/1‑150 (MQ=255)
                                                                                                                  aTTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGACTTCCATCGAt                                     >  1:1914211/1‑150 (MQ=255)
                                                                                                                    tAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGt                                   >  1:3206628/1‑150 (MQ=255)
                                                                                                                     aaaGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGtt                                  >  1:1863493/1‑150 (MQ=255)
                                                                                                                     aaaGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGtt                                  <  1:2240612/150‑1 (MQ=255)
                                                                                                                       aGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGa                                >  2:4196147/1‑150 (MQ=255)
                                                                                                                       aGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGa                                >  2:1998840/1‑150 (MQ=255)
                                                                                                                         cgcCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaa                              >  2:221991/1‑150 (MQ=255)
                                                                                                                             tGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaaaaaa                          >  2:963735/1‑150 (MQ=255)
                                                                                                                             tGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaaaaaa                          >  2:1674092/1‑150 (MQ=255)
                                                                                                                              ggCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAat                         <  2:3595772/150‑1 (MQ=255)
                                                                                                                              ggCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAat                         <  2:1401675/150‑1 (MQ=255)
                                                                                                                               gcgcAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAata                        <  1:1455910/150‑1 (MQ=255)
                                                                                                                                  cAGGCGGTGTACGCATTGACCACCACACTGGCAAGTTACCCGCCTGGCGCGGCTAGGCGGTGATGCTGTATGCCACGGTTTATCGGTATCACCGCATCGCGCTTGATACCAATCCCCGGGGGAATTCCATCGATGATGAAAAAAATATTa                     <  2:2642095/150‑1 (MQ=255)
                                                                                                                                      cGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTagcag                 <  2:2212421/150‑1 (MQ=255)
                                                                                                                                       ggTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGACTTCCATCGATGTTGAAAAAAATATTagcagc                >  2:4011405/1‑150 (MQ=255)
                                                                                                                                        gtgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCg               >  2:2464591/1‑150 (MQ=255)
                                                                                                                                         tgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGt              <  2:2371746/150‑1 (MQ=255)
                                                                                                                                         tgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGt              <  2:1047824/150‑1 (MQ=255)
                                                                                                                                          gtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGtt             <  2:1000026/150‑1 (MQ=255)
                                                                                                                                            aCGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTtgt           <  1:3147256/150‑1 (MQ=255)
                                                                                                                                            aCGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTtgt           <  1:2953315/150‑1 (MQ=255)
                                                                                                                                                  tGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGAc     >  1:1789726/1‑150 (MQ=255)
                                                                                                                                                    aaaaTCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGc   >  2:1239983/1‑150 (MQ=255)
                                                                                                                                                     aaaTCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCa  >  1:2108255/1‑150 (MQ=255)
                                                                                                                                                     aaaTCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCa  >  1:3288239/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCA  >  CP009273/557154‑557452

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