Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,375,260 G→T 100% D387Y (GAT→TAT)  degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,375,2600GT99.6% 953.7 / ‑4.1 269D387Y (GAT→TAT) degQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (125/142);  total (125/144)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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 271 total aligned reads displayed.

GGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAA  >  CP009273/3375113‑3375406
                                                                                                                                                   |                                                                                                                                                  
ggCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAatat                                                                                                                                                  >  1:1216580/1‑150 (MQ=255)
 gCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATg                                                                                                                                                 <  1:902946/150‑1 (MQ=255)
  cACGAAAGTGAAGCTTGGCCTACTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATgg                                                                                                                                                <  2:2370633/150‑1 (MQ=255)
   aCGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggc                                                                                                                                               <  1:1250879/150‑1 (MQ=255)
   aCGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggc                                                                                                                                               <  1:2860101/150‑1 (MQ=255)
    cGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggcg                                                                                                                                              <  2:2398911/150‑1 (MQ=255)
     gAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggcgg                                                                                                                                             >  1:1405564/1‑150 (MQ=255)
      aaaGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggcggc                                                                                                                                            >  2:1579247/1‑150 (MQ=255)
      aaaGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggcggc                                                                                                                                            >  2:457666/1‑150 (MQ=255)
       aaGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCa                                                                                                                                           >  1:100205/1‑150 (MQ=255)
        aGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCaa                                                                                                                                          <  2:339439/150‑1 (MQ=255)
        aGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCaa                                                                                                                                          <  2:357067/150‑1 (MQ=255)
            aaGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGt                                                                                                                                      <  2:1998240/150‑1 (MQ=255)
              gCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTAt                                                                                                                                    <  2:3062451/150‑1 (MQ=255)
                ttGGCCTGCTGCGTAACGTCAAACCACTGGAAGTAGAAGTGACGCTCGAGACCATCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTa                                                                                                                                  <  1:931047/150‑1 (MQ=255)
                 tGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTaa                                                                                                                                 <  1:3999904/150‑1 (MQ=255)
                   gCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTaaaa                                                                                                                               <  1:3663984/150‑1 (MQ=255)
                    ccTGCTTCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCACGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGAAGGTCAGCTAAAATATGGCGGCAAAGGTCTTAAAAt                                                                                                                              >  2:3583171/1‑150 (MQ=255)
                    ccTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAAt                                                                                                                              >  2:2584981/1‑150 (MQ=255)
                    ccTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAAt                                                                                                                              >  1:1813450/1‑150 (MQ=255)
                     ctgctgCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATc                                                                                                                             <  2:33739/150‑1 (MQ=255)
                       gctgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGa                                                                                                                           >  2:2967022/1‑150 (MQ=255)
                        ctgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGAt                                                                                                                          <  2:793482/150‑1 (MQ=255)
                        ctgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGAt                                                                                                                          >  1:2921202/1‑150 (MQ=255)
                         tgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATg                                                                                                                         <  1:2518371/150‑1 (MQ=255)
                         tgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATg                                                                                                                         <  1:1653101/150‑1 (MQ=255)
                            gTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAg                                                                                                                      >  2:552473/1‑150 (MQ=255)
                            gTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAg                                                                                                                      >  2:251754/1‑150 (MQ=255)
                               aCGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTg                                                                                                                   >  2:657488/1‑150 (MQ=255)
                                   cAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCCGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCaa                                                                                                               >  1:105570/1‑150 (MQ=255)
                                   cAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCaa                                                                                                               >  2:2457685/1‑150 (MQ=255)
                                      aCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAgg                                                                                                            <  1:3035589/150‑1 (MQ=255)
                                         aCTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAg                                                                                                         >  2:3841212/1‑150 (MQ=255)
                                         aCTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAg                                                                                                         <  1:2123313/150‑1 (MQ=255)
                                          cTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGc                                                                                                        <  2:3991728/150‑1 (MQ=255)
                                             gAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCa                                                                                                     <  2:3729056/150‑1 (MQ=255)
                                                gTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgct                                                                                                  <  1:681386/150‑1 (MQ=255)
                                                    aaGTGCCGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                                                                                              <  2:2142409/150‑1 (MQ=255)
                                                    aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                                                                                              <  2:3882603/150‑1 (MQ=255)
                                                    aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                                                                                              >  1:2513834/1‑150 (MQ=255)
                                                      gTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAg                                                                                            >  1:3794561/1‑150 (MQ=255)
                                                        gACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggc                                                                                          <  1:3122940/150‑1 (MQ=255)
                                                        gACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggc                                                                                          >  1:829061/1‑150 (MQ=255)
                                                        gACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggc                                                                                          >  2:3672029/1‑150 (MQ=255)
                                                        gACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggc                                                                                          <  2:2283658/150‑1 (MQ=255)
                                                         aCGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggct                                                                                         <  2:3317144/150‑1 (MQ=255)
                                                         aCGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggct                                                                                         <  2:1990945/150‑1 (MQ=255)
                                                          cGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctg                                                                                        <  1:2657099/150‑1 (MQ=255)
                                                           gCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctgg                                                                                       <  1:2692980/150‑1 (MQ=255)
                                                            cTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctggc                                                                                      <  2:593871/150‑1 (MQ=255)
                                                              cGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCtt                                                                                    <  1:3381087/150‑1 (MQ=255)
                                                              cGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCtt                                                                                    <  2:2215205/150‑1 (MQ=255)
                                                               gATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTg                                                                                   >  2:1643888/1‑150 (MQ=255)
                                                                 tACCAGCACCACTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGAATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCa                                                                                 <  1:310266/150‑1 (MQ=255)
                                                                       cACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGa                                                                           <  1:1030694/150‑1 (MQ=255)
                                                                       cACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGa                                                                           <  1:2890116/150‑1 (MQ=255)
                                                                         ccTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACg                                                                         >  1:2419455/1‑150 (MQ=255)
                                                                            ctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtg                                                                      <  1:4143703/150‑1 (MQ=255)
                                                                             ttCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtgt                                                                     >  2:2694934/1‑150 (MQ=255)
                                                                                 tcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATc                                                                 >  2:1914795/1‑150 (MQ=255)
                                                                                   ggCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAt                                                               <  2:2767033/150‑1 (MQ=255)
                                                                                   ggCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAt                                                               >  1:3986409/1‑150 (MQ=255)
                                                                                    gCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAtt                                                              <  1:2669608/150‑1 (MQ=255)
                                                                                      cAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTgg                                                            <  1:3592615/150‑1 (MQ=255)
                                                                                      cAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTgg                                                            <  1:3436206/150‑1 (MQ=255)
                                                                                       aGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGc                                                           >  2:1282945/1‑150 (MQ=255)
                                                                                           cagAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCa                                                       <  2:1723120/148‑1 (MQ=255)
                                                                                             gAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAAc                                                     <  1:2021449/150‑1 (MQ=255)
                                                                                              aacggATCTCGCCTGCGCTGGACGGTGCATCGTTGAGCGTTGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAAcc                                                    <  2:1014350/146‑1 (MQ=255)
                                                                                                aTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcgc                                                  >  2:710669/1‑150 (MQ=255)
                                                                                                       cGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCggg                                           >  1:1565285/1‑150 (MQ=255)
                                                                                                        gCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGt                                          <  1:1342365/150‑1 (MQ=255)
                                                                                                          cAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGa                                        >  1:1714401/1‑150 (MQ=255)
                                                                                                          cAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGa                                        >  1:554612/1‑150 (MQ=255)
                                                                                                           aGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGaa                                       >  2:2196316/1‑150 (MQ=255)
                                                                                                             cgcTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                     >  2:584708/1‑150 (MQ=255)
                                                                                                              gcTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                    >  2:2545024/1‑150 (MQ=255)
                                                                                                              gcTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                    <  2:1835824/150‑1 (MQ=255)
                                                                                                                   aaGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTg                               >  2:1219563/1‑150 (MQ=255)
                                                                                                                   aaGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTg                               <  1:423399/150‑1 (MQ=255)
                                                                                                                   aaGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTg                               <  1:1897833/150‑1 (MQ=255)
                                                                                                                    aGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                              <  2:4119199/150‑1 (MQ=255)
                                                                                                                        gCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGaa                          <  1:30191/150‑1 (MQ=255)
                                                                                                                         cAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGaaa                         >  1:561209/1‑150 (MQ=255)
                                                                                                                              ttGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGt                    <  2:2120372/150‑1 (MQ=255)
                                                                                                                              ttGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGt                    >  2:3083872/1‑150 (MQ=255)
                                                                                                                               tGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTa                   >  1:1266952/1‑150 (MQ=255)
                                                                                                                                 aGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                 >  1:1437579/1‑150 (MQ=255)
                                                                                                                                      tGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCt            <  1:2263795/150‑1 (MQ=255)
                                                                                                                                      tGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCt            <  2:1544521/150‑1 (MQ=255)
                                                                                                                                      tGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATGGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCt            <  2:3415399/150‑1 (MQ=255)
                                                                                                                                      tGGTCAGCTAAAAGAGGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTTATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTACAGTGCt            <  2:1961659/150‑1 (MQ=255)
                                                                                                                                         tCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggc         >  2:4144864/1‑150 (MQ=255)
                                                                                                                                           aGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggcgg       >  2:361095/1‑150 (MQ=255)
                                                                                                                                             cTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCa     >  2:2548264/1‑150 (MQ=255)
                                                                                                                                              tAAAATATGGCGGCAAATGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa    >  1:3250655/1‑150 (MQ=255)
                                                                                                                                               aaaaTATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAATACGATGTGATCATTGGCTTCAACAGCGATCGGTAGAACTCGATTGCTGAAATGAGTAAAGATCTGCCGGCaaa   >  2:1226836/1‑150 (MQ=255)
                                                                                                                                               aaaaTATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaaa   <  2:856657/150‑1 (MQ=255)
                                                                                                                                               aaaaTATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaaa   <  2:938714/150‑1 (MQ=255)
                                                                                                                                                aaaTATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaaaa  >  1:2235647/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAA  >  CP009273/3375113‑3375406

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: