Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,270,401 C→T 8.4% G339G (GGC→GGT PP_4643 → xanthine/uracil permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,270,4010CT 8.4% 599.2 / 31.9 203G339G (GGC→GGTPP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base C (82/104);  new base T (10/7);  total (92/111)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 204 total aligned reads displayed.

CGCGCCGACATGGCATCGAGCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATC  >  NC_002947/5270252‑5270548
                                                                                                                                                     |                                                                                                                                                   
cgcgcCGACATGGCATCGAGCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGc                                                                                                                                                     <  1:2104868/150‑1 (MQ=255)
   gcCGACATGGCATCGAGCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTtg                                                                                                                                                  >  1:4409283/1‑150 (MQ=255)
     cGACATGGCATCGAGCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTtgct                                                                                                                                                >  2:2248985/1‑150 (MQ=255)
      gACATGGCATCGAGCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTtgctg                                                                                                                                               >  1:1685298/1‑150 (MQ=255)
        cATGGCATCGAGCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGcc                                                                                                                                             <  1:5098526/150‑1 (MQ=255)
              aTCGAGCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCAt                                                                                                                                       <  2:3780650/150‑1 (MQ=255)
                cGAGCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATgg                                                                                                                                     <  2:772964/150‑1 (MQ=255)
                cGAGCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATgg                                                                                                                                     <  1:2187315/150‑1 (MQ=255)
                   gCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGcc                                                                                                                                  >  1:4967010/1‑150 (MQ=255)
                     ctCCTGGCGCCGCTTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCgg                                                                                                                                <  2:4679397/149‑1 (MQ=255)
                     aTCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCgg                                                                                                                                <  1:346046/150‑1 (MQ=255)
                     aTCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCgg                                                                                                                                <  1:4869972/150‑1 (MQ=255)
                      tCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCggg                                                                                                                               <  1:2070640/150‑1 (MQ=255)
                      tCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCggg                                                                                                                               <  1:3790724/150‑1 (MQ=255)
                       ccTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGt                                                                                                                              <  2:1595129/150‑1 (MQ=255)
                        cTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTg                                                                                                                             <  2:1967596/150‑1 (MQ=255)
                           gcgcCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGAtt                                                                                                                          <  2:1411136/150‑1 (MQ=255)
                             gcCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTgc                                                                                                                        >  1:2167485/1‑150 (MQ=255)
                                 aTTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTgccgcc                                                                                                                    >  2:5091377/1‑150 (MQ=255)
                                 aTTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTgccgcc                                                                                                                    >  2:3584655/1‑150 (MQ=255)
                                   tttCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGt                                                                                                                  >  1:1812765/1‑150 (MQ=255)
                                    ttCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTg                                                                                                                 <  2:2540034/150‑1 (MQ=255)
                                     tCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGc                                                                                                                <  1:4286384/150‑1 (MQ=255)
                                     tCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGc                                                                                                                <  1:3123290/150‑1 (MQ=255)
                                      cGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGcc                                                                                                               <  2:361665/150‑1 (MQ=255)
                                             ttcACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCg                                                                                                        >  1:4101509/1‑150 (MQ=255)
                                              tcACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCgg                                                                                                       <  2:13671/150‑1 (MQ=255)
                                              tcACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCgg                                                                                                       <  2:13783/150‑1 (MQ=255)
                                                  ccAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACt                                                                                                   >  1:3607570/1‑150 (MQ=255)
                                                   cAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTg                                                                                                  >  2:1467937/1‑150 (MQ=255)
                                                     gagCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTggg                                                                                                >  2:3838920/1‑150 (MQ=255)
                                                       gcgcGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGggcg                                                                                              >  1:2222816/1‑150 (MQ=255)
                                                           gTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCgg                                                                                                                                <  2:1003477/112‑1 (MQ=255)
                                                           gTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCgg                                                                                                                                >  1:1003477/1‑112 (MQ=255)
                                                           gTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGATGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGgcgc                                                                                          >  2:17810/1‑150 (MQ=255)
                                                             tCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCGACCCCGGTACTGGGCGGCGCCg                                                                                        <  2:2370218/150‑1 (MQ=255)
                                                               gCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCggg                                                                                      >  2:4649398/1‑149 (MQ=255)
                                                                 ccAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTtgctg                                                                                                                                               >  1:1916187/1‑91 (MQ=255)
                                                                  cAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATc                                                                                   <  2:4891844/150‑1 (MQ=255)
                                                                  cAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATc                                                                                   >  2:640874/1‑150 (MQ=255)
                                                                  cAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATc                                                                                   <  1:3068244/150‑1 (MQ=255)
                                                                      aCGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGc                                                                               <  2:2359592/150‑1 (MQ=255)
                                                                       cGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTAATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTCTCGTGCt                                                                              >  2:1172604/1‑150 (MQ=255)
                                                                            gCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCg                                                                         <  2:3630786/150‑1 (MQ=255)
                                                                            gCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCg                                                                         >  1:2849575/1‑150 (MQ=255)
                                                                              cTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGc                                                                       >  2:1309880/1‑150 (MQ=255)
                                                                                    gCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTg                                                                 >  1:1550165/1‑150 (MQ=255)
                                                                                    gCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTg                                                                 <  1:2302776/150‑1 (MQ=255)
                                                                                      cGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTTCGGGCGGCGTGTTTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggc                                                               <  2:3325291/150‑1 (MQ=255)
                                                                                       gTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcg                                                              <  2:1204812/150‑1 (MQ=255)
                                                                                       gTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcg                                                              <  2:1718718/150‑1 (MQ=255)
                                                                                       gTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcg                                                              <  2:1922105/150‑1 (MQ=255)
                                                                                       gTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcg                                                              <  2:46966/150‑1 (MQ=255)
                                                                                        tGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcgg                                                             <  1:4798824/150‑1 (MQ=255)
                                                                                            cGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAg                                                         >  1:4428054/1‑150 (MQ=255)
                                                                                             ggggTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGc                                                        >  2:2177345/1‑150 (MQ=255)
                                                                                                 tCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTa                                                    >  2:3593689/1‑150 (MQ=255)
                                                                                                   aaGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATc                                                  >  1:4361661/1‑150 (MQ=255)
                                                                                                    agagCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATcc                                                 <  1:3344071/150‑1 (MQ=255)
                                                                                                    agagCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATcc                                                 >  1:1977743/1‑150 (MQ=255)
                                                                                                    aaagCCGTTATGTGGTGGCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGACATcc                                                 >  1:3899560/3‑150 (MQ=255)
                                                                                                      agCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGc                                               >  1:5425566/1‑150 (MQ=255)
                                                                                                       gCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCa                                              >  2:3895145/1‑150 (MQ=255)
                                                                                                        ccGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCAc                                             <  1:4434460/150‑1 (MQ=255)
                                                                                                         cGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCAcc                                            <  1:1624562/150‑1 (MQ=255)
                                                                                                           cTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCt                                          <  1:133250/150‑1 (MQ=255)
                                                                                                             atgtggtgGCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACttt                                        <  1:4234486/148‑3 (MQ=255)
                                                                                                               gtggtgGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTcc                                      >  2:2456649/1‑150 (MQ=255)
                                                                                                                tggtggCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCa                                     >  1:1435933/1‑150 (MQ=255)
                                                                                                                  gtggCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAg                                   <  2:3106788/150‑1 (MQ=255)
                                                                                                                   tggCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAgg                                  >  1:3381871/1‑150 (MQ=255)
                                                                                                                    ggCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGt                                 >  1:2038246/1‑150 (MQ=255)
                                                                                                                     gcaaCGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTc                                <  2:3043287/147‑1 (MQ=255)
                                                                                                                       tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                              >  2:4279429/1‑150 (MQ=255)
                                                                                                                       tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                              >  1:1698108/1‑150 (MQ=255)
                                                                                                                       tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                              >  2:3205599/1‑150 (MQ=255)
                                                                                                                       tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                              >  1:926912/1‑150 (MQ=255)
                                                                                                                       tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                              >  1:4754118/1‑150 (MQ=255)
                                                                                                                       tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                              >  2:2538218/1‑150 (MQ=255)
                                                                                                                       aaCGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTcga                              >  2:322020/2‑148 (MQ=255)
                                                                                                                       aaCGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTcga                              >  2:2345711/2‑148 (MQ=255)
                                                                                                                         cGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCt                            <  2:3778870/150‑1 (MQ=255)
                                                                                                                          gggCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTa                           >  1:633028/1‑150 (MQ=255)
                                                                                                                           ggcggcGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTAc                          <  1:1528156/150‑1 (MQ=255)
                                                                                                                            gcggcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACa                         >  2:2521601/1‑150 (MQ=255)
                                                                                                                              ggcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAg                       >  2:29862/1‑150 (MQ=255)
                                                                                                                                 gTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaac                    <  1:1736176/150‑1 (MQ=255)
                                                                                                                                    aTTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaacaac                 <  1:582110/150‑1 (MQ=255)
                                                                                                                                         ggtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCaa            >  2:3316111/1‑150 (MQ=255)
                                                                                                                                         ggtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCaa            <  1:5348506/150‑1 (MQ=255)
                                                                                                                                          gtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc           <  2:4428942/150‑1 (MQ=255)
                                                                                                                                          gtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc           <  1:4571253/150‑1 (MQ=255)
                                                                                                                                          gtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc           <  2:1685298/150‑1 (MQ=255)
                                                                                                                                          gtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc           <  1:2392660/150‑1 (MQ=255)
                                                                                                                                           tggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAcc          <  2:2117735/150‑1 (MQ=255)
                                                                                                                                            cacaCTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCt         <  2:4802976/147‑1 (MQ=255)
                                                                                                                                             gtACTCGGCTTGCTGCCAGTCATGAGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTg        >  1:3059547/1‑150 (MQ=255)
                                                                                                                                              tACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGa       <  2:4409283/150‑1 (MQ=255)
                                                                                                                                                cTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatc     <  2:4586888/150‑1 (MQ=255)
                                                                                                                                                cTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatc     >  2:2918699/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGCGCCGACATGGCATCGAGCATCCTGGCGCCGATTTTCGGTTCGTTCACCCAGAGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATC  >  NC_002947/5270252‑5270548

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: