Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,373,985 C→A 100% intergenic (+37/‑117) yhcB → / → degQ DUF1043 family inner membrane‑anchored protein/serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,373,9850CA98.6% 951.5 / 2.4 293intergenic (+37/‑117)yhcB/degQDUF1043 family inner membrane‑anchored protein/serine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base A (118/169);  total (118/175)
Fisher's exact test for biased strand distribution p-value = 8.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 303 total aligned reads displayed.

AATCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTGCA  >  CP009273/3373838‑3374131
                                                                                                                                                   |                                                                                                                                                  
aaTCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTtatc                                                                                                                                                  <  2:2171408/150‑1 (MQ=255)
aaTCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAAGTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTtatc                                                                                                                                                  >  2:1585493/1‑150 (MQ=255)
    gTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccc                                                                                                                                              >  2:601393/1‑150 (MQ=255)
     tCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccct                                                                                                                                             <  1:3776638/150‑1 (MQ=255)
       tGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccctcc                                                                                                                                           >  1:2297518/1‑150 (MQ=255)
        ggCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCtccctccc                                                                                                                                          <  2:77654/150‑1 (MQ=255)
           agagTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGa                                                                                                                                       <  1:390924/150‑1 (MQ=255)
           agagTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGa                                                                                                                                       >  1:1424518/1‑150 (MQ=255)
            gagTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGAc                                                                                                                                      <  1:2438082/150‑1 (MQ=255)
             agTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACt                                                                                                                                     <  1:2889802/150‑1 (MQ=255)
                 tGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCa                                                                                                                                 >  1:1095550/1‑150 (MQ=255)
                    aGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCAttt                                                                                                                              >  2:1235300/1‑150 (MQ=255)
                    aGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCAttt                                                                                                                              >  1:2135401/1‑150 (MQ=255)
                       cAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAAt                                                                                                                           >  1:1106805/1‑150 (MQ=255)
                           aaCGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTgg                                                                                                                       <  1:3219489/150‑1 (MQ=255)
                                tCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGtct                                                                                                                  <  2:377445/150‑1 (MQ=255)
                                 cAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGtctc                                                                                                                 >  1:3815909/1‑150 (MQ=255)
                                  aGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCNCCTCTTATCNCCCTCCCCGACTATCATTTAATCTGGTGTCTCa                                                                                                                >  1:2992640/1‑150 (MQ=255)
                                    gCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCAtt                                                                                                              >  1:1740429/1‑150 (MQ=255)
                                      aCCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGt                                                                                                            >  1:3871821/1‑150 (MQ=255)
                                         ggTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAg                                                                                                         <  1:131338/150‑1 (MQ=255)
                                           tGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGcc                                                                                                       >  2:3683960/1‑150 (MQ=255)
                                             cAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGt                                                                                                     <  1:810355/150‑1 (MQ=255)
                                               gATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCt                                                                                                   <  2:870439/150‑1 (MQ=255)
                                               gATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCt                                                                                                   >  2:2275725/1‑150 (MQ=255)
                                               gATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCt                                                                                                   <  2:2774020/150‑1 (MQ=255)
                                                 tGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGa                                                                                                 >  2:2348223/1‑150 (MQ=255)
                                                   ccTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGaaa                                                                                               <  2:699572/150‑1 (MQ=255)
                                                     tCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAAt                                                                                             >  2:1777309/1‑150 (MQ=255)
                                                     tCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAAt                                                                                             >  1:646424/1‑150 (MQ=255)
                                                       gcgACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTc                                                                                           <  2:2004906/150‑1 (MQ=255)
                                                        cgACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCa                                                                                          <  2:901614/150‑1 (MQ=255)
                                                        cgACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCAAAGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCa                                                                                          >  2:3468008/1‑150 (MQ=255)
                                                           cTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATa                                                                                       >  2:510968/1‑150 (MQ=255)
                                                             aTTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAAc                                                                                     >  1:3272258/1‑150 (MQ=255)
                                                              ttCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACa                                                                                    <  1:3695465/150‑1 (MQ=255)
                                                              ttCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACa                                                                                    <  1:2124418/150‑1 (MQ=255)
                                                              ttCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACa                                                                                    >  1:250653/1‑150 (MQ=255)
                                                                cTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATc                                                                                  <  1:719823/150‑1 (MQ=255)
                                                                       cgcATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGt                                                                           >  1:3458269/1‑150 (MQ=255)
                                                                       cgcATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGt                                                                           >  1:2780572/1‑150 (MQ=255)
                                                                          aTCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGtttt                                                                        <  2:3191026/150‑1 (MQ=255)
                                                                          aTCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGtttt                                                                        <  2:3693407/150‑1 (MQ=255)
                                                                           tCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTg                                                                       >  1:3254726/1‑150 (MQ=255)
                                                                           tCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTg                                                                       <  1:3890220/150‑1 (MQ=255)
                                                                             cGGCCTGCTGCGTACTGGCGCGAAGCGCGTCTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGaa                                                                     >  2:643793/1‑150 (MQ=255)
                                                                              ggCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCTCCCTCCTCTTATCTCCCTCCCCGAATATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAt                                                                    <  1:1753970/150‑1 (MQ=255)
                                                                              ggCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAt                                                                    <  1:3308307/150‑1 (MQ=255)
                                                                               gCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtc                                                                   <  2:3030704/150‑1 (MQ=255)
                                                                                 ctgctgCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtctc                                                                 <  1:2192403/150‑1 (MQ=255)
                                                                                 ctgctgCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtctc                                                                 >  1:900653/1‑150 (MQ=255)
                                                                                   gctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCtt                                                               >  2:4148153/1‑150 (MQ=255)
                                                                                    ctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCttt                                                              <  2:2446694/150‑1 (MQ=255)
                                                                                    ctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCttt                                                              <  1:3297415/150‑1 (MQ=255)
                                                                                      gcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTc                                                            <  1:1064577/150‑1 (MQ=255)
                                                                                      gcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTc                                                            <  1:1021072/150‑1 (MQ=255)
                                                                                        gTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCtt                                                          <  1:3628168/150‑1 (MQ=255)
                                                                                          aCTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTAt                                                        <  1:2666132/150‑1 (MQ=255)
                                                                                           cTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATc                                                       >  1:3200047/1‑150 (MQ=255)
                                                                                            tGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCa                                                      >  1:2071889/1‑150 (MQ=255)
                                                                                             ggCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCAt                                                     <  1:3725026/150‑1 (MQ=255)
                                                                                              gcgcGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGTGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCAtt                                                    <  1:547177/150‑1 (MQ=255)
                                                                                              gcgcGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCAtt                                                    <  2:2071017/150‑1 (MQ=255)
                                                                                               cgcgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTc                                                   >  1:1145738/1‑150 (MQ=255)
                                                                                                gcgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCa                                                  <  2:3812087/150‑1 (MQ=255)
                                                                                                 cgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAg                                                 >  2:3308039/1‑150 (MQ=255)
                                                                                                 cgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAg                                                 <  1:794028/150‑1 (MQ=255)
                                                                                                  gAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTATGAATCTCTTTTCTTATCATTCAgg                                                >  1:4047044/1‑150 (MQ=255)
                                                                                                    aGCGCGTCTAAGTTATTTTGCGGGCGCAGCCTTTGCGCCCTCCTTTTTTCTCCCTCCCCTACGATCTTTTAATCTGGTGTCTGATTGTTAGCCGTCTGAACATTCAATAACATCATACTGTTTTGAATCTCTTTTCTTATCATTCAGGTa                                              <  2:3613560/150‑1 (MQ=255)
                                                                                                     gcgcGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTAc                                             <  2:4041658/150‑1 (MQ=255)
                                                                                                     gcgcGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTAc                                             <  2:1943165/150‑1 (MQ=255)
                                                                                                        cgACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACgag                                          <  2:1185177/150‑1 (MQ=255)
                                                                                                         gACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACgaga                                         <  2:3401873/150‑1 (MQ=255)
                                                                                                          aCTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGagag                                        >  2:1412210/1‑150 (MQ=255)
                                                                                                              atttatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAgg                                    >  1:3067953/1‑150 (MQ=255)
                                                                                                              atttatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAgg                                    >  1:4089076/1‑150 (MQ=255)
                                                                                                                 tatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaat                                 >  1:476850/1‑150 (MQ=255)
                                                                                                                   tttttCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaataa                               <  1:2226585/150‑1 (MQ=255)
                                                                                                                    ttttCGGGCGCAGCCATTGGGCCCTGCTCTTATCTCCCTCCCCGACTATAATTTAATCTGGTGTCTAATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTTCGAGAGCAGGaataat                              <  2:439931/150‑1 (MQ=255)
                                                                                                                    ttttCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaataat                              <  2:4160013/150‑1 (MQ=255)
                                                                                                                            cgcAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaac                      <  1:4165054/150‑1 (MQ=255)
                                                                                                                              cAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaa                    >  1:3102599/1‑150 (MQ=255)
                                                                                                                               aGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaa                   >  1:197019/1‑150 (MQ=255)
                                                                                                                               aGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaa                   <  2:719801/150‑1 (MQ=255)
                                                                                                                               aGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTGTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaa                   >  2:1522585/1‑150 (MQ=255)
                                                                                                                                gCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaac                  >  1:3720078/1‑150 (MQ=255)
                                                                                                                                 ccATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAAcc                 >  1:344096/1‑150 (MQ=255)
                                                                                                                                 ccATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAAcc                 >  1:3634582/1‑150 (MQ=255)
                                                                                                                                  cATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAAccc                <  2:3738766/150‑1 (MQ=255)
                                                                                                                                  cATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAAccc                <  2:1379238/150‑1 (MQ=255)
                                                                                                                                      gcgcCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCt            <  1:1977569/150‑1 (MQ=255)
                                                                                                                                      gcgcCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCt            <  1:969563/150‑1 (MQ=255)
                                                                                                                                       cgcCCTCCTCTTAGCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTg           <  2:1055045/150‑1 (MQ=255)
                                                                                                                                        gcCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGt          >  2:3900112/1‑150 (MQ=255)
                                                                                                                                          cctcctCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTg        <  1:2576009/150‑1 (MQ=255)
                                                                                                                                          cctcctCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTg        >  1:2604183/1‑150 (MQ=255)
                                                                                                                                          cctcctCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTg        >  2:3103987/1‑150 (MQ=255)
                                                                                                                                             cctcTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGt     <  2:1106805/150‑1 (MQ=255)
                                                                                                                                               tcttatCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCGGTTGAGTGc   >  2:1276991/1‑150 (MQ=255)
                                                                                                                                                cttatCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTGCa  <  1:2399621/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTGCA  >  CP009273/3373838‑3374131

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: