Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,363,680 C→G 5.5% G29A (GGC→GCC)  glmM ← phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,363,6800CG 5.5% 238.8 / 3.3 72G29A (GGC→GCC) glmMphosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base G (2/2);  total (36/36)
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

GAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTACCACATCAAC  >  NC_002947/5363541‑5363826
                                                                                                                                           |                                                                                                                                                   
gAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCa                                                                                                                                            <  6:160182/149‑1 (MQ=255)
     gcaTCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCAt                                                                                                                                        >  4:164061/1‑148 (MQ=255)
          acaTCGGCACCTGCCGCGGACAGGCCTGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTc                                                                                                                                  <  2:174184/149‑1 (MQ=255)
          acaTCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGt                                                                                                                                   <  6:246325/148‑1 (MQ=255)
            aTCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGGCAAGCTTCAGCATGAAGTCAg                                                                                                                                >  4:1345/1‑149 (MQ=255)
            aTCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGGCAAGCTTCAGCATGAAGTCAg                                                                                                                                >  4:6612/1‑149 (MQ=255)
            aTCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAg                                                                                                                                >  2:159881/1‑149 (MQ=255)
             tCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAgg                                                                                                                               >  1:337822/1‑149 (MQ=255)
             tCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCATGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAgg                                                                                                                               >  2:258066/1‑149 (MQ=255)
               ggCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGc                                                                                                                              <  3:279945/148‑1 (MQ=255)
                        cgcgGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCggg                                                                                                                    >  6:240503/1‑148 (MQ=255)
                                 gCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCg                                                                                                              <  8:128233/145‑1 (MQ=255)
                                 gCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCg                                                                                                              >  7:128233/1‑145 (MQ=255)
                                  ccGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATg                                                                                                                                       >  2:344089/1‑120 (MQ=255)
                                  ccGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATg                                                                                                                                       <  1:344089/120‑1 (MQ=255)
                                             tGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGaa                                                                                              >  5:8935/1‑149 (MQ=255)
                                               cGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGGCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATg                                                                                            <  3:6612/149‑1 (MQ=255)
                                               cGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGGCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATg                                                                                            <  3:1345/149‑1 (MQ=255)
                                                   cTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGt                                                                                        <  5:194516/149‑1 (MQ=255)
                                                    tCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTc                                                                                       >  3:151533/1‑149 (MQ=255)
                                                          aTGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTAc                                                                                  >  4:114560/1‑148 (MQ=255)
                                                              aGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCgg                                                                                                     >  2:356773/1‑125 (MQ=255)
                                                              aGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCgg                                                                                                     <  1:356773/125‑1 (MQ=255)
                                                               gCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGt                                                                            <  1:159881/149‑1 (MQ=255)
                                                                ccGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGCGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGt                                                                            <  1:258066/149‑1 (MQ=255)
                                                                  ggAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTAtt                                                                         >  5:145448/1‑149 (MQ=255)
                                                                   gagaTACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGAt                                                                                                                         >  3:138120/1‑101 (MQ=255)
                                                                   gagaTACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGAt                                                                                                                         <  4:138120/101‑1 (MQ=255)
                                                                    agagACGCGTGGCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTAtttt                                                                       <  3:276434/149‑1 (MQ=255)
                                                                    agaTACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTAtttt                                                                       <  3:164061/149‑1 (MQ=255)
                                                                    agaTACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTAtttt                                                                       <  4:354843/149‑1 (MQ=255)
                                                                     gaTACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTc                                                                      >  5:373434/1‑149 (MQ=255)
                                                                          gcgTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCt                                                                  >  4:156108/1‑148 (MQ=255)
                                                                             tgtCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATa                                                              <  5:240503/149‑1 (MQ=255)
                                                                               tCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAg                                                             <  3:156108/148‑1 (MQ=255)
                                                                                ccTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAgg                                                            >  4:52916/1‑148 (MQ=255)
                                                                                 cTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGa                                                          <  4:151533/149‑1 (MQ=255)
                                                                                     cccACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCTGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTc                                                       <  5:120219/148‑1 (MQ=255)
                                                                                     cccACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTc                                                       <  6:8935/148‑1 (MQ=255)
                                                                                     cccACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTc                                                       <  3:114560/148‑1 (MQ=255)
                                                                                      ccaccaGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCaaa                                                                              >  8:10393/1‑124 (MQ=255)
                                                                                      ccaccaGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCaaa                                                                              <  7:10393/124‑1 (MQ=255)
                                                                                      ccaccaGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTcc                                                      >  5:166663/1‑148 (MQ=255)
                                                                                       caccaGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCaa                                                    >  6:106408/1‑149 (MQ=255)
                                                                                            gCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGctc                                               <  3:52916/149‑1 (MQ=255)
                                                                                            gCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGctc                                               >  2:348642/1‑149 (MQ=255)
                                                                                             cACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGctct                                              >  3:285842/1‑149 (MQ=255)
                                                                                             cACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGctct                                              >  5:23871/1‑149 (MQ=255)
                                                                                              aCGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAAt                                                                                             <  3:265898/101‑1 (MQ=255)
                                                                                              aCGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAAt                                                                                             >  4:265898/1‑101 (MQ=255)
                                                                                                gcgGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTAt                                           >  4:363091/1‑149 (MQ=255)
                                                                                                gcgGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTAt                                           >  8:188656/1‑149 (MQ=255)
                                                                                                  ggCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTc                                         <  8:355133/149‑1 (MQ=255)
                                                                                                    cAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAgg                                                            <  8:24334/128‑1 (MQ=255)
                                                                                                      aTGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCg                                     <  1:348642/149‑1 (MQ=255)
                                                                                                         gCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGttt                                  <  3:363091/149‑1 (MQ=255)
                                                                                                          ccTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGtttt                                 <  6:145448/149‑1 (MQ=255)
                                                                                                               cttgcGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCAcc                            >  4:193697/1‑149 (MQ=255)
                                                                                                                     gAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGc                                      <  1:56642/133‑1 (MQ=255)
                                                                                                                     gAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGc                                      >  2:56642/1‑133 (MQ=255)
                                                                                                                       aGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGAt                    <  4:285842/149‑1 (MQ=255)
                                                                                                                           cATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATg                >  2:184797/1‑149 (MQ=255)
                                                                                                                                cggccgCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTAc           >  5:22469/1‑149 (MQ=255)
                                                                                                                                 ggccgCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTAcc          <  6:23871/149‑1 (MQ=255)
                                                                                                                                 ggccgCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTAcc          <  2:337822/149‑1 (MQ=255)
                                                                                                                                 ggccgCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTAc           >  2:290404/1‑148 (MQ=255)
                                                                                                                                   ccgccCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTACcac        <  6:166663/149‑1 (MQ=255)
                                                                                                                                     gccCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTACCACAt      >  4:250110/1‑149 (MQ=255)
                                                                                                                                      cccAGCCAAGCTTCAGCATGAAGTAAGGCGTGATCGGGTAGTCGCCGACGCGGCCAAGAATGCCGTTGGTACCAAAGTATTTTCTGATCATAGGGACTCCAAGGCTCTTATTCGGCGTTTTGCACAGCAGCGATCATGCGTCCCACAt      >  5:140765/1‑148 (MQ=255)
                                                                                                                                       ccAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATCCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTACCACATc     >  6:119398/1‑148 (MQ=255)
                                                                                                                                         aGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTACCACATCAAc  <  6:22469/149‑1 (MQ=255)
                                                                                                                                         aGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTACCACATCAAc  >  4:340809/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATGCGTACCACATCAAC  >  NC_002947/5363541‑5363826

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