Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,653,241 G→C 5.1% intergenic (+41/+29) PP_2323 → / ← aroF‑I hypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,653,2410GC 5.1% 675.5 / 19.7 197intergenic (+41/+29)PP_2323/aroF‑Ihypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (71/115);  new base C (7/3);  total (78/119)
Fisher's exact test for biased strand distribution p-value = 9.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 199 total aligned reads displayed.

CGAGCGGGCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACTGCAAATCG  >  NC_002947/2653092‑2653390
                                                                                                                                                     |                                                                                                                                                     
    cGGGCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGcc                                                                                                                                                   <  2:938378/150‑1 (MQ=255)
      ggCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCg                                                                                                                                                 >  1:2373372/1‑150 (MQ=255)
         tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCg                                                                                                                                              >  1:2687912/1‑150 (MQ=255)
         tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCg                                                                                                                                              >  2:1669995/1‑150 (MQ=255)
         tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCg                                                                                                                                              >  1:1747145/1‑150 (MQ=255)
           tGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGtt                                                                                                                                            <  1:189232/150‑1 (MQ=255)
            ggAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGttt                                                                                                                                           <  2:4496060/150‑1 (MQ=255)
            ggAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGttt                                                                                                                                           >  1:2154591/1‑150 (MQ=255)
            ggAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGttt                                                                                                                                           <  2:1163840/150‑1 (MQ=255)
            ggAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttt                                                                                                                                           <  2:1348415/150‑1 (MQ=255)
               gTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTc                                                                                                                                        <  2:2334585/150‑1 (MQ=255)
               gTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTc                                                                                                                                        >  2:1939935/1‑150 (MQ=255)
                tACGCCGAGCAGTTGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGAAAGGCCTGAGCTGACGACCGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCa                                                                                                                                       <  1:3068491/150‑1 (MQ=255)
                tACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCa                                                                                                                                       >  1:3734230/1‑150 (MQ=255)
                   gCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTc                                                                                                                                    >  1:3843462/1‑150 (MQ=255)
                    ccGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAATCGCCGGGCAGTGCCCGGCGTTTTTGAttgc                                                                                                                                   >  1:619301/1‑148 (MQ=255)
                     cGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTccc                                                                                                                                  >  1:4488214/1‑150 (MQ=255)
                      gAGCAGTTGGGTTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGATTTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCg                                                                                                                                 <  1:1191516/150‑1 (MQ=255)
                          aGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAAt                                                                                                                             <  2:686369/150‑1 (MQ=255)
                             tGGGGTACACGGTGTTCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGtt                                                                                                                          <  2:3939143/150‑1 (MQ=255)
                             tGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGtt                                                                                                                          <  1:3728669/150‑1 (MQ=255)
                                ggTACACGGTGATCCCGAATTGCTCATATGTGGAGCGTTACATGGAGCCCCAGCAGCCCCAATCGACCAAGCCCTCAGCTGACGACAGGCACCAGACACAACAAAAAACGCCCGGCAGTGCCCGGCGAATTTCATTCCCGCAATGTTCAg                                                                                                                       <  2:619301/150‑1 (MQ=255)
                                         ttaTCCCGTCTTGCTAATTTGTGGAGCGTTACTTGGAGCGCCAGCAGCGGCATTCGAGAATGGCCTGAGCTGACGACAGGCACCACACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTa                                                                                                              <  1:1110524/148‑1 (MQ=255)
                                          gATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTAc                                                                                                             <  2:3181159/150‑1 (MQ=255)
                                                   ttGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgca                                                                                                    >  2:876713/1‑150 (MQ=255)
                                                   ttGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgca                                                                                                    >  1:2109511/1‑150 (MQ=255)
                                                    tGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgcag                                                                                                   <  1:1841373/150‑1 (MQ=255)
                                                       tCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCa                                                                                                                                       <  1:1352813/111‑1 (MQ=255)
                                                          tatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACACCAAACAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCACTGTTCAGCCGCGCTTACGCTGCGGCAGCACATc                                                                                             <  2:2626078/150‑1 (MQ=255)
                                                           atGTGGAGCGTTACCTGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACACCAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATcc                                                                                            <  1:316250/150‑1 (MQ=255)
                                                             gtgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                          <  1:479172/150‑1 (MQ=255)
                                                             gtgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                          >  2:299489/1‑150 (MQ=255)
                                                             gtgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                          <  2:1542133/150‑1 (MQ=255)
                                                              tgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTg                                                                                         <  1:116796/150‑1 (MQ=255)
                                                              gggAGCGTTACCTGGCGCGCCAGCCGCGGCATTCGAGCACGGCCTGAGCTGACGACCGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTg                                                                                         <  1:2357450/149‑1 (MQ=255)
                                                               ggAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGa                                                                                        >  2:1375231/1‑150 (MQ=255)
                                                               ggAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGa                                                                                        <  1:1387812/150‑1 (MQ=255)
                                                               ggAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGa                                                                                        >  2:2543323/1‑150 (MQ=255)
                                                                    gTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTg                                                                                   <  1:467154/150‑1 (MQ=255)
                                                                        cATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcat                                                                               <  2:2941348/150‑1 (MQ=255)
                                                                         aTGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatg                                                                              <  2:4418133/150‑1 (MQ=255)
                                                                         aTGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatg                                                                              >  2:1083986/1‑150 (MQ=255)
                                                                          tGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGCGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatgc                                                                             <  1:322197/150‑1 (MQ=255)
                                                                                gacaggagCGGGATTCGTGGTAGGCCTGATCTGACGAGAGGCACCAGACACAACAAAAATCGCCGGGCATTGCCCGGAGATTTTCATTCCCGCAATGTTCAGCCGAGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCAtgctg                                                                       <  2:3964353/148‑1 (MQ=255)
                                                                                 ccAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCAtgctgc                                                                      >  1:3530354/1‑150 (MQ=255)
                                                                                  cagcagCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTgcg                                                                     >  1:3289955/1‑150 (MQ=255)
                                                                                  cagcagCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTgcg                                                                     >  1:678469/1‑150 (MQ=255)
                                                                                      agcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAgg                                                                 <  2:2154591/150‑1 (MQ=255)
                                                                                      agcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAgg                                                                 <  2:178728/150‑1 (MQ=255)
                                                                                       gcGGCTTTCGAGCAAGGCCTGAGCTGACGACCGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAggg                                                                <  1:71591/150‑1 (MQ=255)
                                                                                        cGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGt                                                               <  1:755459/150‑1 (MQ=255)
                                                                                         ggCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGtt                                                              <  2:1253909/150‑1 (MQ=255)
                                                                                           cATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGtttt                                                            <  2:2202305/150‑1 (MQ=255)
                                                                                             ttCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtct                                                          >  2:4105705/1‑150 (MQ=255)
                                                                                             ttCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtct                                                          >  1:4243876/1‑150 (MQ=255)
                                                                                              tCGAGCAAGGCCTGCGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtctc                                                         <  2:1755179/150‑1 (MQ=255)
                                                                                              tCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtctc                                                         <  1:3172304/150‑1 (MQ=255)
                                                                                              tCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtctc                                                         <  2:64662/150‑1 (MQ=255)
                                                                                               cGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCg                                                        <  2:574438/150‑1 (MQ=255)
                                                                                                gAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCgg                                                       >  2:1830027/1‑150 (MQ=255)
                                                                                                gAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCgg                                                       <  1:4492610/150‑1 (MQ=255)
                                                                                                 aGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGt                                                      >  2:1333580/1‑150 (MQ=255)
                                                                                                  gCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTa                                                     <  1:1669995/150‑1 (MQ=255)
                                                                                                  gCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTa                                                     <  1:1516615/150‑1 (MQ=255)
                                                                                                    aaGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGc                                                   >  1:4487428/1‑150 (MQ=255)
                                                                                                     aGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGGGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGcc                                                  <  1:2554810/150‑1 (MQ=255)
                                                                                                      ggCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCg                                                 >  1:2207771/1‑150 (MQ=255)
                                                                                                       gCCTGAGCTGACGACTGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGa                                                <  1:96366/150‑1 (MQ=255)
                                                                                                       gCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGa                                                <  1:3638060/150‑1 (MQ=255)
                                                                                                       gCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGa                                                <  1:1372352/150‑1 (MQ=255)
                                                                                                          tGAGCTGACGACCGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCa                                             <  2:1630465/150‑1 (MQ=255)
                                                                                                          tGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCa                                             <  1:1584659/150‑1 (MQ=255)
                                                                                                          tGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCa                                             <  1:1585137/150‑1 (MQ=255)
                                                                                                                   gacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTATCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCGGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCa                                    >  2:1404156/1‑150 (MQ=255)
                                                                                                                    acAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                                   <  1:4109902/150‑1 (MQ=255)
                                                                                                                    acAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                                   <  1:2761295/150‑1 (MQ=255)
                                                                                                                    acAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                                   <  2:3933526/150‑1 (MQ=255)
                                                                                                                    acAGCCACCACACACAACAAAAAACCCCGGGCAGTCCCCGGCGTTTTCCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                                   <  1:2462964/150‑1 (MQ=255)
                                                                                                                      aGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGc                                 >  1:1943517/1‑150 (MQ=255)
                                                                                                                        gCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCGTGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGt                               >  1:2061565/1‑150 (MQ=255)
                                                                                                                         cACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTc                              <  2:2807125/150‑1 (MQ=255)
                                                                                                                         cACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTc                              <  2:6516/150‑1 (MQ=255)
                                                                                                                          aCCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTcc                             <  1:1317876/150‑1 (MQ=255)
                                                                                                                           ccAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCg                            >  1:2117256/1‑150 (MQ=255)
                                                                                                                            cAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGt                           <  1:3058333/150‑1 (MQ=255)
                                                                                                                             aGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTg                          >  2:802683/1‑150 (MQ=255)
                                                                                                                                      aaaaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACg                 >  1:4329113/1‑150 (MQ=255)
                                                                                                                                      aaaaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACg                 >  2:54644/1‑150 (MQ=255)
                                                                                                                                       aaaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATACTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGc                >  1:3972852/1‑150 (MQ=255)
                                                                                                                                        aaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGcc               >  2:48873/1‑150 (MQ=255)
                                                                                                                                        aaaaCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGAAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGcc               >  2:4038248/1‑150 (MQ=255)
                                                                                                                                         aaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCa              >  1:1017770/1‑150 (MQ=255)
                                                                                                                                         aaaCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgc                                                                                                     <  1:3396148/63‑1 (MQ=255)
                                                                                                                                         aaaCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgc                                                                                                     >  2:3396148/1‑63 (MQ=255)
                                                                                                                                           aCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATa            >  2:3478811/1‑150 (MQ=255)
                                                                                                                                           aCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATa            <  2:1051595/150‑1 (MQ=255)
                                                                                                                                           aCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATa            >  2:2216740/1‑150 (MQ=255)
                                                                                                                                            cGCCGGTCAGTGCCCGTCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATAc           <  2:4488201/150‑1 (MQ=255)
                                                                                                                                             gCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACt          <  2:3069224/150‑1 (MQ=255)
                                                                                                                                                  gCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACTGCaaa     >  2:2346088/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAGCGGGCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACTGCAAATCG  >  NC_002947/2653092‑2653390

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