Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,595,828 T→A 9.3% intergenic (+166/+133) mccA → / ← PP_4068 methylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha/Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,595,8280TA 9.3% 478.0 / 4.2 187intergenic (+166/+133)mccA/PP_4068methylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha/Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (80/78);  new base A (17/5);  total (102/85)
Fisher's exact test for biased strand distribution p-value = 2.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

TTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCA  >  NC_002947/4595680‑4595977
                                                                                                                                                    |                                                                                                                                                     
tttGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATc                                                                                                                                                      <  2:1440562/150‑1 (MQ=255)
attGCGGCCCTAACCCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATc                                                                                                                                                      <  2:2897563/149‑1 (MQ=255)
 tggCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATcc                                                                                                                                                     <  1:500813/148‑1 (MQ=255)
   gCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGa                                                                                                                                                   >  2:2337368/1‑150 (MQ=255)
         cTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCGCCGACATACGAGCgggg                                                                                                                                             >  1:172385/1‑150 (MQ=255)
         cTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGGAAAACCCCGCCGACAACCGATCgggg                                                                                                                                             >  1:1036741/1‑150 (MQ=21)
         cTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATGCGCGGGTAAACCCGCGCCGACAGCCCATCgggg                                                                                                                                             >  2:304282/1‑150 (MQ=21)
         cTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACAGTATTCGCGGGTAAACACGCTCCAACATCCGATCgggg                                                                                                                                             >  1:3017804/1‑150 (MQ=255)
          tAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCGACAACCGATCggggg                                                                                                                                            >  1:3127323/1‑150 (MQ=255)
           aaCGCAACACACGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAACGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGt                                                                                                                                           <  1:1423246/150‑1 (MQ=255)
            aCGCAACACAAGGCCGCCCCTGCAACCGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTg                                                                                                                                          <  1:2410783/150‑1 (MQ=255)
            aCGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTg                                                                                                                                          <  2:3386272/150‑1 (MQ=255)
              gCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGa                                                                                                                                        >  2:5304893/1‑150 (MQ=255)
                aaCACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGaca                                                                                                                                      <  2:4628431/150‑1 (MQ=255)
                     aaGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGt                                                                                                                                 >  1:773741/1‑150 (MQ=255)
                      aGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGtt                                                                                                                                <  2:2145287/150‑1 (MQ=255)
                       ggCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTg                                                                                                                               <  1:2945966/150‑1 (MQ=255)
                         ccgccCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGt                                                                                                                             <  2:2015320/150‑1 (MQ=255)
                         ccgccCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGt                                                                                                                             <  1:1644614/150‑1 (MQ=255)
                           gccCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGa                                                                                                                           <  1:1224211/150‑1 (MQ=255)
                            ccccTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGAt                                                                                                                          <  1:2440765/150‑1 (MQ=255)
                            ccccTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGAt                                                                                                                          <  2:4321862/150‑1 (MQ=255)
                               cTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAg                                                                                                                       <  2:254683/150‑1 (MQ=255)
                               cTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAACCGATCGGGGGTGGACACATGTTGGTGATTAg                                                                                                                       >  2:1058409/1‑150 (MQ=255)
                               cTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCCCGCCGACATCCGATCGGGGGTGGACACATGTTGGTGATTAg                                                                                                                       >  2:1058384/1‑150 (MQ=255)
                                tGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCAACAACCCAGCGGGGGTGGACACAAGTTGGTGATTAgg                                                                                                                      >  2:129481/1‑150 (MQ=18)
                                   aaCAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAAACTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAACCGCACCTACAGCCGAGCGGGGGTGGACACAAGTTGGTGATTAGgcgc                                                                                                                   >  2:1670560/1‑150 (MQ=18)
                                      aGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGGTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAgg                                                                                                                <  1:1018790/150‑1 (MQ=255)
                                        cccGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGcc                                                                                                              <  1:3733059/150‑1 (MQ=255)
                                        cccGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAACCGATCGGGGGTGGACACAAGTTGGTGATTAGGCGCCGGcc                                                                                                              >  1:4044112/1‑150 (MQ=255)
                                         ccGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGAGTAGGCGCAGGCcg                                                                                                             >  2:4905868/1‑150 (MQ=255)
                                          cGTAAGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCcgc                                                                                                            <  1:4819492/150‑1 (MQ=255)
                                            tAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGt                                                                                                          >  2:773741‑M2/1‑149 (MQ=255)
                                              ggAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgt                                                                                                        >  1:351702‑M2/1‑147 (MQ=255)
                                               gAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgta                                                                                                       <  2:1532312‑M2/150‑5 (MQ=255)
                                               gAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgta                                                                                                       >  2:4032805‑M2/1‑146 (MQ=255)
                                                aGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgtag                                                                                                      >  2:552577‑M2/1‑145 (MQ=255)
                                                 gCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAACCGCGCCGACATCCCATCGGGGGTGGTCACAAGTTGGTGATTAGGCGCAGGCCCCGtgtagg                                                                                                     <  1:3904969‑M2/150‑7 (MQ=255)
                                                   gggTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGt                                                                                                                             <  2:4118958/124‑1 (MQ=255)
                                                   gggTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGt                                                                                                                             >  1:4118958/1‑124 (MQ=255)
                                                      tttACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgtaggagcgg                                                                                                >  1:1532312‑M2/1‑139 (MQ=255)
                                                       ttACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCGAAAACCGAGCGGGGGGGGACCCATGTTGGTGATTAGGCGCAGGCCGCGtggaggagcggg                                                                                               <  2:448141‑M2/150‑13 (MQ=255)
                                                        gttgtgcgaaaggGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAACCGCGCCTACAGCCGAGCGGGGGGGGACCCATGTTGGTGCGTTGGCGCAGGCCGCCggtcggagcggcg                                                                                              <  2:3195984‑M2/145‑14 (MQ=255)
                                                          cccGCGAACGGGCCGGGCCTGCGGGCCAATCTCTCCCCGTCGCTCCCCGTTGCCGCCGCTACCGTATTCGCGTGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgtaggagcgggttt                                                                                            >  1:1162873‑M2/1‑135 (MQ=255)
                                                           ccGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCGCCTACATCCGATCGGGGGTGGACACAGGTTGGTGATTAGGCGCAGGCCGCGgggaggagcgggttta                                                                                           <  2:4581084‑M2/150‑17 (MQ=255)
                                                           ccGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCGCCGAACACCGAGCGGGGGTGGACCCAAGTTGGTGATTAGGCGCCCGCCGCGtgtaggctcgggtttg                                                                                           <  2:3911654‑M2/150‑17 (MQ=255)
                                                           ccGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACATCCGAACGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgtaggagcgggttta                                                                                           <  2:3572517‑M2/150‑17 (MQ=255)
                                                             gcgAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACATCCCATCGGGGGTGGACACAAGTTGGTGATTAGGCGCAGGCCGCGtgtaggagcgggtttacc                                                                                         <  2:1893671‑M2/150‑19 (MQ=255)
                                                                aaaGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTg                                                                                                                               >  1:231575/1‑109 (MQ=255)
                                                                        gggccTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGggc                                                                              >  2:1332092/1‑150 (MQ=255)
                                                                             tgctgGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTgg                                                                                                                              >  2:5011559/1‑97 (MQ=255)
                                                                               ctgGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAACCCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCct                                                                       <  1:4701629/150‑1 (MQ=255)
                                                                               ctgGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCCGGTTTACCCGCGAAAGGGCCGGCCct                                                                       >  2:2700592/1‑150 (MQ=255)
                                                                                 ggCCAATCCCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTTTTCGCGGGTAACCCCGCTCCTACATCCGATCGGGGGTGGACACTTGTTGGTGATTAGGCGCCGGCCGCGTGTTGGAGGGGGTTTACCCGCGAAAGGGCCGGCCctgc                                                                     <  1:1332092/150‑1 (MQ=37)
                                                                                   ccAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCCCGCCTACATCCGATCGGGGGTGGACACAAGTTGGTGATTAGGCGCAGGCCGCGGGTTGGAGCGGGTTTACCCGCGAAAGGGCCGGCCcagctg                                                                   >  1:1281074/1‑150 (MQ=21)
                                                                                    cagtaactACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCCGGTAAAACCCCTCCTACCTCCGATAGGGGGTGGACACATGTTGGTTATTAGGCGCACGCCGCGTGGAGGAGCGGGTTTAACCGCGAAAGGGCCGGCCCAGCTgg                                                                  >  2:2484987/7‑150 (MQ=25)
                                                                                    cAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTTAACCCCCTCCTACCACCCATCGGGGGTGGGCACATGTTGGTGTTTATGGGCAGGCCGCGTGTAGGAGCGGGGTTTCCCGCGCAAGGGCCGGCCCTGCTgg                                                                  >  1:759412/1‑150 (MQ=25)
                                                                                    cAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCCGGTAAACCCCCTCCAACATCCCATCGGGGGTGGAAACAAGTTGGTGATTAGGCGCAGGCCGCGGGTAGGAGCGGGTTTAACCGCGAAAGGGACCGCCCAGCTgg                                                                  >  2:460021/1‑150 (MQ=25)
                                                                                      aTCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgtaggagcgggtaaacccgctc                                                                                    >  2:1230759‑M2/1‑107 (MQ=255)
                                                                                      aTCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgtaggagcgggtaaacccgctc                                                                                    <  1:1230759‑M2/130‑24 (MQ=255)
                                                                                        ctctACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGccca                                                              >  2:3496313/1‑148 (MQ=255)
                                                                                         tctACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGCCACATGTTGGTGATTAGGCGCAGGCCGCGTGGAGGAGCGGGTGTACCCGCGAAAGGGCCGGCCCTGCTGGCCAAc                                                             >  2:5250464/1‑150 (MQ=255)
                                                                                         actACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACAACCCATCGGGGGTGGACACATGTTGGTGAATAGGCGCAGGCCGCGTGTAGGAGCGGGGTTACCCGCCAAAGGGCCGGCCCAGCTGGCCAAc                                                             >  2:366681/2‑150 (MQ=21)
                                                                                           tACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACtt                                                           <  2:2280077/150‑1 (MQ=255)
                                                                                              ccGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGGAAACCCCCGCCGACCACCGAGCGGGGGTGGGCACATGTTGGGGCTTAGGCGCCGGCCGCGGGGAGGGGCGGGGGTACCCCCCCCAGGGCCCGCCCCGCGGGCCCAcggccc                                                        >  1:2564132/1‑145 (MQ=12)
                                                                                               cGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAACCCCGCTCCTACATCCGATGGGGGGTGGACACATGTTGGTGATTAGGCGGATGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCAc                                                       <  2:3795731/150‑1 (MQ=255)
                                                                                               cGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCAc                                                       <  2:2729140/150‑1 (MQ=255)
                                                                                                     tCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCGCCGACCGCCGCGCGGGGGGGGGCCCCGGGTGGGGGTTAGGCGCAGGCCGCGGGGGGGGGCGGGGGTACCCGCGCAAGGGCCGGCCCGGCGGGCCCACAGCCACcccccg                                                 >  1:4373687/1‑150 (MQ=12)
                                                                                                       tgcGTGGCCGCCTCTACGGTTCTCGGGGGTAATCCGGCTCCTCCATCCGTTGGGGGGTGTACCCATGTTGGTGATTCGGCGCACGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGtt                                               <  1:4140192/148‑1 (MQ=21)
                                                                                                           tggccgccGCTACCGCTGTCGCGGGTAACCGCGCTCCGTCAACCGATGGGGCGTGGACCCAAGATGCTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCAcc                                           <  2:681050/148‑1 (MQ=18)
                                                                                                               cgccgcTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCCTTg                                       >  1:765284/1‑150 (MQ=255)
                                                                                                                     tACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTATGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTgccgccc                                 >  2:2014273/1‑149 (MQ=255)
                                                                                                                     tACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCa                                 >  2:4752399/1‑150 (MQ=255)
                                                                                                                     tACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACAAGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCa                                 >  2:2638699/1‑150 (MQ=255)
                                                                                                                        cGTATTCGCGGGTAACCCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCacta                              <  1:3612448/150‑1 (MQ=255)
                                                                                                                           aTTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTc                           >  2:5028610/1‑150 (MQ=255)
                                                                                                                             tCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCaa                         >  1:4428100/1‑150 (MQ=255)
                                                                                                                               gcgGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAAc                       >  1:4535636/1‑150 (MQ=255)
                                                                                                                                cgtgTAACCCCTCTCCTTCATCCGATCGGGGGTTGTCAGATGTGGGTGTTTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTAGCCGCGAAAGGGCCTGCCCTGCTGGCCAACTTCCACCCACCGTTCACCGTTGCCGCCACTACTCaaagc                      <  2:5108844/147‑3 (MQ=18)
                                                                                                                                 gggTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAAccc                     >  2:3563363/1‑150 (MQ=255)
                                                                                                                                 gggTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACGTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAAccc                     <  2:3529173/150‑1 (MQ=255)
                                                                                                                                  gttAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGGGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCt                    <  1:9118/148‑1 (MQ=255)
                                                                                                                                   gTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCtt                   >  2:3159364/1‑150 (MQ=255)
                                                                                                                                    ttaccccGCTCCTTCATGCTATGGGGGGTTTACTCATGTTGGTGCTTCGGCGCATGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCttt                  <  1:5160900/146‑1 (MQ=18)
                                                                                                                                    tAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCttt                  <  1:5326881/150‑1 (MQ=255)
                                                                                                                                      aaCCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGGAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCa                <  2:4037162/150‑1 (MQ=255)
                                                                                                                                        cccGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgtaggagcgggtttacccgcgaatacggtagcggcggcaacggtgagcgacgggtagagattggccagca                                    <  1:4084932‑M2/128‑72 (MQ=255)
                                                                                                                                        cccGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGtgtaggagcgggtttacccgcgaatacggtagcggcggcaacggtgagcgacgggtagagattggccagca                                    >  2:4084932‑M2/1‑57 (MQ=255)
                                                                                                                                        cccGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCaaa              >  2:1199010/1‑150 (MQ=255)
                                                                                                                                         ccGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCaaaa             >  1:2689981/1‑150 (MQ=255)
                                                                                                                                          ggCTCCTACATGCGTTGGGGGGTGGACTCATGTTGGTGTTTAGGCGGAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCaaaaa            <  1:4672394/149‑1 (MQ=21)
                                                                                                                                          cGCTCCTACATCCGATCGGGGGTGGACAAATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCaaaaa            <  1:1918410/150‑1 (MQ=255)
                                                                                                                                           gCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAg           <  2:1020558/150‑1 (MQ=255)
                                                                                                                                             tCCTACATCCGATCGGGGGTGTACACATGTTGGTGATTAGGCGCAGGCCGCGtgtaggagcgggtttaccctcgaatactgtagaggcggcaacggtgagcgacgtgtatagattggccagcaggcacggcccattcgcgggtaaacccg         <  2:2837784‑M2/150‑99 (MQ=255)
                                                                                                                                             gccTACATCGGCTCGGGGGTGGACACTTGTTGGTGATTTGGCGCCGGCCGCGTGTCGGTGGGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGcc         <  2:1611916/149‑1 (MQ=21)
                                                                                                                                              ccTACATGCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGccc        <  1:1505884/150‑1 (MQ=255)
                                                                                                                                               cTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGcccc       >  1:1519426/1‑150 (MQ=255)
                                                                                                                                                 aCATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCAt     <  2:1484586/150‑1 (MQ=255)
                                                                                                                                                  cATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATg    >  1:4747206/1‑150 (MQ=255)
                                                                                                                                                    tCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCa  >  1:2071899/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCA  >  NC_002947/4595680‑4595977

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