Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,915,998 C→A 100% 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,9980CA98.5% 858.0 / 3.0 263A116S (GCT→TCT) yqcAshort‑chain flavodoxin, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (118/141);  total (120/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 271 total aligned reads displayed.

ATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGCACAATG  >  CP009273/2915852‑2916143
                                                                                                                                                  |                                                                                                                                                 
attcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAgagc                                                                                                                                                 <  1:3457910/149‑1 (MQ=255)
 ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGct                                                                                                                                                >  1:828483/1‑149 (MQ=255)
 ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGct                                                                                                                                                >  2:2775710/1‑149 (MQ=255)
 ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGct                                                                                                                                                >  2:1217368/1‑149 (MQ=255)
    aGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTg                                                                                                                                             >  1:361313/1‑149 (MQ=255)
     gCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGt                                                                                                                                            >  1:1056841/1‑149 (MQ=255)
      cACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGtt                                                                                                                                           <  1:3553120/149‑1 (MQ=255)
      cACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGtt                                                                                                                                           <  1:2652410/149‑1 (MQ=255)
           cTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGc                                                                                                                                      >  2:3280000/1‑149 (MQ=255)
             aaaaCGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGcaa                                                                                                                                    >  2:590857/1‑149 (MQ=255)
              aaaCGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGcaac                                                                                                                                   >  2:3497075/1‑149 (MQ=255)
               aaCGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGcaaca                                                                                                                                  >  2:4024579/1‑149 (MQ=255)
                aCGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGcaacaa                                                                                                                                 <  1:547246/149‑1 (MQ=255)
                   aCATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGc                                                                                                                              >  2:150795/1‑149 (MQ=255)
                   aCATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGc                                                                                                                              <  1:2626927/149‑1 (MQ=255)
                          ccGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAAc                                                                                                                       <  2:1024035/149‑1 (MQ=255)
                          ccGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAAc                                                                                                                       <  1:3037280/149‑1 (MQ=255)
                              tcgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGtt                                                                                                                   <  2:3736449/149‑1 (MQ=255)
                               cgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGttt                                                                                                                  <  1:3597311/149‑1 (MQ=255)
                                 ctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgc                                                                                                                <  2:2080913/149‑1 (MQ=255)
                                     aCATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgccgcc                                                                                                            >  1:2326168/1‑149 (MQ=255)
                                        tatCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCAtt                                                                                                         <  1:1115186/149‑1 (MQ=255)
                                        tatCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCAtt                                                                                                         <  1:2939605/149‑1 (MQ=255)
                                           cAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCa                                                                                                      <  2:1670856/149‑1 (MQ=255)
                                             ggACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGa                                                                                                    >  2:3987671/1‑149 (MQ=255)
                                             ggACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGa                                                                                                    >  2:773830/1‑149 (MQ=255)
                                              gACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGaa                                                                                                   <  1:1299429/149‑1 (MQ=255)
                                               acaacaGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGaaa                                                                                                  >  2:1458127/1‑149 (MQ=255)
                                               acaacaGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGaaa                                                                                                  <  1:326596/149‑1 (MQ=255)
                                               acaacaGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGaaa                                                                                                  <  1:1703660/149‑1 (MQ=255)
                                                   caGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTca                                                                                              <  1:237198/149‑1 (MQ=255)
                                                         gCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAac                                                                                        <  2:1429086/149‑1 (MQ=255)
                                                         gCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAac                                                                                        <  1:1673523/149‑1 (MQ=255)
                                                         gCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAac                                                                                        <  2:3189650/149‑1 (MQ=255)
                                                           cccACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAacta                                                                                      <  1:6022/149‑1 (MQ=255)
                                                           cccACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAacta                                                                                      <  1:1622186/149‑1 (MQ=255)
                                                              aCTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTg                                                                                   >  2:512974/1‑149 (MQ=255)
                                                              aCTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTg                                                                                   >  2:2466853/1‑149 (MQ=255)
                                                              aCTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTg                                                                                   <  1:3091639/149‑1 (MQ=255)
                                                              aCTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTg                                                                                   >  2:2315740/1‑149 (MQ=255)
                                                                tGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTc                                                                                 >  1:1710029/1‑149 (MQ=255)
                                                                     gACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGa                                                                            >  2:1277088/1‑149 (MQ=255)
                                                                     gACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGa                                                                            >  2:3326743/1‑149 (MQ=255)
                                                                      aCCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAg                                                                           <  2:370561/149‑1 (MQ=255)
                                                                         cACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAgcgc                                                                        <  1:644020/149‑1 (MQ=255)
                                                                          aCGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTGCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCa                                                                       >  1:3423051/1‑149 (MQ=255)
                                                                           cGGATTTGAGTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCaa                                                                      <  1:3820827/149‑1 (MQ=255)
                                                                            ggATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAAt                                                                     <  2:2859827/149‑1 (MQ=255)
                                                                               tttGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCAc                                                                  >  2:448888/1‑149 (MQ=255)
                                                                                 tGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGc                                                                <  2:3493879/149‑1 (MQ=255)
                                                                                 tGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGc                                                                <  2:1435496/149‑1 (MQ=255)
                                                                                  gATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGcc                                                               >  2:3159672/1‑149 (MQ=255)
                                                                                   aTTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCa                                                              >  1:2708847/1‑149 (MQ=255)
                                                                                   aTTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCa                                                              >  1:813653/1‑149 (MQ=255)
                                                                                       cGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGAATCGAACTGTTTGCCGCCATTGCATAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAAc                                                          <  1:2786618/149‑1 (MQ=255)
                                                                                         tttCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGc                                                        >  1:1903739/1‑149 (MQ=255)
                                                                                         tttCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGc                                                        >  2:184456/1‑149 (MQ=255)
                                                                                         tttCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGc                                                        >  2:3304190/1‑149 (MQ=255)
                                                                                           tccggttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAg                                                      >  2:2554643/1‑149 (MQ=255)
                                                                                              ggttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGAtt                                                   >  2:1646567/1‑149 (MQ=255)
                                                                                               gttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTc                                                  <  1:2189884/149‑1 (MQ=255)
                                                                                                ttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCg                                                 <  2:3088015/149‑1 (MQ=255)
                                                                                                ttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCg                                                 <  2:551295/149‑1 (MQ=255)
                                                                                                 tccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCgg                                                <  2:3412118/149‑1 (MQ=255)
                                                                                                  ccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGc                                               <  2:745223/149‑1 (MQ=255)
                                                                                                   cggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCt                                              <  1:67694/149‑1 (MQ=255)
                                                                                                     ggTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTgg                                            <  2:2641026/149‑1 (MQ=255)
                                                                                                     ggTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTgg                                            <  2:2288297/149‑1 (MQ=255)
                                                                                                         ctcACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGc                                        >  1:1004109/1‑149 (MQ=255)
                                                                                                          tcACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGcc                                       <  1:1466411/149‑1 (MQ=255)
                                                                                                             cTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAg                                    >  2:298440/1‑149 (MQ=255)
                                                                                                                gCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACt                                 >  1:1318268/1‑149 (MQ=255)
                                                                                                                 cGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTa                                <  2:379308/149‑1 (MQ=255)
                                                                                                                 cGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTa                                <  1:3497075/149‑1 (MQ=255)
                                                                                                                     aatcaaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCtt                            >  2:1706787/1‑149 (MQ=255)
                                                                                                                      atcaaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCttt                           <  1:1301368/149‑1 (MQ=255)
                                                                                                                       tcaaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTg                          >  1:2600566/1‑149 (MQ=255)
                                                                                                                        caaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGa                         <  2:3048603/149‑1 (MQ=255)
                                                                                                                        caaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGa                         <  2:1056841/149‑1 (MQ=255)
                                                                                                                         aaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGAt                        <  2:678661/149‑1 (MQ=255)
                                                                                                                         aaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGAt                        <  2:2741692/149‑1 (MQ=255)
                                                                                                                         aaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGAt                        <  2:1925721/149‑1 (MQ=255)
                                                                                                                            aGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTcc                     <  1:1963965/149‑1 (MQ=255)
                                                                                                                             gCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTccc                    >  2:2088131/1‑149 (MQ=255)
                                                                                                                             gCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTccc                    <  1:78309/149‑1 (MQ=255)
                                                                                                                             gCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTccc                    >  1:1399304/1‑149 (MQ=255)
                                                                                                                              cATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCt                   <  2:908064/149‑1 (MQ=255)
                                                                                                                                tttCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGa                 <  2:3905862/149‑1 (MQ=255)
                                                                                                                                  tCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGaaa               <  1:1444118/149‑1 (MQ=255)
                                                                                                                                      cAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGt           <  1:4036968/149‑1 (MQ=255)
                                                                                                                                       aaCCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTg          >  1:3500199/1‑149 (MQ=255)
                                                                                                                                        aCCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTgg         >  2:382846/1‑149 (MQ=255)
                                                                                                                                         cccGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGc        <  2:3220918/149‑1 (MQ=255)
                                                                                                                                         cccGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGc        >  2:1792020/1‑149 (MQ=255)
                                                                                                                                         cccGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGc        <  2:1552170/149‑1 (MQ=255)
                                                                                                                                           cGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGcac      >  1:3400462/1‑149 (MQ=255)
                                                                                                                                           cGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGcac      >  1:3400305/1‑149 (MQ=255)
                                                                                                                                            gCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGcaca     >  1:4088546/1‑149 (MQ=255)
                                                                                                                                              tGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGCACAAt   <  2:2371028/149‑1 (MQ=255)
                                                                                                                                               gagaGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGCACAATg  >  2:1325192/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGCACAATG  >  CP009273/2915852‑2916143

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: