Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,915,998 C→A 29.3% A116S (GCT→TCT)  yqcA ← short‑chain flavodoxin, FMN‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,915,9980CA29.3% 655.5 / 242.2 320A116S (GCT→TCT) yqcAshort‑chain flavodoxin, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (109/117);  new base A (41/53);  total (150/170)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 325 total aligned reads displayed.

GATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGCACAATGCT  >  CP009273/2915851‑2916145
                                                                                                                                                   |                                                                                                                                                   
  ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGct                                                                                                                                                  <  2:1712074/149‑1 (MQ=255)
   tcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGctc                                                                                                                                                 <  1:2415348/149‑1 (MQ=255)
      gCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGt                                                                                                                                              <  1:3034674/149‑1 (MQ=255)
      gCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGt                                                                                                                                              <  2:3470650/149‑1 (MQ=255)
      gCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGt                                                                                                                                              >  2:886563/1‑149 (MQ=255)
      gCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGt                                                                                                                                              <  2:711990/149‑1 (MQ=255)
         ccGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTcc                                                                                                                                           <  2:642834/149‑1 (MQ=255)
          cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCt                                                                                                                                          >  1:3326564/1‑149 (MQ=255)
            cTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGc                                                                                                                                        >  1:927718/1‑149 (MQ=255)
            cTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGc                                                                                                                                        <  1:647523/149‑1 (MQ=255)
             tAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGca                                                                                                                                       >  2:2870284/1‑149 (MQ=255)
             tAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGca                                                                                                                                       >  1:1706409/1‑149 (MQ=255)
              aaaaCGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGcaa                                                                                                                                      >  1:3760882/1‑149 (MQ=255)
                 aCGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGcaacaa                                                                                                                                   >  1:4067623/1‑149 (MQ=255)
                  cGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAg                                                                                                                                  <  2:4770078/149‑1 (MQ=255)
                  cGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAg                                                                                                                                  >  2:69232/1‑149 (MQ=255)
                  cGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAg                                                                                                                                  >  2:1853149/1‑149 (MQ=255)
                     cATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCATCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCa                                                                                                                               <  1:3937116/149‑1 (MQ=255)
                      aTTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCAt                                                                                                                              >  2:1244192/1‑149 (MQ=255)
                      aTTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCAt                                                                                                                              >  2:3849201/1‑149 (MQ=255)
                       tttACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATc                                                                                                                             >  1:2514259/1‑149 (MQ=255)
                        ttACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCg                                                                                                                            >  2:1385037/1‑149 (MQ=255)
                        ttACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCg                                                                                                                            <  2:3264898/149‑1 (MQ=255)
                        ttACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCg                                                                                                                            >  1:1800320/1‑149 (MQ=255)
                          aCCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGaa                                                                                                                          <  1:3272447/149‑1 (MQ=255)
                              ctcgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGt                                                                                                                      <  1:189747/149‑1 (MQ=255)
                               tcgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGtt                                                                                                                     <  1:2013754/149‑1 (MQ=255)
                               tcgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGtt                                                                                                                     <  1:4465771/149‑1 (MQ=255)
                                cgctGAACATATCAGGACAACCGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGttt                                                                                                                    <  2:844419/149‑1 (MQ=255)
                                  ctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgc                                                                                                                  >  1:1358479/1‑149 (MQ=255)
                                  ctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgc                                                                                                                  <  2:1119074/149‑1 (MQ=255)
                                   tGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgcc                                                                                                                 >  2:1809015/1‑149 (MQ=255)
                                   tGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgcc                                                                                                                 >  2:1424688/1‑149 (MQ=255)
                                   tGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgcc                                                                                                                 >  1:4733829/1‑149 (MQ=255)
                                      aCATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGAGTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgccgcc                                                                                                              <  2:4625302/149‑1 (MQ=255)
                                       cATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCa                                                                                                             >  2:4043188/1‑149 (MQ=255)
                                       cATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCa                                                                                                             >  1:1615298/1‑149 (MQ=255)
                                       cATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCa                                                                                                             <  1:867284/149‑1 (MQ=255)
                                        atatCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGTGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCAt                                                                                                            >  1:1026925/1‑149 (MQ=255)
                                             aGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAg                                                                                                       <  1:1159207/149‑1 (MQ=255)
                                             aGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAg                                                                                                       >  1:819551/1‑149 (MQ=255)
                                                acaacaGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGaaa                                                                                                    <  2:3004354/149‑1 (MQ=255)
                                                  aacaGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAAtt                                                                                                  <  1:4121959/149‑1 (MQ=255)
                                                   acaGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTc                                                                                                 >  2:1036664/1‑149 (MQ=255)
                                                   acaGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTc                                                                                                 >  1:4376641/1‑149 (MQ=255)
                                                           ccccACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAACTTCACATAact                                                                                         >  2:4353599/1‑149 (MQ=255)
                                                           ccccACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAact                                                                                         <  2:2247028/149‑1 (MQ=255)
                                                             ccACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGGTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAactac                                                                                       <  2:3094398/149‑1 (MQ=255)
                                                              cACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAactact                                                                                      <  1:1559524/149‑1 (MQ=255)
                                                               aCTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTg                                                                                     <  2:2484220/149‑1 (MQ=255)
                                                                cggTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGt                                                                                    <  2:1175202/147‑1 (MQ=255)
                                                                cTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGt                                                                                    >  1:2314700/1‑149 (MQ=255)
                                                                 tGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACGGTc                                                                                   <  2:3897666/149‑1 (MQ=255)
                                                                  gTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCg                                                                                  >  2:2740575/1‑149 (MQ=255)
                                                                  gTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCg                                                                                  >  2:3942714/1‑149 (MQ=255)
                                                                       aCCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAg                                                                             >  1:2867447/1‑149 (MQ=255)
                                                                        cccACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAgc                                                                            >  1:2679682/1‑149 (MQ=255)
                                                                          cACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAgcgc                                                                          >  2:632649/1‑149 (MQ=255)
                                                                           aCGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCa                                                                         >  2:1281223/1‑149 (MQ=255)
                                                                           aCGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCa                                                                         >  2:3189646/1‑149 (MQ=255)
                                                                           aCGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCa                                                                         >  1:2958158/1‑149 (MQ=255)
                                                                            cGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCaa                                                                        <  2:2081846/149‑1 (MQ=255)
                                                                            cGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCaa                                                                        <  2:3242119/149‑1 (MQ=255)
                                                                              gATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATc                                                                      <  1:106268/149‑1 (MQ=255)
                                                                              gATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATc                                                                      <  1:2810073/149‑1 (MQ=255)
                                                                               aTTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCa                                                                     >  1:3544260/1‑149 (MQ=255)
                                                                                  tGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGc                                                                  <  2:1603067/149‑1 (MQ=255)
                                                                                   gATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGcc                                                                 >  2:4584343/1‑149 (MQ=255)
                                                                                          tttCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGc                                                          >  1:1050341/1‑149 (MQ=255)
                                                                                          tttCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGc                                                          >  1:3828991/1‑149 (MQ=255)
                                                                                              cggttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGAt                                                      >  2:1870778/1‑149 (MQ=255)
                                                                                                 ttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCg                                                   >  2:3601566/1‑149 (MQ=255)
                                                                                                 ttccggGTACTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCg                                                   >  2:1663076/1‑149 (MQ=255)
                                                                                                    cggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGGTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCt                                                <  2:4690816/149‑1 (MQ=255)
                                                                                                        ttCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGaa                                            >  1:499926/1‑149 (MQ=255)
                                                                                                        ttCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGaa                                            >  1:1274172/1‑149 (MQ=255)
                                                                                                          ctcACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGc                                          >  2:251944/1‑149 (MQ=255)
                                                                                                           tcACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGcc                                         <  2:2268947/149‑1 (MQ=255)
                                                                                                               tGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGa                                     >  2:2595708/1‑149 (MQ=255)
                                                                                                                ggCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGAc                                    <  2:1085042/149‑1 (MQ=255)
                                                                                                                 gCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACt                                   <  2:3155818/149‑1 (MQ=255)
                                                                                                                  cGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACTCAGATTCGGCTGGAAGCCCAGACTa                                  <  2:4253493/149‑1 (MQ=255)
                                                                                                                   gTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTAt                                 >  2:3424727/1‑149 (MQ=255)
                                                                                                                    tcaatcaaCAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATc                                >  1:463187/1‑149 (MQ=255)
                                                                                                                    tcaatcaaCAGCATATCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATc                                >  1:2238998/1‑149 (MQ=255)
                                                                                                                      aatcaaCAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCtt                              >  2:3065410/1‑149 (MQ=255)
                                                                                                                        tcaaCAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTg                            <  1:1373045/149‑1 (MQ=255)
                                                                                                                          aaCAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGAt                          >  1:4535648/1‑149 (MQ=255)
                                                                                                                           ccAGCATTTCGCCAACCCGCTCAGCGCTCTGTTCCTGCAACAAGGCAGCGAACTGGTTGCCGCCATAGCAGAACTTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGAGTCGGCTGGACGCCCAGACTATCTTTGAta                         >  1:3079102/2‑148 (MQ=255)
                                                                                                                             aGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTcc                       >  1:4249109/1‑149 (MQ=255)
                                                                                                                             aGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTcc                       >  1:3592255/1‑149 (MQ=255)
                                                                                                                               cATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCt                     <  2:2672550/149‑1 (MQ=255)
                                                                                                                                aTTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTg                    <  1:4042808/149‑1 (MQ=255)
                                                                                                                                aTTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTg                    <  2:2503365/149‑1 (MQ=255)
                                                                                                                                 tttCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACCTAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGa                   <  1:1294781/149‑1 (MQ=255)
                                                                                                                                 tttCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGa                   <  2:2992508/149‑1 (MQ=255)
                                                                                                                                     gCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAaga               >  1:545874/1‑149 (MQ=255)
                                                                                                                                      ccAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAagag              <  2:2297307/149‑1 (MQ=255)
                                                                                                                                        aaCCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTg            >  2:2766627/1‑149 (MQ=255)
                                                                                                                                         aCCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTgg           >  2:3664858/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGCACAATGCT  >  CP009273/2915851‑2916145

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: