Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,303 A→T 30.8% 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,3030AT30.8% 479.3 / 306.6 305E717D (GAA→GAT) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base A (86/125);  new base T (46/48);  total (132/173)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 310 total aligned reads displayed.

GGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGAC  >  CP009273/557154‑557449
                                                                                                                                                     |                                                                                                                                                  
  gTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaaa                                                                                                                                                  <  1:1347378/150‑1 (MQ=255)
   tataACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAAt                                                                                                                                                 <  1:1882706/150‑1 (MQ=255)
   tataACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAAt                                                                                                                                                 >  2:394065/1‑150 (MQ=255)
    ataACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATc                                                                                                                                                >  1:1023148/1‑150 (MQ=255)
         tACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACt                                                                                                                                           <  1:3099827/150‑1 (MQ=255)
             aCCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCa                                                                                                                                       <  2:4089875/150‑1 (MQ=255)
                 cgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAtttt                                                                                                                                   >  2:3146706/1‑150 (MQ=255)
                 cgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAtttt                                                                                                                                   <  1:2371962/150‑1 (MQ=255)
                    tCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAAc                                                                                                                                >  2:1416470/1‑150 (MQ=255)
                     cAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAAcc                                                                                                                               >  1:3671113/1‑150 (MQ=255)
                     cAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAAcc                                                                                                                               <  2:4437401/150‑1 (MQ=255)
                     cAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAAcc                                                                                                                               >  1:2624258/1‑150 (MQ=255)
                        cATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGAc                                                                                                                            >  2:2010520/1‑150 (MQ=255)
                          tGACGTTATCTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTg                                                                                                                          <  2:1389762/150‑1 (MQ=255)
                            aCGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGgc                                                                                                                        >  1:4045741/1‑150 (MQ=255)
                            aCGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGgc                                                                                                                        <  1:2792987/150‑1 (MQ=255)
                            aCGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGgc                                                                                                                        <  2:2168920/150‑1 (MQ=255)
                             cGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGgcg                                                                                                                       <  2:3688687/150‑1 (MQ=255)
                                   cTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTa                                                                                                                 >  2:575017/1‑150 (MQ=255)
                                    tGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTAt                                                                                                                <  2:3462796/150‑1 (MQ=255)
                                    tGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTAt                                                                                                                <  2:536263/150‑1 (MQ=255)
                                     ggCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATg                                                                                                               >  2:1456935/1‑150 (MQ=255)
                                         gCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGt                                                                                                           <  1:259665/150‑1 (MQ=255)
                                           tGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGa                                                                                                         >  2:1267703/1‑150 (MQ=255)
                                            gTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGAt                                                                                                        <  2:1859422/150‑1 (MQ=255)
                                             tCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATg                                                                                                       <  2:550346/150‑1 (MQ=255)
                                             tCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATg                                                                                                       <  2:492839/150‑1 (MQ=255)
                                              cTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGc                                                                                                      <  2:3360011/150‑1 (MQ=255)
                                                ggCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTg                                                                                                    >  2:2642635/1‑150 (MQ=255)
                                                   ggTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTAt                                                                                                 <  2:1563854/150‑1 (MQ=255)
                                                   ggTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTAt                                                                                                 >  2:220543/1‑150 (MQ=255)
                                                   ggTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTAt                                                                                                 <  2:4279631/150‑1 (MQ=255)
                                                    gtgtGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATg                                                                                                <  2:4606833/150‑1 (MQ=255)
                                                     tgtgGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGc                                                                                               <  1:4268834/150‑1 (MQ=255)
                                                        ggTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAGCCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCAc                                                                                            <  2:3172072/150‑1 (MQ=255)
                                                            gggCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGtt                                                                                        >  2:3813321/1‑150 (MQ=255)
                                                             ggCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGttt                                                                                       <  2:1629829/150‑1 (MQ=255)
                                                              gCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTa                                                                                      <  1:3348658/150‑1 (MQ=255)
                                                               cATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTAt                                                                                     <  2:1445729/150‑1 (MQ=255)
                                                                aTGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATc                                                                                    >  2:3308590/1‑150 (MQ=255)
                                                                 tGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCg                                                                                   >  2:2496756/1‑150 (MQ=255)
                                                                  gAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACTCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCgg                                                                                  <  2:4732555/150‑1 (MQ=255)
                                                                        ggCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAAc                                                                            <  1:3914711/150‑1 (MQ=255)
                                                                           aTAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGt                                                                         >  2:63221/1‑150 (MQ=255)
                                                                             aaCGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTAt                                                                       <  1:1142023/150‑1 (MQ=255)
                                                                             aaCGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTAt                                                                       >  1:348492/1‑150 (MQ=255)
                                                                             aaCGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTAt                                                                       >  1:1502565/1‑150 (MQ=255)
                                                                              aCGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATc                                                                      <  2:1780436/150‑1 (MQ=255)
                                                                               cGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATcg                                                                     <  1:2492431/150‑1 (MQ=255)
                                                                                 cTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATcgcg                                                                   >  1:2899116/1‑150 (MQ=255)
                                                                                  tGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCgcgc                                                                  >  1:3531852/1‑150 (MQ=255)
                                                                                   tagccagccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGAATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCt                                                                 >  2:2487288/2‑150 (MQ=255)
                                                                                   gAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCt                                                                 <  1:3895988/150‑1 (MQ=255)
                                                                                    agccagccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCtt                                                                <  1:4134392/150‑1 (MQ=255)
                                                                                     gccagccTGTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTg                                                               >  2:3548324/1‑150 (MQ=255)
                                                                                      ccagccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGa                                                              >  1:88414/1‑150 (MQ=255)
                                                                                      ccagccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGa                                                              <  1:1693952/150‑1 (MQ=255)
                                                                                       cagccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGAt                                                             >  2:3028764/1‑150 (MQ=255)
                                                                                       cagccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGAt                                                             <  2:2340160/150‑1 (MQ=255)
                                                                                        agccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATa                                                            >  2:1286199/1‑150 (MQ=255)
                                                                                          ccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCAGTATAACCGTATCGCGCTTGATAcc                                                          >  2:3248311/1‑150 (MQ=255)
                                                                                          ccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATAcc                                                          >  2:4466579/1‑150 (MQ=255)
                                                                                              tAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATa                                                      <  1:3090670/150‑1 (MQ=255)
                                                                                              tAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATa                                                      <  2:1023833/150‑1 (MQ=255)
                                                                                               atgggATAACAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACTAATTGATAATCAATCTGGCATTTTAACAGTCTGGCGCGGCTATGCGGTGATTCTGTTTGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATAc                                                     <  2:3554181/148‑1 (MQ=255)
                                                                                                ggggATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACg                                                    >  2:768041/1‑150 (MQ=255)
                                                                                                   gATACCAATGTTTTGATTAAAGCGCCTGGCGCTGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATg                                                 <  1:470028/150‑1 (MQ=255)
                                                                                                      aCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATgggg                                              <  2:4379381/150‑1 (MQ=255)
                                                                                                       ccAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGa                                             <  1:312073/150‑1 (MQ=255)
                                                                                                        cAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGaa                                            >  2:2779421/1‑150 (MQ=255)
                                                                                                         aaTGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAAt                                           <  2:2718824/150‑1 (MQ=255)
                                                                                                          aTGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAAtt                                          <  2:420419/150‑1 (MQ=255)
                                                                                                          aTGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAAtt                                          <  1:3481903/150‑1 (MQ=255)
                                                                                                           tGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTc                                         <  1:1448390/150‑1 (MQ=255)
                                                                                                            gTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTcc                                        <  1:1991327/150‑1 (MQ=255)
                                                                                                                 gATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGa                                   <  2:3886677/150‑1 (MQ=255)
                                                                                                                   ttAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATg                                 >  2:3065553/1‑150 (MQ=255)
                                                                                                                   ttAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATg                                 >  2:2083468/1‑150 (MQ=255)
                                                                                                                    tAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGt                                <  2:1217256/150‑1 (MQ=255)
                                                                                                                    tAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGt                                >  2:448312/1‑150 (MQ=255)
                                                                                                                      aaGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTg                              <  1:2691373/150‑1 (MQ=255)
                                                                                                                      aaGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTg                              >  1:4304910/1‑150 (MQ=255)
                                                                                                                      aaGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTg                              >  2:4137368/1‑150 (MQ=255)
                                                                                                                       aGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATACCATCGATGttca                             >  2:4011379/1‑148 (MQ=255)
                                                                                                                        gcgcCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaa                            <  1:2525521/150‑1 (MQ=255)
                                                                                                                         cgcCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaa                           <  2:1979243/150‑1 (MQ=255)
                                                                                                                           ccTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaaaa                         >  1:1258808/1‑150 (MQ=255)
                                                                                                                                     gCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTagca               <  1:4466355/150‑1 (MQ=255)
                                                                                                                                     gCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTagca               >  1:7438/1‑150 (MQ=255)
                                                                                                                                      cGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTagcag              >  2:540091/1‑150 (MQ=255)
                                                                                                                                       ggTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTagcagc             >  2:2136962/1‑150 (MQ=255)
                                                                                                                                       ggTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTagcagc             >  2:4233704/1‑150 (MQ=255)
                                                                                                                                        gtgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCg            >  2:1097875/1‑150 (MQ=255)
                                                                                                                                         tgtACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGACTTCCATCGATGTTGAAAAAAATATTAGCAGCGt           >  1:1347488/1‑150 (MQ=255)
                                                                                                                                         tgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGt           <  1:2779421/150‑1 (MQ=255)
                                                                                                                                         tgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGt           <  2:690150/150‑1 (MQ=255)
                                                                                                                                          gtACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGtt          <  1:4716395/150‑1 (MQ=255)
                                                                                                                                          gtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGtt          >  1:2157067/1‑150 (MQ=255)
                                                                                                                                          gtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGAGGTTGAAAAAAATATTAGCAGCGtt          >  1:4558804/1‑150 (MQ=255)
                                                                                                                                              gCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGc      >  2:2393491/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGAC  >  CP009273/557154‑557449

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: ‑