Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 458,241 C→G 6.9% G181A (GGC→GCC)  pqqC ← pyrroloquinoline‑quinone synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947458,2410CG 6.9% 299.4 / 9.6 97G181A (GGC→GCC) pqqCpyrroloquinoline‑quinone synthase
Reads supporting (aligned to +/- strand):  ref base C (47/43);  new base G (4/3);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TAGGCCATGCTCATGGCATCGAGCATACTCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCAGCTGGCGCGGCGG  >  NC_002947/458097‑458387
                                                                                                                                                |                                                                                                                                                   
tAGGCCATGCTCATGGCATCGAGCATACTCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAgg                                                                                                                                                <  2:1779451/150‑1 (MQ=255)
    ccATGCTCATGGCATCGAGCATACTCCAGAGGATATCCGGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCTAGGCCTCGGTTTCCACATCGCGTCGCGCCTGGGCCAGGCggg                                                                                                                                            >  1:1331768/1‑150 (MQ=255)
     cATGCTCATGGCATCGAGCATACTCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCCCGCCTGGGCCAGGCgcga                                                                                                                                           <  2:2110508/150‑4 (MQ=255)
      aTGCACATGCCATCGAGAATACTACAGTGGATCTCTAGCCTGAACTGAAGCATTTCCAGTATCCGCTACTCGCCAGCACGGGTGGTGTAGCGCTGCAGGGTAAACGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGGCCAGGCgcgac                                                                                                                                          <  2:1331768/150‑5 (MQ=255)
         cTCATGGCATCGAGCATACTCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACgg                                                                                                                                       <  1:488125/150‑1 (MQ=255)
             tGGCATCGAGCATACTCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAGTCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAAt                                                                                                                                   >  2:1275966/1‑150 (MQ=255)
             tGGCATCGAGCATACTCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAAt                                                                                                                                   >  1:137828/1‑150 (MQ=255)
                 aTCGAGCATACTCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTc                                                                                                                               <  1:1076502/150‑1 (MQ=255)
                    gAGCATACTCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATa                                                                                                                            <  1:688447/150‑1 (MQ=255)
                     aGCATACTCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAg                                                                                                                           <  2:2223245/150‑1 (MQ=255)
                            tcgagagGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCcggcgg                                                                                                                    >  1:2405540/4‑150 (MQ=255)
                            tCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCcggcgg                                                                                                                    >  1:157329/1‑150 (MQ=255)
                              cAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTAc                                                                                                                                         <  1:58085/127‑1 (MQ=255)
                              cAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTAc                                                                                                                                         >  2:58085/1‑127 (MQ=255)
                              cAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGCGTAGTGCTGCAGGTTAATCGCCCGGCCCTCAACCAAATCCGGCCGGAGCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCgggg                                                                                                                  >  1:542756/1‑150 (MQ=255)
                               agagGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGt                                                                                                                 <  2:2060749/150‑1 (MQ=255)
                               agagGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGt                                                                                                                 <  2:157329/150‑1 (MQ=255)
                                   gATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGAGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGAt                                                                                                             >  2:252581/1‑150 (MQ=255)
                                        ccAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACg                                                                                                        >  2:2138604/1‑150 (MQ=255)
                                        ccAGCTCGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACg                                                                                                        >  1:1583165/1‑150 (MQ=255)
                                         cAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACgg                                                                                                       <  2:362525/150‑1 (MQ=255)
                                         cAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACgg                                                                                                       <  2:1674200/150‑1 (MQ=255)
                                          aGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACggg                                                                                                      <  1:1406064/150‑1 (MQ=255)
                                           gCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACACCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGt                                                                                                     <  2:208995/150‑1 (MQ=255)
                                            cTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTa                                                                                                    <  1:1383787/150‑1 (MQ=255)
                                                   tGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTAc                                                                                                                                         >  1:1616785/1‑106 (MQ=255)
                                                   tGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTAc                                                                                                                                         <  2:1616785/106‑1 (MQ=255)
                                                   tGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGtgctgc                                                                                             <  1:2138604/150‑1 (MQ=255)
                                                   tGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCATGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGtgctgc                                                                                             >  1:1889211/1‑150 (MQ=255)
                                                      aGTATCTCCAGCATACGCGCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGc                                                                                          >  2:1678379/1‑150 (MQ=255)
                                                       gTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTa                                                                                                                                          <  2:1940622/101‑1 (MQ=255)
                                                       gTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTa                                                                                                                                          >  1:1940622/1‑101 (MQ=255)
                                                        tATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCGACGGGGAGTTGTGCGGccg                                                                                        >  2:1345958/1‑149 (MQ=255)
                                                            tcCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATAGTCCTAGCCGGCGGGGACGGTCGACGGGTAGTGCTGCGGCCACCTg                                                                                    >  1:644427/1‑150 (MQ=255)
                                                            tcCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTg                                                                                    >  1:650153/1‑150 (MQ=255)
                                                              cAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCgg                                                                                                                                             >  2:138272/1‑91 (MQ=255)
                                                              cAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCgg                                                                                                                                             <  1:138272/91‑1 (MQ=255)
                                                              cAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTc                                                                                  >  1:1503034/1‑150 (MQ=255)
                                                              cAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTc                                                                                  <  1:343974/150‑1 (MQ=255)
                                                              cAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTc                                                                                  <  1:1345958/150‑1 (MQ=255)
                                                                gCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATAGTCATAGCCGGCGGGGTAGCTCCACGGGTCGTGCTGCGGCCACCTGTGCa                                                                                >  2:198332/1‑150 (MQ=255)
                                                                gCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGCGGAGTTGCGGCCAGCTGTCCa                                                                                >  2:794747/1‑150 (MQ=255)
                                                                   tACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATc                                                                                                            <  1:501715/119‑1 (MQ=255)
                                                                   tACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATc                                                                                                            >  2:501715/1‑119 (MQ=255)
                                                                      gCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTc                                                                                  <  2:1379570/142‑1 (MQ=255)
                                                                      gCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTc                                                                                  >  1:1379570/1‑142 (MQ=255)
                                                                        tCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGGCCAGGCgcga                                                                                                                                           >  2:421022/1‑80 (MQ=255)
                                                                        tCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGAc                                                                        >  1:518694/1‑150 (MQ=255)
                                                                        tCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCCCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGGCCAGGCgcga                                                                                                                                           <  1:421022/83‑4 (MQ=255)
                                                                          cTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGCGCTGCGGCCAGCTGTCCAGGCGTGACtg                                                                      <  1:794747/150‑1 (MQ=255)
                                                                           tGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTCCGCAAATACTCATAGCCGGCGGGGTCGATCCACGGGGAGTGCTGCGGCCCGCTGTCCAGGCGTGACggg                                                                     >  2:2257541/1‑150 (MQ=255)
                                                                                cGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGtgctgc                                                                                             >  1:1428301/1‑121 (MQ=255)
                                                                                cGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGtgctgc                                                                                             <  2:1428301/121‑1 (MQ=255)
                                                                                    cGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGc                                                           <  1:198332/150‑1 (MQ=255)
                                                                                     gggTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTa                                                                                                    >  1:1324902/1‑109 (MQ=255)
                                                                                     gggTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTa                                                                                                    <  2:1324902/109‑1 (MQ=255)
                                                                                         ggtgTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGc                                                      <  2:650153/150‑1 (MQ=255)
                                                                                          gtgtAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCCTa                                                                                                                            <  2:2211634/80‑1 (MQ=255)
                                                                                          gtgtAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATa                                                                                                                            >  1:2211634/1‑80 (MQ=255)
                                                                                             tAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGGAGATCCACGGGTAGTGCTGCGGTCAGCGGTCCAGGCGTGCATGGTGGATTTGCGGGGCGaaa                                                  >  1:1797979/1‑149 (MQ=255)
                                                                                               gTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAg                                                <  1:252581/150‑1 (MQ=255)
                                                                                                   tgcAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCcg                                                                                                                        <  2:1747037/75‑1 (MQ=255)
                                                                                                   tgcAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCcg                                                                                                                        >  1:1747037/1‑75 (MQ=255)
                                                                                                   tgcAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCg                                            >  2:534456/1‑150 (MQ=255)
                                                                                                   tgcAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCg                                            >  1:2043325/1‑150 (MQ=255)
                                                                                                    gcAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCgg                                           <  2:1503034/150‑1 (MQ=255)
                                                                                                    gcAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCgg                                           >  1:619453/1‑150 (MQ=255)
                                                                                                      aGGGGAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACAGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTc                                         <  1:2257541/150‑1 (MQ=255)
                                                                                                       gggTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCa                                        <  2:2405540/150‑1 (MQ=255)
                                                                                                         gTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACCGCTCGGTCAGc                                      >  1:646728/1‑150 (MQ=255)
                                                                                                             tCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGACGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAAc                                  >  1:2290889/1‑150 (MQ=255)
                                                                                                               gCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCgggg                                                       >  2:1755305/1‑127 (MQ=255)
                                                                                                               gCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCgggg                                                       <  1:1755305/127‑1 (MQ=255)
                                                                                                                ccAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCgg                                                         <  2:207692/124‑1 (MQ=255)
                                                                                                                ccAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCgg                                                         >  1:207692/1‑124 (MQ=255)
                                                                                                                 cAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAActgct                              >  2:2044528/1‑150 (MQ=255)
                                                                                                                    gCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACtgg                                                                     <  1:1800984/109‑1 (MQ=255)
                                                                                                                    gCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACtgg                                                                     >  2:1800984/1‑109 (MQ=255)
                                                                                                                     ccGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTggcg                          >  1:973791/1‑150 (MQ=255)
                                                                                                                          ttCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTc                     >  1:1217562/1‑150 (MQ=255)
                                                                                                                            ccACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATGTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCt                   <  1:460449/150‑1 (MQ=255)
                                                                                                                              acaTCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGc                 >  2:1014289/1‑150 (MQ=255)
                                                                                                                              acaTCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGATGAGGCCAGCGGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGc                 >  2:2374130/1‑150 (MQ=255)
                                                                                                                               caTCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCg                                     >  2:1415519/1‑129 (MQ=255)
                                                                                                                               caTCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCg                                     <  1:1415519/129‑1 (MQ=255)
                                                                                                                                aTCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCa               <  1:1688167/150‑1 (MQ=255)
                                                                                                                                aTCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCa               <  1:631293/150‑1 (MQ=255)
                                                                                                                                   gcgtcgcgCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCAGCt            >  2:1993653/1‑150 (MQ=255)
                                                                                                                                     tacccgcCTGGGCCAGGCGGGTACGGGAACAAGCATATCCGATGGCTGCGACCCACGGGCACTGCTGCGGCCCCAACCCCAGCCTGCCTCCCGCGATTGCCGGGGCGAACAGCTCGGTCAGCGAACTGCTGACGGCTTCCCGCCAGCTgg          >  1:550628/5‑150 (MQ=255)
                                                                                                                                     gtcgcgCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCAGCTgg          >  1:554957/1‑150 (MQ=255)
                                                                                                                                        gcgcCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCAGCTGgcgc       <  2:1583165/150‑1 (MQ=255)
                                                                                                                                        ccgcCTGGGCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCAGCTGgcgc       >  2:1845449/2‑150 (MQ=255)
                                                                                                                                          gcCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCAGCTGGCgcgg     >  2:1972322/1‑150 (MQ=255)
                                                                                                                                           ccTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCAGCTGGCgcggc    <  1:1993653/150‑1 (MQ=255)
                                                                                                                                             tGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGtg                                                                                                 <  1:770771/56‑1 (MQ=255)
                                                                                                                                             tGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGtg                                                                                                 >  2:770771/1‑56 (MQ=255)
                                                                                                                                             tGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCAGCTGGCGcggcgg  <  2:1217562/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TAGGCCATGCTCATGGCATCGAGCATACTCCAGAGGATATCCAGCTTGAACTGCAGTATCTCCAGCATACGCTCCTGGCCCGCACGGGTGGTGTAGTGCTGCAGGGTAATCGCCAGGCCGTGTTCCACATCGCGTCGCGCCTGGCCCAGGCGGGTACGGAAATACTCATAGCCGGCGGGGTCGATCCACGGGTAGTGCTGCGGCCAGCTGTCCAGGCGTGACTGGTGGATTTGCGGGGCGAACAGCTCGGTCAGCGAACTGCTGGCGGCTTCCTGCCAGCTGGCGCGGCGG  >  NC_002947/458097‑458387

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: