Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,712,122 A→T 6.0% intergenic (+71/‑187) ppc → / → adk phosphoenolpyruvate carboxylase/adenylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,712,1220AT 6.0% 210.3 / 4.9 67intergenic (+71/‑187)ppc/adkphosphoenolpyruvate carboxylase/adenylate kinase
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base T (2/2);  total (31/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 = 6.53e-01

CGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCC  >  NC_002947/1711973‑1712268
                                                                                                                                                     |                                                                                                                                                  
cGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCa                                                                                                                                                    <  2:937666/150‑1 (MQ=255)
  ctctGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAgg                                                                                                                                                  <  1:21844/150‑1 (MQ=255)
        aCAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAAc                                                                                                                                            <  1:540344/150‑1 (MQ=255)
          aGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCg                                                                                                                                          <  1:483480/150‑1 (MQ=255)
            cccGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGcc                                                                                                                                        >  1:956248/1‑150 (MQ=255)
             ccGTTGGAGCAGGCCTTGCTGGTGCCCGTGGCCGGTCAAATTCTCGAGCGTCGCAACACCGGCTAAGGGCAAACCCGTGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCt                                                                                                                                       <  1:917500/150‑1 (MQ=255)
             ccGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGGGTAACACGGGCTAAGGGCCAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCt                                                                                                                                       <  2:751097/150‑1 (MQ=255)
             ccGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCt                                                                                                                                       <  1:567890/150‑1 (MQ=255)
             ccGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCAGGGGTCATGTCTGGGTTTTCCGATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCt                                                                                                                                       <  1:725524/150‑1 (MQ=255)
               gTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGc                                                                                                                                     >  1:935800/1‑150 (MQ=255)
                     gCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTc                                                                                                                               >  1:899447/1‑150 (MQ=255)
                        ggCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGc                                                                                                                            >  1:922550/1‑150 (MQ=255)
                            ttGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCa                                                                                                                        <  1:1682092/150‑1 (MQ=255)
                            ttGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACAAAGGGACCGGACAGATCCGGGGCATCGTGCATGACGCGCGTACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCa                                                                                                                        >  1:1692341/1‑150 (MQ=255)
                               cttgTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCACGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGt                                                                                                                     <  2:1719776/147‑1 (MQ=255)
                                tGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTg                                                                                                                    >  2:1837239/1‑150 (MQ=255)
                                          ggCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCt                                                                                                          <  1:167120/150‑1 (MQ=255)
                                             gggTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCACATTCCt                                                                                                          <  1:785542/147‑1 (MQ=255)
                                             gggTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCACATTCCt                                                                                                          >  2:785542/1‑147 (MQ=255)
                                               gTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAAc                                                                                                     <  1:264770/150‑1 (MQ=255)
                                                   tGCCGCCGGTCGGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTg                                                                                                 <  2:756502/150‑1 (MQ=255)
                                                          gggCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCt                                                                                          >  2:160429/1‑150 (MQ=255)
                                                              tGCGCAACTCCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTc                                                                                      <  1:1366207/150‑1 (MQ=255)
                                                                   aaCACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCGTACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCt                                                                                 >  1:1017414/1‑150 (MQ=255)
                                                                     cacCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCTGGCGCAACCGcc                                                                                                                                        >  1:749145/1‑93 (MQ=255)
                                                                     cacCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCTGGCGCAACCGcc                                                                                                                                        <  2:749145/93‑1 (MQ=255)
                                                                      acCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAgg                                                                                                                                                  <  1:937143/82‑1 (MQ=255)
                                                                      acCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAgg                                                                                                                                                  >  2:937143/1‑82 (MQ=255)
                                                                        cGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCgg                                                                            <  1:1656650/150‑1 (MQ=255)
                                                                          gCTAACGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGgc                                                                          <  1:610659/150‑1 (MQ=255)
                                                                                ggCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCtgtg                                                                    >  1:1608473/1‑150 (MQ=255)
                                                                                ggCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCtgtg                                                                    >  1:8433/1‑150 (MQ=255)
                                                                                     gCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAGCCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCt                                                               >  2:177381/1‑150 (MQ=255)
                                                                                               tGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCt                                                     >  1:1445980/1‑150 (MQ=255)
                                                                                                gTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGc                                                                                           <  1:1094373/111‑1 (MQ=255)
                                                                                                gTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGc                                                                                           >  2:1094373/1‑111 (MQ=255)
                                                                                                   tGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTg                                                 >  2:1015435/1‑150 (MQ=255)
                                                                                                     gTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACgaat                                                                                         <  1:1037693/108‑4 (MQ=255)
                                                                                                     gTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCtt                                                                                         >  2:1037693/1‑108 (MQ=255)
                                                                                                     gTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGAc                                               >  1:1677961/1‑150 (MQ=255)
                                                                                                          ggCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTg                                                 >  1:1128569/1‑143 (MQ=255)
                                                                                                          ggCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTg                                                 <  2:1128569/143‑1 (MQ=255)
                                                                                                          ggCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCGGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGAtt                                          >  1:991091/1‑150 (MQ=255)
                                                                                                          ggCCATACAGGCCCGGCACAAGGAAGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGAtt                                          >  1:154328/1‑150 (MQ=255)
                                                                                                              aTACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAGCCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATg                                      <  1:177381/150‑1 (MQ=255)
                                                                                                              aTACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGAcc                                              <  1:806150/142‑1 (MQ=255)
                                                                                                              aTACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGAcc                                              >  2:806150/1‑142 (MQ=255)
                                                                                                               tACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATgg                                     <  2:956248/150‑1 (MQ=255)
                                                                                                               tACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATgg                                     >  2:1006256/1‑150 (MQ=255)
                                                                                                                  aGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCa                                  <  2:1703175/150‑1 (MQ=255)
                                                                                                                       cGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATAccc                             >  1:1771199/1‑150 (MQ=255)
                                                                                                                         gCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGa                           <  2:8433/150‑1 (MQ=255)
                                                                                                                           acaAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTc                                   >  1:601897/1‑140 (MQ=255)
                                                                                                                           acaAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTc                                   <  2:601897/140‑1 (MQ=255)
                                                                                                                            caAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGAttt                        <  2:1677961/150‑1 (MQ=255)
                                                                                                                              aGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTc                      <  1:1280187/150‑1 (MQ=255)
                                                                                                                              aGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTc                      <  2:1017414/150‑1 (MQ=255)
                                                                                                                               ggACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTcc                     >  2:939315/1‑150 (MQ=255)
                                                                                                                                gACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCg                    <  2:922550/150‑1 (MQ=255)
                                                                                                                                  cGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGa                  <  1:1837239/150‑1 (MQ=255)
                                                                                                                                   ggggAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGAc                 <  1:1015435/150‑1 (MQ=255)
                                                                                                                                      gAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTg              <  2:1692341/150‑1 (MQ=255)
                                                                                                                                            ggTTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTgtgt                                                                   <  2:555178/91‑1 (MQ=255)
                                                                                                                                            ggTTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTgtgt                                                                   >  1:555178/1‑91 (MQ=255)
                                                                                                                                              ttGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGt      >  1:1365552/1‑150 (MQ=255)
                                                                                                                                              ttGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGt      >  2:1315425/1‑150 (MQ=255)
                                                                                                                                                 cgcCAGGCGCAACCGCCTGCCCGGTCAGCTGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGc   <  2:935800/150‑1 (MQ=255)
                                                                                                                                                  gcCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGcc  >  2:1431578/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCC  >  NC_002947/1711973‑1712268

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