Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,033,430 G→T 17.3% D73Y (GAC→TAC)  fldB → flavodoxin 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,033,4300GT17.3% 896.8 / 138.2 288D73Y (GAC→TAC) fldBflavodoxin 2
Reads supporting (aligned to +/- strand):  ref base G (101/137);  new base T (28/22);  total (129/159)
Fisher's exact test for biased strand distribution p-value = 8.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

TATTATCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACGGA  >  CP009273/3033282‑3033578
                                                                                                                                                    |                                                                                                                                                    
tattatCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTcg                                                                                                                                                      >  2:4488067/1‑149 (MQ=255)
 attatCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTcga                                                                                                                                                     >  1:3295052/1‑149 (MQ=255)
     tCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCgacgac                                                                                                                                                 >  2:1101568/1‑149 (MQ=255)
     tCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCgacgac                                                                                                                                                 >  2:3399778/1‑149 (MQ=255)
      cGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGAcc                                                                                                                                                <  1:2732510/149‑1 (MQ=255)
      cGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGAcc                                                                                                                                                <  1:2309694/149‑1 (MQ=255)
       ggCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCt                                                                                                                                               <  2:3867344/149‑1 (MQ=255)
            aGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAAcc                                                                                                                                          <  1:3911338/149‑1 (MQ=255)
               aCTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTg                                                                                                                                       >  2:4387118/1‑149 (MQ=255)
                cTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGa                                                                                                                                      <  1:5057658/149‑1 (MQ=255)
                  ggTGACCTTACATAACCTCAAGGACGACTCACCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGTTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAg                                                                                                                                    <  2:2837235/149‑1 (MQ=255)
                   gTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAgg                                                                                                                                   <  2:5123594/149‑1 (MQ=255)
                           aCATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAAttg                                                                                                                           >  2:2599774/1‑149 (MQ=255)
                             aTAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTCAAAttgtt                                                                                                                         >  1:4224244/1‑149 (MQ=255)
                             aTAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAAttgtt                                                                                                                         >  1:1224978/1‑149 (MQ=255)
                              tAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAAttgttg                                                                                                                        <  2:4642263/149‑1 (MQ=255)
                               aaCCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTgc                                                                                                                       <  2:3986734/149‑1 (MQ=255)
                               aaCCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTgc                                                                                                                       >  1:3545143/1‑149 (MQ=255)
                                aCCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTgcg                                                                                                                      >  1:943210/1‑149 (MQ=255)
                                aCCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTgcg                                                                                                                      >  1:1417474/1‑149 (MQ=255)
                                    cacggACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGt                                                                                                                  <  1:2952095/146‑1 (MQ=255)
                                     aaGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTa                                                                                                                 <  2:3253906/149‑1 (MQ=255)
                                     aaGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTa                                                                                                                 <  2:3305074/149‑1 (MQ=255)
                                     aaGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTa                                                                                                                 >  2:5069363/1‑149 (MQ=255)
                                     aaGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGACCTTGAAGGTAAAATTGTTGCGCTGTAt                                                                                                                >  1:3034874/1‑149 (MQ=255)
                                       ggACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATg                                                                                                               <  2:1376905/149‑1 (MQ=255)
                                        gacgacTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATgg                                                                                                              <  1:1927446/149‑1 (MQ=255)
                                           gacTCCCCGAAATTCATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCt                                                                                                           <  1:1170717/149‑1 (MQ=255)
                                           gacTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCt                                                                                                           >  2:4594635/1‑149 (MQ=255)
                                           gacTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCt                                                                                                           <  2:2359261/149‑1 (MQ=255)
                                              tCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTgg                                                                                                        <  2:4277230/149‑1 (MQ=255)
                                               ccccGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGc                                                                                                       >  2:4418351/1‑149 (MQ=255)
                                                   gAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATc                                                                                                   <  1:4694842/149‑1 (MQ=255)
                                                     aaTTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCaa                                                                                                 >  1:3965310/1‑149 (MQ=255)
                                                     aaTTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCaa                                                                                                 <  2:4829223/149‑1 (MQ=255)
                                                     aaTTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCaa                                                                                                 >  1:1369789/1‑149 (MQ=255)
                                                     aaTTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCaa                                                                                                 >  1:5042173/1‑149 (MQ=255)
                                                     aaTTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCaa                                                                                                 >  1:2607992/1‑149 (MQ=255)
                                                     aaTTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCaa                                                                                                 >  2:3762060/1‑149 (MQ=255)
                                                      aTTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAAc                                                                                                <  1:4798780/149‑1 (MQ=255)
                                                        tAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTg                                                                                              >  1:3580690/1‑149 (MQ=255)
                                                        tAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTg                                                                                              >  1:4168788/1‑149 (MQ=255)
                                                          aTGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTggg                                                                                            >  2:498934/1‑149 (MQ=255)
                                                          aTGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTggg                                                                                            >  1:3007544/1‑149 (MQ=255)
                                                          aTGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTggg                                                                                            >  2:2988592/1‑149 (MQ=255)
                                                            ggATCAGTACGATGTGCTCATTCTGGGTATACCGACCTGGGATTTTGGTGAAATCCAGGAAGAATGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGAt                                                                                          <  2:205323/149‑1 (MQ=255)
                                                            ggAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGAt                                                                                          <  2:3672295/149‑1 (MQ=255)
                                                            ggAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGAt                                                                                          <  2:1075253/149‑1 (MQ=255)
                                                            ggAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAAGCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGAt                                                                                          <  2:2294134/149‑1 (MQ=255)
                                                              aGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATAc                                                                                        <  1:1978264/149‑1 (MQ=255)
                                                                cAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACgg                                                                                      <  2:3547541/149‑1 (MQ=255)
                                                                 aGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGc                                                                                     <  2:2915291/149‑1 (MQ=255)
                                                                   tACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGa                                                                                   <  2:4115049/149‑1 (MQ=255)
                                                                   tACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGa                                                                                   <  2:1231423/149‑1 (MQ=255)
                                                                   tACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGa                                                                                   >  1:4895117/1‑149 (MQ=255)
                                                                        tgtgCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGt                                                                              >  2:2749869/1‑149 (MQ=255)
                                                                        tgtgCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGt                                                                              >  1:2915872/1‑149 (MQ=255)
                                                                        tgtgCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGt                                                                              >  1:2786303/1‑149 (MQ=255)
                                                                        tgtgCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGt                                                                              <  1:1758415/149‑1 (MQ=255)
                                                                          tgCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTc                                                                            <  1:4677727/149‑1 (MQ=255)
                                                                             tCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTc                                                                         <  1:4509426/149‑1 (MQ=255)
                                                                             tCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTc                                                                         <  1:172590/149‑1 (MQ=255)
                                                                               aTTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCGGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGa                                                                       <  1:4412873/149‑1 (MQ=255)
                                                                                   tGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATgcg                                                                   >  1:1285245/1‑149 (MQ=255)
                                                                                    gggTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATgcgc                                                                  >  2:1404843/1‑149 (MQ=255)
                                                                                    gggTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATgcgc                                                                  >  1:460405/1‑149 (MQ=255)
                                                                                       tATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCg                                                               >  1:808857/1‑149 (MQ=255)
                                                                                         tCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGt                                                             >  1:1073009/1‑149 (MQ=255)
                                                                                          cccGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTa                                                            <  1:2655230/149‑1 (MQ=255)
                                                                                          cccGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTa                                                            <  1:1457927/149‑1 (MQ=255)
                                                                                                cTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTAtgctgc                                                      >  1:4624384/1‑149 (MQ=255)
                                                                                                 tGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCa                                                     <  1:5038930/149‑1 (MQ=255)
                                                                                                 tGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCa                                                     <  1:2438542/149‑1 (MQ=255)
                                                                                                  gggATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCAt                                                    <  2:1834728/149‑1 (MQ=255)
                                                                                                  gggATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCAt                                                    >  1:2182478/1‑149 (MQ=255)
                                                                                                    gATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGa                                                  <  2:1895207/149‑1 (MQ=255)
                                                                                                    gATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGa                                                  <  1:4284468/149‑1 (MQ=255)
                                                                                                                 tCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGAcc                                     <  2:2194906/149‑1 (MQ=255)
                                                                                                                    aGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCaaa                                  <  1:160244/149‑1 (MQ=255)
                                                                                                                        aGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCg                              >  2:2102357/1‑149 (MQ=255)
                                                                                                                        aGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCg                              >  2:3963361/1‑149 (MQ=255)
                                                                                                                         gACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGt                             >  2:644574/1‑149 (MQ=255)
                                                                                                                         gACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGt                             >  2:1618713/1‑149 (MQ=255)
                                                                                                                         gACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGt                             <  1:1032931/149‑1 (MQ=255)
                                                                                                                         gACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGt                             >  1:4535629/1‑149 (MQ=255)
                                                                                                                             gggAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAg                         <  1:2920629/149‑1 (MQ=255)
                                                                                                                             gggAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAg                         >  1:4071847/1‑149 (MQ=255)
                                                                                                                                aaGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtc                      <  2:4184602/149‑1 (MQ=255)
                                                                                                                                   ccGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtcgtc                   <  1:3138180/149‑1 (MQ=255)
                                                                                                                                      tCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGc                >  2:855238/1‑149 (MQ=255)
                                                                                                                                       cTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCt               <  1:4123869/149‑1 (MQ=255)
                                                                                                                                         gggATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTAc             <  2:314192/149‑1 (MQ=255)
                                                                                                                                           gATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTg           >  1:1865136/1‑149 (MQ=255)
                                                                                                                                              cAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGcc        <  1:37820/149‑1 (MQ=255)
                                                                                                                                              cAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGcc        <  2:2805768/149‑1 (MQ=255)
                                                                                                                                                gCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCaa      >  1:3717829/1‑149 (MQ=255)
                                                                                                                                                gCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCaa      >  1:2085180/1‑149 (MQ=255)
                                                                                                                                                gCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCaa      <  2:2461528/149‑1 (MQ=255)
                                                                                                                                                gCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCGGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCaa      <  1:4387118/149‑1 (MQ=255)
                                                                                                                                                  tCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACg    <  2:2553952/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATTATCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACGGA  >  CP009273/3033282‑3033578

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: