Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,455,122 C→G 7.4% A68G (GCC→GGC)  nicA → nicotinate dehydrogenase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,455,1220CG 7.4% 762.3 / 39.1 229A68G (GCC→GGC) nicAnicotinate dehydrogenase subunit A
Reads supporting (aligned to +/- strand):  ref base C (93/119);  new base G (6/11);  total (99/130)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 244 total aligned reads displayed.

CATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCG  >  NC_002947/4454973‑4455269
                                                                                                                                                     |                                                                                                                                                   
   gCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGcgccgc                                                                                                                                                  <  2:3652178/150‑1 (MQ=255)
     cTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGccgccg                                                                                                                                                <  1:5500669/150‑1 (MQ=255)
       gACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccggc                                                                                                                                              <  1:3941617/150‑1 (MQ=255)
        acacCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccggcc                                                                                                                                             <  2:1096105/150‑1 (MQ=255)
         cacCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccggccg                                                                                                                                            <  2:3898747/150‑1 (MQ=255)
          acCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTGGGCCGGCGGCgccggccaa                                                                                                                                           <  2:408249/150‑3 (MQ=255)
            cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaagg                                                                                                                                         <  2:4509851/150‑5 (MQ=255)
            cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCaa                                                                                                                                         >  1:2597973/1‑150 (MQ=255)
            cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCaa                                                                                                                                         <  2:2697257/150‑1 (MQ=255)
            cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCaa                                                                                                                                         <  1:1964562/150‑1 (MQ=255)
            cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTCGGCCGGCGGCgccggccaagg                                                                                                                                         <  2:1404910/150‑5 (MQ=255)
             ccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaaggg                                                                                                                                        <  2:739173/150‑6 (MQ=255)
               gctgctgctGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAatat                                                                                                                                      <  2:5524831/150‑1 (MQ=255)
                 tgctgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATca                                                                                                                                    <  1:2680711/150‑1 (MQ=255)
                    tgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcacca                                                                                                                                 <  2:2531193/150‑1 (MQ=255)
                            cTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGaa                                                                                                                         >  1:5573151/1‑150 (MQ=255)
                             tGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAg                                                                                                                        >  1:4505796/1‑150 (MQ=255)
                              gcgcAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAgg                                                                                                                       >  1:3338002/1‑150 (MQ=255)
                              gcgcAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAACATCACCACCCTCGAAgg                                                                                                                       >  2:3726358/1‑150 (MQ=255)
                               cgcAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAggg                                                                                                                      <  1:4474717/150‑1 (MQ=255)
                                gcAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGc                                                                                                                     <  1:5407694/150‑1 (MQ=255)
                                 cAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCt                                                                                                                    >  2:2368123/1‑150 (MQ=255)
                                  aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTg                                                                                                                   <  2:3583985/150‑1 (MQ=255)
                                   aCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTgg                                                                                                                  <  2:4853414/150‑1 (MQ=255)
                                      aCCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                               <  1:1423208/150‑1 (MQ=255)
                                        cTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcagc                                                                                                             <  2:1873244/150‑1 (MQ=255)
                                        cTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcagc                                                                                                             <  2:299458/150‑1 (MQ=255)
                                           tgCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCCGCAAg                                                                                                          >  1:1055667/1‑150 (MQ=255)
                                             ccTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGATGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGc                                                                                                        >  1:4090504/1‑150 (MQ=255)
                                             ccTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggc                                                                                                                                              >  1:4742757/1‑112 (MQ=255)
                                              cTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGcc                                                                                                       >  2:3098223/1‑150 (MQ=255)
                                               tGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccg                                                                                                      >  1:4276093/1‑150 (MQ=255)
                                                gAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGATGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgc                                                                                                     >  2:2960450/1‑150 (MQ=255)
                                                  aCGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgcac                                                                                                   <  2:5256824/150‑1 (MQ=255)
                                                      ccccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgcaccgca                                                                                               >  2:2000863/1‑150 (MQ=255)
                                                      ccccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgcaccgca                                                                                               >  1:2867462/1‑150 (MQ=255)
                                                        ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCgcaccgcacc                                                                                             >  2:2391152/1‑150 (MQ=255)
                                                         cAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCAccc                                                                                            <  2:2297586/150‑1 (MQ=255)
                                                         cAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCAccc                                                                                            >  2:1511439/1‑150 (MQ=255)
                                                            gTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGt                                                                                         <  1:514105/150‑1 (MQ=255)
                                                             tACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTg                                                                                        <  1:4344672/150‑1 (MQ=255)
                                                              aCGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGc                                                                                       <  2:4448157/150‑1 (MQ=255)
                                                                  cTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACa                                                                                   <  1:1775099/150‑1 (MQ=255)
                                                                    gCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAgg                                                                                 <  2:623950/150‑1 (MQ=255)
                                                                     cGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGc                                                                                <  2:1597373/150‑1 (MQ=255)
                                                                        ccTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGtt                                                                             <  2:3654120/150‑1 (MQ=255)
                                                                         cTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTc                                                                            <  1:1666311/150‑1 (MQ=255)
                                                                         cTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTc                                                                            <  2:1449989/150‑1 (MQ=255)
                                                                           gggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCAt                                                                          <  1:5290436/150‑1 (MQ=255)
                                                                            ggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATc                                                                         <  1:2651340/150‑1 (MQ=255)
                                                                            ggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATc                                                                         <  1:1371049/150‑1 (MQ=255)
                                                                             gCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcg                                                                        >  2:2022655/1‑150 (MQ=255)
                                                                               gAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgac                                                                      <  2:1060289/150‑1 (MQ=255)
                                                                                 aTGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgacga                                                                    <  1:5614199/150‑1 (MQ=255)
                                                                                     ggggCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAg                                                                >  1:2469786/1‑150 (MQ=255)
                                                                                          gtgtACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAggcggc                                                           <  1:2265404/150‑1 (MQ=255)
                                                                                           tgtACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGcc                                                          <  1:2893358/150‑1 (MQ=255)
                                                                                             tACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCa                                                        <  1:2362164/150‑1 (MQ=255)
                                                                                             tACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCa                                                        >  1:4361279/1‑150 (MQ=255)
                                                                                               cGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAgt                                                      <  2:533410/150‑1 (MQ=255)
                                                                                               cGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAgt                                                      >  1:3713128/1‑150 (MQ=255)
                                                                                                ggTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCATCAGGCGTTCATc                                                                         <  2:4023072/130‑1 (MQ=255)
                                                                                                ggTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATc                                                                         >  1:4023083/1‑130 (MQ=255)
                                                                                                       atcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTAc                                              >  1:1302272/1‑150 (MQ=255)
                                                                                                           aCGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCa                                          <  1:3800502/150‑1 (MQ=255)
                                                                                                            cGGCGTGGCCGCCCGCTCCTGCGTCCTCCCCGTGGCCGGCGCCGCCGGCCGCAAGATCCCCCCCCGCGACGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCAt                                         <  2:4621779/150‑1 (MQ=255)
                                                                                                            cGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCAt                                         <  1:676496/150‑1 (MQ=255)
                                                                                                               cGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGaa                                      <  2:810073/150‑1 (MQ=255)
                                                                                                                gTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAAc                                     >  1:4940014/1‑150 (MQ=255)
                                                                                                                   gccgccCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGc                                  <  2:4598129/150‑1 (MQ=255)
                                                                                                                    ccgcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCa                                 >  1:2082342/1‑150 (MQ=255)
                                                                                                                    ccgcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCa                                 <  1:3265874/150‑1 (MQ=255)
                                                                                                                     cgcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCAt                                >  2:2727998/1‑150 (MQ=255)
                                                                                                                      gcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATg                               >  2:3611078/1‑150 (MQ=255)
                                                                                                                      gcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATg                               >  1:2239748/1‑150 (MQ=255)
                                                                                                                        ccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGACCGGCATGAt                             >  1:3032358/1‑150 (MQ=255)
                                                                                                                           cTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCAt                          <  1:4997710/150‑1 (MQ=255)
                                                                                                                             ccTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGa                        <  2:3050772/150‑1 (MQ=255)
                                                                                                                             ccTGCGTCATCCCCATGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGa                        <  2:4320102/150‑1 (MQ=255)
                                                                                                                              cTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGAc                       >  1:3943735/1‑150 (MQ=255)
                                                                                                                               tGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAg                                                                                  >  1:76868/1‑90 (MQ=255)
                                                                                                                                   tCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCa                  <  2:4683545/150‑1 (MQ=255)
                                                                                                                                      tCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAgg               <  1:1896242/150‑1 (MQ=255)
                                                                                                                                         ttgcggCCGGCGGCGCCGGCCACAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATc                                                                         <  1:2022497/85‑1 (MQ=255)
                                                                                                                                         ccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTCCTGCATGAACGGCATGATCATGACCGCCAAGGccc            >  2:3323832/1‑150 (MQ=255)
                                                                                                                                          cTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCt           >  2:3762945/1‑150 (MQ=255)
                                                                                                                                          cTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCt           <  1:2783287/150‑1 (MQ=255)
                                                                                                                                          cTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCt           >  2:2435380/1‑150 (MQ=255)
                                                                                                                                          cTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTGACTGCATGAACGGCATGATCATGACCGCCAAGGCCCt           >  1:4644240/1‑150 (MQ=255)
                                                                                                                                            tGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGc         <  2:2120716/150‑1 (MQ=255)
                                                                                                                                            tGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGc         >  1:2441159/1‑150 (MQ=255)
                                                                                                                                            cggCCGGCGGCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGc         >  2:3500768/2‑150 (MQ=255)
                                                                                                                                            cggCCGGCGGCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGc         >  2:5089316/2‑150 (MQ=255)
                                                                                                                                              gCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTc       >  2:1308752/1‑150 (MQ=255)
                                                                                                                                               ccGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCg      >  1:3291885/1‑150 (MQ=255)
                                                                                                                                                cGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGa     <  1:2816746/150‑1 (MQ=255)
                                                                                                                                                cGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGa     >  1:3073812/1‑150 (MQ=255)
                                                                                                                                                  gcgcCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATc   <  2:1801899/150‑1 (MQ=255)
                                                                                                                                                  gcgcCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATc   <  1:1457094/150‑1 (MQ=255)
                                                                                                                                                   cgccgcCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCg  >  2:2614731/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCG  >  NC_002947/4454973‑4455269

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: