Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,654,931 G→C 6.0% intergenic (+37/‑221) 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,9310GC 6.0% 308.6 / 9.2 103intergenic (+37/‑221)aceA/PP_4118isocitrate lyase/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (43/52);  new base C (4/2);  total (48/55)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 104 total aligned reads displayed.

AGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGC  >  NC_002947/4654785‑4655080
                                                                                                                                                  |                                                                                                                                                      
agcagGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGaaa                                                                                                                                                    <  2:1559055/150‑1 (MQ=255)
agcagGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTNCCACTGNTAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGnaa                                                                                                                                                    <  1:1935981/150‑1 (MQ=255)
 gcagGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAg                                                                                                                                                   <  1:348017/150‑1 (MQ=255)
   agGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCt                                                                                                                                                 >  2:1975453/1‑150 (MQ=255)
     gAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTcc                                                                                                                                               <  2:234605/150‑1 (MQ=255)
      aaGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTATGAACGcgcta                                                                                                                                              >  2:1957575/1‑146 (MQ=255)
       aGTTGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCgg                                                                                                                                             <  1:957503/150‑1 (MQ=255)
       aGTGGGCACTGGCTACTTCGACGACATGCCCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCgg                                                                                                                                             <  2:1073681/150‑1 (MQ=255)
       aGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCgg                                                                                                                                             <  1:445965/150‑1 (MQ=255)
           ggCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCt                                                                                                                                         >  2:326625/1‑150 (MQ=255)
             cACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCttt                                                                                                                                       <  2:1929418/150‑1 (MQ=255)
               cTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGt                                                                                                                                     <  2:1049917/150‑1 (MQ=255)
                tGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGtt                                                                                                                                    <  1:326625/150‑1 (MQ=255)
                   cTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGttttt                                                                                                                                 >  2:1476261/1‑150 (MQ=255)
                    tACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTa                                                                                                                                <  1:2069704/150‑1 (MQ=255)
                     aCTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAc                                                                                                                               <  2:1340493/150‑1 (MQ=255)
                      cTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAcc                                                                                                                              <  1:1975453/150‑1 (MQ=255)
                       ttCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCGTGTGAAAGCTCGGGGCCTTTGTTTTTAccc                                                                                                                             >  2:582657/1‑150 (MQ=255)
                        tCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCg                                                                                                                            <  1:1954069/150‑1 (MQ=255)
                                aTGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGTACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAggg                                                                                                                    >  2:106601/1‑150 (MQ=255)
                                  gACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggcc                                                                                                                  <  1:2241202/150‑1 (MQ=255)
                                  gACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTAACCGGTTCGACCGCAGGAGAGTAGTTCTAATGATAGTGGGTGGCCGTGGGACAAATGCCCGGCGCTGGTGAAGGCTCCGGGCGTGTGTGTTTCGCCGAGCAAgggcc                                                                                                                  >  2:226706/1‑150 (MQ=255)
                                    ccaccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgg                                                                                                                <  2:1828319/150‑1 (MQ=255)
                                    ccaccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCAATGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgg                                                                                                                <  1:961243/150‑1 (MQ=255)
                                     caccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccggg                                                                                                               >  1:1422290/1‑150 (MQ=255)
                                     caccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccggg                                                                                                               >  1:1394175/1‑150 (MQ=255)
                                      accGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgggc                                                                                                              >  1:1411844/1‑150 (MQ=255)
                                        cGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCttt                                                                                                                                       <  2:610016/123‑1 (MQ=37)
                                        cGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCttt                                                                                                                                       >  1:610016/1‑123 (MQ=37)
                                          tGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCttt                                                                                                                                       <  2:279112/121‑1 (MQ=37)
                                          tGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCttt                                                                                                                                       >  1:279112/1‑121 (MQ=37)
                                            aTCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgctg                                                                                                        <  2:2142068/150‑1 (MQ=255)
                                              ccAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGc                                                                                                      <  2:2163325/150‑1 (MQ=255)
                                              ccAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGc                                                                                                      <  2:487568/150‑1 (MQ=255)
                                                aGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCa                                                                                                    >  1:2037843/1‑150 (MQ=255)
                                                      gcgcTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAttttt                                                                                              <  2:1422290/150‑1 (MQ=255)
                                                         cTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTAcc                                                                                           <  1:2437414/150‑1 (MQ=255)
                                                                  tGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCa                                                                                  <  1:1196542/150‑1 (MQ=255)
                                                                     cTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCgg                                                                                                                                             <  1:1708815/88‑1 (MQ=255)
                                                                     cTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCgg                                                                                                                                             >  2:1708815/1‑88 (MQ=255)
                                                                          cTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCCGCTGGCCAATTTTTACCAGCCGTTCACCGTTTcc                                                                          <  1:2008143/150‑1 (MQ=255)
                                                                           tGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTccg                                                                         >  2:971112/1‑150 (MQ=255)
                                                                             aCCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAc                                                                                                                               >  1:1007419/1‑94 (MQ=255)
                                                                             aCCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAc                                                                                                                               <  2:1007419/94‑1 (MQ=255)
                                                                              ccGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTccgccg                                                                      >  1:1171821/1‑150 (MQ=255)
                                                                                     gACCAAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATGTTTACCAGCCGTTCACCGTTTCCGCGGCTACCGt                                                               >  2:1309290/1‑150 (MQ=255)
                                                                                      aCCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTa                                                              >  1:1714287/1‑150 (MQ=255)
                                                                                       ccGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCt                                                                                                                                         <  2:2276340/74‑1 (MQ=255)
                                                                                       ccGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCt                                                                                                                                         >  1:2276340/1‑74 (MQ=255)
                                                                                        cGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGATTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTAtt                                                            >  1:158695/1‑150 (MQ=255)
                                                                                           agaagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcgc                                                         >  1:991850/1‑150 (MQ=255)
                                                                                            gaagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGTCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGATAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcgcg                                                        <  2:262058/150‑1 (MQ=255)
                                                                                            gaagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcgcg                                                        <  1:106601/150‑1 (MQ=255)
                                                                                             aagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCgg                                                       <  1:335853/150‑1 (MQ=255)
                                                                                                     ttCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGtt                                                                                                                                    >  1:124817/1‑65 (MQ=255)
                                                                                                     ttCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGtt                                                                                                                                    <  2:124817/65‑1 (MQ=255)
                                                                                                     ttCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAAccc                                               >  1:967304/1‑150 (MQ=255)
                                                                                                       ccACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGc                                             <  1:1470179/150‑1 (MQ=255)
                                                                                                        cACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCa                                            <  2:967304/150‑1 (MQ=255)
                                                                                                        cACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCa                                            >  1:1721492/1‑150 (MQ=255)
                                                                                                         aCTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTTAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCAc                                           >  1:458934/1‑150 (MQ=255)
                                                                                                         aCTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCa                                                                                          <  2:164572/103‑1 (MQ=255)
                                                                                                         aCTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCa                                                                                          >  1:164572/1‑103 (MQ=255)
                                                                                                         aCTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCGAGTCACCGTTTCCGCCGCTACCGAATTCGCGGGTAAACCCGCAc                                           >  1:1095448/1‑150 (MQ=255)
                                                                                                             aTAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTaaa                                                  >  1:1754221/1‑139 (MQ=255)
                                                                                                             aTAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTaaa                                                  <  2:1754221/139‑1 (MQ=255)
                                                                                                               aGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTacac                                     <  1:971112/150‑1 (MQ=255)
                                                                                                                     ttGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCg                               <  2:1394175/150‑1 (MQ=255)
                                                                                                                          gTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgg                                                                                                                >  2:1767176/1‑64 (MQ=255)
                                                                                                                           ttGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggagga                         >  1:120132/1‑150 (MQ=255)
                                                                                                                             gAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAtttt                                                                                               <  2:592609/78‑1 (MQ=255)
                                                                                                                             gAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAtttt                                                                                               >  1:592609/1‑78 (MQ=255)
                                                                                                                             gAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                       <  1:1476261/150‑1 (MQ=255)
                                                                                                                             gAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                       <  2:1714287/150‑1 (MQ=255)
                                                                                                                             aaaCAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                       >  2:1412509/2‑150 (MQ=25)
                                                                                                                             aaaCAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                       >  1:1866027/2‑150 (MQ=21)
                                                                                                                              aaCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTc                                                                           <  1:1453579/97‑1 (MQ=255)
                                                                                                                              aaCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTc                                                                           >  2:1453579/1‑97 (MQ=255)
                                                                                                                              aaCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCgg                                                       <  2:1734271/117‑1 (MQ=255)
                                                                                                                              aaCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCgg                                                       >  1:1734271/1‑117 (MQ=255)
                                                                                                                              aaCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                       >  2:1378352/1‑149 (MQ=255)
                                                                                                                              aaCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGgaga                       <  1:1378352/149‑1 (MQ=255)
                                                                                                                               aCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGaca                     >  1:763233/1‑150 (MQ=255)
                                                                                                                               aCAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGaca                     <  1:1412509/150‑1 (MQ=25)
                                                                                                                                 aaaGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTAc                                                                                            >  1:109169/1‑77 (MQ=255)
                                                                                                                                 aaaGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTAc                                                                                            <  2:109169/77‑1 (MQ=255)
                                                                                                                                  aaGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAcc                                                                                >  1:266931/1‑88 (MQ=255)
                                                                                                                                  aaGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAcc                                                                                <  2:266931/88‑1 (MQ=255)
                                                                                                                                   aGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGAcacaca                 >  2:467352/1‑150 (MQ=255)
                                                                                                                                     gCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAc                                                                                 <  2:1480631/84‑1 (MQ=255)
                                                                                                                                     gCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAc                                                                                 >  1:1480631/1‑84 (MQ=255)
                                                                                                                                     gCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATc               <  1:467352/150‑1 (MQ=255)
                                                                                                                                       ccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggc                                                                                                                   >  2:353139/1‑48 (MQ=255)
                                                                                                                                       ccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggc                                                                                                                   <  1:353139/48‑1 (MQ=255)
                                                                                                                                       ccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCgg             >  1:2394817/1‑150 (MQ=255)
                                                                                                                                       ccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCgg             >  2:485347/1‑150 (MQ=255)
                                                                                                                                             cTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTaca                                      >  1:1028338/1‑119 (MQ=255)
                                                                                                                                             cTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTaca                                      <  2:1028338/119‑1 (MQ=255)
                                                                                                                                             cTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTg       >  1:2237427/1‑150 (MQ=255)
                                                                                                                                                  gAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGgcgc  >  1:779857/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
AGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGC  >  NC_002947/4654785‑4655080

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: