Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,112,450 C→T 5.3% A76T (GCC→ACC)  aroP‑II ← aromatic amino acid transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,112,4500CT 5.3% 534.1 / 12.0 169A76T (GCC→ACC) aroP‑IIaromatic amino acid transport protein
Reads supporting (aligned to +/- strand):  ref base C (89/71);  new base T (6/3);  total (95/74)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 170 total aligned reads displayed.

GATCTCCGGCCACCAGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCACC  >  NC_002947/5112301‑5112597
                                                                                                                                                     |                                                                                                                                                   
gATCTCCGGCCACCAGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGc                                                                                                                                                     <  1:5179317/150‑1 (MQ=255)
          caccaGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGctct                                                                                                                                           <  2:2070504/150‑1 (MQ=255)
             caGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGACCACCGGCTCTTcg                                                                                                                                        >  1:5049547/1‑150 (MQ=255)
              agtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcga                                                                                                                                       <  1:3114457/150‑1 (MQ=255)
              agtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAAGAGCCGGTGGCCGGctcgttcag                                                                                                                                      <  1:1450446/150‑7 (MQ=255)
               gtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcgac                                                                                                                                      <  1:3847638/150‑1 (MQ=255)
               gtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcgac                                                                                                                                      <  1:2353175/150‑1 (MQ=255)
                 agtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCtttcagcc                                                                                                                                    <  1:5096037/150‑7 (MQ=255)
                    gAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCa                                                                                                                                 >  2:3756901/1‑150 (MQ=255)
                     aaCATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCAt                                                                                                                                >  2:1220677/1‑150 (MQ=255)
                     aaCATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCAt                                                                                                                                >  2:4497072/1‑150 (MQ=255)
                        aTACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGAt                                                                                                                                   <  2:1432243/144‑1 (MQ=255)
                         tACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCa                                                                                                                            <  1:1368807/150‑1 (MQ=255)
                         tACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCa                                                                                                                            <  1:1285361/150‑1 (MQ=255)
                             tGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCa                                                                                                                        <  1:2549399/150‑1 (MQ=255)
                                 gACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTg                                                                                                                    >  2:3480904/1‑150 (MQ=255)
                                  aCCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTgg                                                                                                                   <  1:4356440/150‑1 (MQ=255)
                                   ccgccgAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGgc                                                                                                                  <  2:4068934/150‑1 (MQ=255)
                                    cgccgAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGgcg                                                                                                                 >  1:4522297/1‑150 (MQ=255)
                                        gAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATg                                                                                                             <  1:241792/150‑1 (MQ=255)
                                              tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCTGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                       >  1:1149231/1‑150 (MQ=255)
                                              tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCCTGATCagg                                                                                                       >  1:4747811/1‑148 (MQ=255)
                                              tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                       >  1:4679590/1‑150 (MQ=255)
                                              tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                       >  1:4701364/1‑150 (MQ=255)
                                              tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                       >  2:2197598/1‑150 (MQ=255)
                                              tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                       >  1:3384544/1‑150 (MQ=255)
                                                cGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGaa                                                                                                     <  2:1115777/150‑1 (MQ=255)
                                                  aCATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGAGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAgg                                                                                                   <  2:4676730/150‑1 (MQ=255)
                                                   cATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGc                                                                                                  <  1:1002277/150‑1 (MQ=255)
                                                    ctCCCCACCAGGATATCCAGCCCCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCg                                                                                                 <  2:1217635/149‑1 (MQ=255)
                                                     tCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGa                                                                                                <  1:4839071/150‑1 (MQ=255)
                                                      ccccACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGAt                                                                                               <  1:261894/150‑1 (MQ=255)
                                                       cccACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATg                                                                                              >  2:3915751/1‑150 (MQ=255)
                                                          accaGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAg                                                                                           >  2:5127193/1‑150 (MQ=255)
                                                          accaGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAg                                                                                           >  2:5128064/1‑150 (MQ=255)
                                                            caGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGcc                                                                                         >  1:3944218/1‑150 (MQ=255)
                                                             aGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcgac                                                                                                                                      <  1:2363729/104‑1 (MQ=255)
                                                               gATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCa                                                                                      >  1:3432061/1‑150 (MQ=255)
                                                               gATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCa                                                                                      >  2:2641235/1‑150 (MQ=255)
                                                                  ataCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGAt                                                                                   >  2:1991413/1‑150 (MQ=255)
                                                                  ataCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGAt                                                                                   >  2:937373/1‑150 (MQ=255)
                                                                   taCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATc                                                                                  <  2:2343491/150‑1 (MQ=255)
                                                                    aCAGCACCCAGCAGTTCCAGCCCGACAGGCAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCCTGATCATGACGGCGATGAAGCCGCAGAGcg                                                                                 >  2:2982514/1‑150 (MQ=255)
                                                                    aCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAACTGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGACGCCGCAGATcg                                                                                 >  2:2982557/1‑150 (MQ=255)
                                                                    aCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCAc                                                                                                                                                  <  1:1370857/85‑1 (MQ=255)
                                                                    aCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCAc                                                                                                                                                  >  2:1370857/1‑85 (MQ=255)
                                                                    aCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATc                                                                                                          >  1:5425199/1‑125 (MQ=255)
                                                                    aCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcg                                                                                 >  1:2024595/1‑150 (MQ=255)
                                                                    aCAGCACCCAGCAGGTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTTGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcg                                                                                 <  2:2368564/150‑1 (MQ=255)
                                                                     cAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcgc                                                                                >  2:239164/1‑150 (MQ=255)
                                                                     cAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcgc                                                                                >  2:239104/1‑150 (MQ=255)
                                                                     cAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcgc                                                                                >  1:272546/1‑150 (MQ=255)
                                                                      aGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAGGGCGATGAAGCCGCAGATcgcg                                                                               >  2:3531182/1‑150 (MQ=255)
                                                                        cACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTa                                                                             <  1:1823395/150‑1 (MQ=255)
                                                                        cACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTa                                                                             <  1:2182949/150‑1 (MQ=255)
                                                                          cccAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGc                                                                           <  1:2787334/150‑1 (MQ=255)
                                                                           ccAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGcc                                                                          >  2:532357/1‑150 (MQ=255)
                                                                             cggagTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGTCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGccta                                                                        <  2:4066248/147‑3 (MQ=255)
                                                                             agcagTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCaa                                                                        >  2:689807/1‑150 (MQ=255)
                                                                                    ccAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGatcat                                                                 <  2:1533598/150‑1 (MQ=255)
                                                                                        cccGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGAc                                                             <  1:1049277/150‑1 (MQ=255)
                                                                                          cGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACgg                                                           >  2:4967256/1‑150 (MQ=255)
                                                                                           gACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACggg                                                          >  2:4913638/1‑150 (MQ=255)
                                                                                                gAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCcggc                                                     >  2:4795455/1‑150 (MQ=255)
                                                                                                  aGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCcggcgg                                                   <  2:4202115/150‑1 (MQ=255)
                                                                                                  aGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGCCAGAACGAGCCGGCCACCGGCTCTACGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCAACGACGGGCcggcgg                                                   >  1:1478475/1‑150 (MQ=255)
                                                                                                     cGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGAtt                                                <  2:3829675/150‑1 (MQ=255)
                                                                                                      ggCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGAttt                                               >  2:3717272/1‑150 (MQ=255)
                                                                                                        cGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTtca                                             <  1:377231/150‑1 (MQ=255)
                                                                                                         gAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTtcat                                            >  2:1674659/1‑150 (MQ=255)
                                                                                                          cgcccacccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTtcatc                                           >  1:2491796/3‑150 (MQ=255)
                                                                                                             ccacccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACg                                        >  1:5234232/1‑150 (MQ=255)
                                                                                                               acccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGcc                                      >  2:5156113/1‑150 (MQ=255)
                                                                                                                ccccAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCTACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCt                                     >  1:3605023/1‑150 (MQ=255)
                                                                                                                 cccAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTg                                    >  2:357721/1‑150 (MQ=255)
                                                                                                                         gtgtGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTgccgagccg                            >  1:2985729/1‑150 (MQ=255)
                                                                                                                          tgtgGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTgccgagccga                           <  1:1736113/150‑1 (MQ=255)
                                                                                                                           gtgGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGccgagccgag                          <  1:3717272/150‑1 (MQ=255)
                                                                                                                            tgGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAgg                         <  1:544302/150‑1 (MQ=255)
                                                                                                                            tgGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAgg                         <  1:81107/150‑1 (MQ=255)
                                                                                                                              ggCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGaa                       >  1:951246/1‑150 (MQ=255)
                                                                                                                              ggCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGaa                       >  2:4343890/1‑150 (MQ=255)
                                                                                                                                  aaaTGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAgg                   <  1:1665967/150‑1 (MQ=255)
                                                                                                                                    aTGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGcc                 <  2:310529/150‑1 (MQ=255)
                                                                                                                                       gCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGatcat                                                                 >  1:5141220/1‑99 (MQ=255)
                                                                                                                                       gCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGatcat                                                                 <  2:5141224/99‑1 (MQ=255)
                                                                                                                                        cTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTa             <  2:1860410/150‑1 (MQ=255)
                                                                                                                                        cTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTa             <  2:638023/150‑1 (MQ=255)
                                                                                                                                        cTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTa             <  2:638254/150‑1 (MQ=255)
                                                                                                                                         tcgaagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGACCAGGCCGGTAc            >  1:4575211/6‑150 (MQ=38)
                                                                                                                                          cgaagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAgg                         >  2:4378286/5‑136 (MQ=255)
                                                                                                                                          cgaagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAgg                         <  1:4378286/132‑1 (MQ=255)
                                                                                                                                           gaagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCa          >  2:3623266/4‑150 (MQ=255)
                                                                                                                                           aaCGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCa          >  2:4607737/1‑150 (MQ=255)
                                                                                                                                           aaCGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCa          >  1:2878439/1‑150 (MQ=255)
                                                                                                                                           aaCGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGACCGGTACCa          >  2:4517984/1‑150 (MQ=255)
                                                                                                                                            aagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCaa         >  2:5122370/3‑150 (MQ=255)
                                                                                                                                            aagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCaa         >  1:255970/3‑150 (MQ=255)
                                                                                                                                              gAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATc       >  2:3267678/1‑150 (MQ=255)
                                                                                                                                                 ccGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCa    >  2:1920911/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GATCTCCGGCCACCAGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCACC  >  NC_002947/5112301‑5112597

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