Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,654,937 T→C 5.2% intergenic (+43/‑215) aceA → / → PP_4118 isocitrate lyase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,654,9370TC 5.2% 284.5 / 2.8 94intergenic (+43/‑215)aceA/PP_4118isocitrate lyase/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (41/48);  new base C (3/2);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

AGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAGGC  >  NC_002947/4654788‑4655084
                                                                                                                                                         |                                                                                                                                                   
agGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGTCGTGTGAAAGCt                                                                                                                                                     >  1:1018645/1‑150 (MQ=255)
 ggAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTc                                                                                                                                                    <  2:1681050/150‑1 (MQ=255)
  gAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTcc                                                                                                                                                   <  1:460458/150‑1 (MQ=255)
    aGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCgg                                                                                                                                                 <  2:1218282/150‑1 (MQ=255)
    aGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCgg                                                                                                                                                 <  2:1358894/150‑1 (MQ=255)
            cTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGCCGGTGACTTCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGTTCGCGTTGGCCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGt                                                                                                                                         <  1:1870094/150‑1 (MQ=255)
             tGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGtt                                                                                                                                        <  1:213903/150‑1 (MQ=255)
             tGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGtt                                                                                                                                        <  1:945716/150‑1 (MQ=255)
              ggCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGttt                                                                                                                                       >  1:1069111/1‑150 (MQ=255)
              ggCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGttt                                                                                                                                       <  2:1888060/150‑1 (MQ=255)
              ggCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGttt                                                                                                                                       <  2:554869/150‑1 (MQ=255)
                cTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GAGCTTTCACAAGCCCCGGGCCTTTGttcaa                                                                                                                                     <  1:364543/150‑4 (MQ=38)
                 tACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGCCCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTa                                                                                                                                    <  1:616340/150‑1 (MQ=255)
                    ttCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTGGACCGCGGAAGAGCAGGTCCACTGACAGTGGGTTGCCGTTGAACAATGTCCTG‑‑GGCCTTGTGAAAGCCCCGCGCCCTGGCCTTTAccc                                                                                                                                 >  2:855241/1‑150 (MQ=255)
                       gacgacATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCgag                                                                                                                              >  1:308453/1‑150 (MQ=255)
                       gacgacATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCgag                                                                                                                              >  1:503456/1‑150 (MQ=255)
                       gacgacATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCgag                                                                                                                              >  2:1928275/1‑150 (MQ=255)
                          gacATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGaaa                                                                                                                           <  1:1907860/150‑1 (MQ=255)
                          gacATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGaaa                                                                                                                           >  1:704428/1‑150 (MQ=255)
                          gacATGACCACCGTGATCCAGGGCGGCGCTGCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTGCCACTGATAGTGAG‑‑CCAGTGGGGGTAGGGTGCGGTGGGCTTGTGAAAGCTACGGGCCTTTGCGTTTACCGGAGaaa                                                                                                                           >  1:2089585/1‑150 (MQ=255)
                            cATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTa                                                                                                                                    <  1:1594017/139‑1 (MQ=255)
                            cATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTa                                                                                                                                    >  2:1594017/1‑139 (MQ=255)
                             aTGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAggg                                                                                                                        >  1:1971252/1‑150 (MQ=255)
                              tGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggc                                                                                                                       <  2:255461/150‑1 (MQ=255)
                                accaccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GAGCTTTCACAAGCCCCGGGCCTTTg                                                                                                                                          <  1:1534362/129‑1 (MQ=37)
                                accaccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GAGCTTTCACAAGCCCCGGGCCTTTg                                                                                                                                          >  2:1534362/1‑129 (MQ=37)
                                  caccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccggg                                                                                                                   >  2:2016736/1‑150 (MQ=255)
                                   accGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgggc                                                                                                                  <  2:1854215/150‑1 (MQ=255)
                                     cGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCct                                                                                                                >  2:1096411/1‑150 (MQ=255)
                                       tGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgc                                                                                                              >  2:1964956/1‑150 (MQ=255)
                                        gctCCCGGGCGGCGCTTACGCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgct                                                                                                             <  2:1018645/148‑1 (MQ=255)
                                         aTCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgctg                                                                                                            <  1:96873/150‑1 (MQ=255)
                                             aGGGCGACGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCAGGCCTTGGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCa                                                                                                        >  1:754384/1‑150 (MQ=255)
                                                gcggcgCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAtt                                                                                                     <  1:1928275/150‑1 (MQ=255)
                                                gcggcgCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGACAGTTGGCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAtt                                                                                                     <  1:1353009/150‑1 (MQ=255)
                                                  ggcgCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAtttt                                                                                                   >  1:1813081/1‑150 (MQ=255)
                                                    cgcTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTa                                                                                                 >  1:1778099/1‑150 (MQ=255)
                                                     gcTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTAc                                                                                                >  1:2213998/1‑150 (MQ=255)
                                                           tcgGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAc                                                                                                                                   >  1:1957131/1‑109 (MQ=255)
                                                           tcgGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAc                                                                                                                                   <  2:1957131/109‑1 (MQ=255)
                                                              gTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTc                                                                                       >  2:767532/1‑150 (MQ=255)
                                                              gTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTGGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTc                                                                                       >  2:1399598/1‑150 (MQ=255)
                                                                  cTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGtt                                                                                                                                        >  1:1125576/1‑97 (MQ=255)
                                                                  cTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGtt                                                                                                                                        <  2:1125576/97‑1 (MQ=255)
                                                                   tGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGTTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTTCTTTCCAATTTTTACCAGCCGTTCACCGt                                                                                  <  1:63653/150‑1 (MQ=255)
                                                                   tGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGt                                                                                  >  2:2163996/1‑150 (MQ=255)
                                                                   tGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGt                                                                                  <  1:872033/150‑1 (MQ=255)
                                                                     cACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGt                                                                                                                                         >  2:1816401/1‑93 (MQ=255)
                                                                     cACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGt                                                                                                                                         <  1:1816401/93‑1 (MQ=255)
                                                                     cACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGcc                                                                                           >  2:1654680/1‑139 (MQ=255)
                                                                     cACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGcc                                                                                           <  1:1654680/139‑1 (MQ=255)
                                                                     cACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACAAGCCGTTCACCGttt                                                                                >  1:647880/1‑150 (MQ=255)
                                                                          aCCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTccgcc                                                                           <  2:994683/150‑1 (MQ=255)
                                                                             ggTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCt                                                                        >  2:333529/1‑150 (MQ=255)
                                                                               ttCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTAc                                                                      >  2:335567/1‑150 (MQ=255)
                                                                                        agaagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTcc                                                                                                                                                   >  2:2095992/1‑64 (MQ=255)
                                                                                        agaagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTcc                                                                                                                                                   <  1:2095992/64‑1 (MQ=255)
                                                                                        agaagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcgc                                                             <  1:147194/150‑1 (MQ=255)
                                                                                         gaagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcgcg                                                            <  2:1488640/150‑1 (MQ=255)
                                                                                          aagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCgg                                                           <  2:1069111/150‑1 (MQ=255)
                                                                                                aGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAAc                                                     <  2:970282/150‑1 (MQ=255)
                                                                                                 gTTCCACTGATAGTGGGGTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTGCCCGAGAACGCGCCGGGCCTGCTGGACAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAAcc                                                    <  1:1748074/150‑1 (MQ=37)
                                                                                                      cACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCt                                                                                                                                             >  2:509924/1‑57 (MQ=255)
                                                                                                      cACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCt                                                                                                                                             <  1:509924/57‑1 (MQ=255)
                                                                                                      cACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCa                                                >  1:2283502/1‑150 (MQ=255)
                                                                                                       aCTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCAc                                               >  1:861924/1‑150 (MQ=255)
                                                                                                            tAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTaca                                          >  2:1511076/1‑150 (MQ=255)
                                                                                                                 ggTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGAt                                     <  1:767532/150‑1 (MQ=255)
                                                                                                                       cGTTTAGCAAAAGCCCG‑‑GGGCTTGGTAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAGGGGCCGGGCCTGCTGGCCAATTCTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                               <  1:1100736/150‑1 (MQ=25)
                                                                                                                       cGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                               <  2:27430/150‑1 (MQ=255)
                                                                                                                       cGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                               <  1:2016736/150‑1 (MQ=255)
                                                                                                                         ttGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTCCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGCGTAAACCCGCACCTACACTGCTCGggggga                             >  2:919144/1‑150 (MQ=255)
                                                                                                                           gAAAAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                           <  2:754384/150‑1 (MQ=255)
                                                                                                                            aaCAAAGGCCCG‑‑GAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGac                          >  1:1806316/1‑150 (MQ=25)
                                                                                                                             aCAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCCATCTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGaca                         <  1:1964956/150‑1 (MQ=255)
                                                                                                                              cAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTAt                                                                 >  2:1432306/1‑109 (MQ=255)
                                                                                                                              cAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTAt                                                                 <  1:1432306/109‑1 (MQ=255)
                                                                                                                                aaGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGTGAAAGGGCCGGGCCTGGTGTACATTTTTTACCAGCCGTTCACGGTTTCCGCCGCTACCTTATTCGCGGGTAGACCCGCACCCACACTGATCGTAAGTAGacacac                      >  2:176869/1‑150 (MQ=255)
                                                                                                                                 aGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGAcacaca                     <  2:2089585/150‑1 (MQ=255)
                                                                                                                                     ccG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCgg                 >  1:1349716/1‑150 (MQ=255)
                                                                                                                                      cG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGc                <  1:1399598/150‑1 (MQ=255)
                                                                                                                                       g‑‑gggCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgct                                                                                                             >  1:1147557/1‑56 (MQ=255)
                                                                                                                                       g‑‑gggCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgct                                                                                                             <  2:1147557/56‑1 (MQ=255)
                                                                                                                                       g‑‑gggCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCg               >  1:2019579/1‑150 (MQ=255)
                                                                                                                                          gggCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGa              <  2:2213998/150‑1 (MQ=255)
                                                                                                                                           ggCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGAt             >  1:2238316/1‑150 (MQ=255)
                                                                                                                                            gCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCATCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGAtt            <  1:639994/150‑1 (MQ=255)
                                                                                                                                            gCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGAtt            <  2:1813081/150‑1 (MQ=255)
                                                                                                                                              ttGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTgg          <  2:861924/150‑1 (MQ=255)
                                                                                                                                              ttGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTgg          >  1:22300/1‑150 (MQ=255)
                                                                                                                                                gtgAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGgc        <  1:2163996/150‑1 (MQ=255)
                                                                                                                                                  gAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGgcgc      >  1:491053/1‑150 (MQ=255)
                                                                                                                                                   aaaGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCa     >  1:735833/1‑150 (MQ=255)
                                                                                                                                                     aaGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGTCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAg    <  1:1096411/150‑1 (MQ=255)
                                                                                                                                                       gCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAGGc  >  2:171668/1‑150 (MQ=255)
                                                                                                                                                         |                                                                                                                                                   
AGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCG‑‑GGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAGGC  >  NC_002947/4654788‑4655084

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: