Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,344,121 A→C 7.9% C8G (TGC→GGC)  PP_1168 ← TRAP dicarboxylate transporter subunit DctQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,344,1210AC 7.9% 206.2 / 10.7 76C8G (TGC→GGC) PP_1168TRAP dicarboxylate transporter subunit DctQ
Reads supporting (aligned to +/- strand):  ref base A (34/36);  new base C (3/3);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

GCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCGAC  >  NC_002947/1343972‑1344265
                                                                                                                                                     |                                                                                                                                                
gCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCa                                                                                                                                                  <  1:1008253/150‑1 (MQ=255)
    gCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGc                                                                                                                                              <  1:1577649/150‑1 (MQ=255)
      gAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCg                                                                                                                                            >  1:1328821/1‑150 (MQ=255)
        aGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAg                                                                                                                                          >  1:1713687/1‑150 (MQ=255)
          aTGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGt                                                                                                                                        <  2:203422/150‑1 (MQ=255)
          aTGCCCCACGAGATGACCAGTGTTATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGt                                                                                                                                        <  1:1880503/150‑1 (MQ=255)
             ccccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGa                                                                                                                                     >  1:1363665/1‑150 (MQ=255)
               ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGc                                                                                                                                              <  1:97560/139‑1 (MQ=255)
               ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGc                                                                                                                                              >  2:97560/1‑139 (MQ=255)
               ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGact                                                                                                                                   >  2:1470569/1‑148 (MQ=255)
               ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAtt                                                                                                                                   >  1:1563518/1‑150 (MQ=255)
                          tCAGCGTCATGGCCAGGTTCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcg                                                                                                                        >  2:993634/1‑150 (MQ=255)
                           cAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcgg                                                                                                                       >  2:893073/1‑150 (MQ=255)
                             gCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATCGGCATCACCCAAGCGTCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCt                                                                                                                     <  2:272446/150‑1 (MQ=255)
                               gTCAGGGCCAGGATGGACAGGCCCGCCATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATACGCATCACCAGAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGt                                                                                                                   <  2:1559617/150‑1 (MQ=255)
                                      ccAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTg                                                                                                              >  1:1341522/1‑148 (MQ=255)
                                      ccAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTg                                                                                                              <  2:1341522/148‑1 (MQ=255)
                                       cAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGAc                                                                                                           >  2:1423699/1‑150 (MQ=255)
                                        aGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACg                                                                                                          >  2:410944/1‑150 (MQ=255)
                                        aGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACg                                                                                                          <  1:631773/150‑1 (MQ=255)
                                            ccGAAAGGCCCGCGATCCAGATGAAACTGCGGTACAGCGTGTCGTTCACGCTCACGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGCGCGCCGCACTGCCGAGGTCGATTGTGCAAGGGCGGCTGTTACTGGACGGCGt                                                                                                      <  2:153087/149‑1 (MQ=255)
                                               aCAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTAc                                                                                                                >  2:456367/1‑137 (MQ=255)
                                               aCAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTAc                                                                                                                <  1:456367/137‑1 (MQ=255)
                                               aCAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGa                                                                                                            >  2:1175433/1‑141 (MQ=255)
                                               aCAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGa                                                                                                            <  1:1175433/141‑1 (MQ=255)
                                               aCAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGa                                                                                                   >  2:438667/1‑150 (MQ=255)
                                               aCAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGa                                                                                                   >  1:384589/1‑150 (MQ=255)
                                                 aGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGAACAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATa                                                                                                 <  1:859296/150‑1 (MQ=255)
                                                   gCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAGGGGCGGCTGTTACTGGACGGCGTCGATa                                                                                                 <  1:296679/148‑1 (MQ=255)
                                                   gCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATa                                                                                                 >  2:296679/1‑148 (MQ=255)
                                                    cccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                              >  1:1932510/1‑150 (MQ=255)
                                                    cccGCGATACAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGACTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                              <  1:1470569/150‑1 (MQ=255)
                                                       gcgATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTtc                                                                                           >  1:1013038/1‑150 (MQ=255)
                                                              aGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTgg                                                                                                             <  1:1630050/125‑1 (MQ=255)
                                                              aGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTgg                                                                                                             >  2:1630050/1‑125 (MQ=255)
                                                                aTGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCg                                                                                  >  2:433500/1‑150 (MQ=255)
                                                                           gTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGTTTGGGCAAGGGCGGCTGTTAc                                                                                                                >  2:1628654/1‑109 (MQ=255)
                                                                           gTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGTTTGGGCAAGGGCGGCTGTTAc                                                                                                                <  1:1628654/109‑1 (MQ=255)
                                                                            tACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTg                                                                                                                                  <  1:1327131/90‑1 (MQ=255)
                                                                            tACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTg                                                                                                                                  >  2:1327131/1‑90 (MQ=255)
                                                                             aCAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgttc                                                                     <  2:384589/150‑4 (MQ=255)
                                                                                  gtgtCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCt                                                                                             >  1:1174090/1‑121 (MQ=255)
                                                                                  gtgtCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCt                                                                                             <  2:1174090/121‑1 (MQ=255)
                                                                                      cGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTAttttt                                                            <  2:1563518/150‑1 (MQ=255)
                                                                                       gTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTc                                                           <  1:893073/150‑1 (MQ=255)
                                                                                        ttCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAtt                                                                                                                                   <  2:1684729/77‑1 (MQ=255)
                                                                                        ttCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAtt                                                                                                                                   >  1:1684729/1‑77 (MQ=255)
                                                                                        ttCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcg                                                          >  1:187980/1‑150 (MQ=255)
                                                                                             gCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTAc                                                     <  2:1932510/150‑1 (MQ=255)
                                                                                                 aGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGaa                                                 <  2:187980/150‑1 (MQ=255)
                                                                                                       aGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTg                                           >  1:836231/1‑150 (MQ=255)
                                                                                                            tttCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgtc                                      <  2:1013038/150‑1 (MQ=255)
                                                                                                               cATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTa                                   <  2:1713687/150‑1 (MQ=255)
                                                                                                                 tAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGaa                                                 <  1:467416/134‑1 (MQ=255)
                                                                                                                 tAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGaa                                                 >  2:467416/1‑134 (MQ=255)
                                                                                                                  aGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGcc                                              <  1:1755071/136‑1 (MQ=255)
                                                                                                                  aGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGcc                                              >  2:1755071/1‑136 (MQ=255)
                                                                                                                   ggCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCATTGGTTTCGCGTACCGAAGCCGTGGCGTCGTAGAAc                               <  1:1423699/150‑1 (MQ=255)
                                                                                                                       tCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcgc                                                         <  1:1818717/120‑1 (MQ=255)
                                                                                                                       tCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcgc                                                         >  2:1818717/1‑120 (MQ=255)
                                                                                                                           cAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttctt                       >  2:1653216/1‑150 (MQ=255)
                                                                                                                           cAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttctt                       <  1:1653216/150‑1 (MQ=255)
                                                                                                                           cAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGACGTACTGCACGCGGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttctt                       >  1:1573032/1‑150 (MQ=255)
                                                                                                                               tcGGCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACg                                                                                                          <  2:1624564/62‑1 (MQ=18)
                                                                                                                               tcGGCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACg                                                                                                          >  1:1624564/2‑63 (MQ=18)
                                                                                                                               tcGGCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttctt                       <  2:604756/145‑1 (MQ=255)
                                                                                                                               tcGGCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttctt                       >  1:604756/2‑146 (MQ=255)
                                                                                                                                cGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCCAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCgc                  <  2:1328821/150‑2 (MQ=255)
                                                                                                                                  gCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACg                <  1:410944/150‑1 (MQ=255)
                                                                                                                                  gCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACg                <  1:993634/150‑1 (MQ=255)
                                                                                                                                    aGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGc                            <  2:377363/136‑1 (MQ=255)
                                                                                                                                    aGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGc                            >  1:377363/1‑136 (MQ=255)
                                                                                                                                     gggCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTa             >  2:269159/1‑150 (MQ=255)
                                                                                                                                      ggcggcGGGCGCCGCACTGCCCAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTGCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTaa            >  2:1289847/1‑150 (MQ=255)
                                                                                                                                             ggcgCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTc     <  2:1363665/150‑1 (MQ=255)
                                                                                                                                               cgccgcACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCGa   <  1:438667/150‑1 (MQ=255)
                                                                                                                                                gccgcACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGAGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGGTTCTTGTTGACGGTAATGAAGTAGac  >  1:207473/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCGAC  >  NC_002947/1343972‑1344265

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