Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,112,452 A→G 5.2% V75A (GTG→GCG) ‡ 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,4520AG 5.2% 556.2 / 11.6 172V75A (GTG→GCG) ‡aroP‑IIaromatic amino acid transport protein
Reads supporting (aligned to +/- strand):  ref base A (91/72);  new base G (6/3);  total (97/75)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 174 total aligned reads displayed.

TCCGGCCACCAGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCACCGCCG  >  NC_002947/5112305‑5112601
                                                                                                                                                   |                                                                                                                                                     
tcCGGCCACCAGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCAcc                                                                                                                                                     >  1:2388047/1‑150 (MQ=255)
  cGGCCACCAGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCgg                                                                                                                                                   >  1:2664590/1‑150 (MQ=255)
       accaGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGTGGCCGGctcgt                                                                                                                                              <  2:2056443/150‑3 (MQ=255)
         caGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGACCACCGGCTCTTcg                                                                                                                                            >  1:5049547/1‑150 (MQ=255)
          agtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcga                                                                                                                                           <  2:5284936/150‑1 (MQ=255)
          agtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcga                                                                                                                                           <  1:3114457/150‑1 (MQ=255)
          agtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAAGAGCCGGTGGCCGGctcgttcag                                                                                                                                          <  1:1450446/150‑7 (MQ=255)
           gtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcgac                                                                                                                                          <  1:2353175/150‑1 (MQ=255)
                gAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCa                                                                                                                                     >  2:3756901/1‑150 (MQ=255)
                 aaCATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCAt                                                                                                                                    >  2:1220677/1‑150 (MQ=255)
                 aaCATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCAt                                                                                                                                    >  2:4497072/1‑150 (MQ=255)
                     tACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCa                                                                                                                                <  1:1368807/150‑1 (MQ=255)
                     tACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCa                                                                                                                                <  1:1285361/150‑1 (MQ=255)
                       cTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCAcc                                                                                                                              >  2:5063378/1‑150 (MQ=255)
                        ttGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCAccc                                                                                                                             >  2:3750431/1‑150 (MQ=255)
                             gACCGCCGAAAGCTCCGAAATCACCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACACACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTg                                                                                                                        >  2:3480985/1‑150 (MQ=255)
                              aCCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTgg                                                                                                                       <  1:2156335/150‑1 (MQ=255)
                              aCCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTgg                                                                                                                       >  2:1665967/1‑150 (MQ=255)
                              aCCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTgg                                                                                                                       <  1:4356627/150‑1 (MQ=255)
                               ccgccgAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGgc                                                                                                                      <  2:4068934/150‑1 (MQ=255)
                                cgccgAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGgcg                                                                                                                     >  1:4522297/1‑150 (MQ=255)
                                 gccgAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGgcgc                                                                                                                    <  2:1008298/150‑1 (MQ=255)
                                   cgAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCAt                                                                                                                  >  1:3822728/1‑150 (MQ=255)
                                       aGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATc                                                                                                              <  1:2591732/150‑1 (MQ=255)
                                          tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCTGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                           >  1:1149231/1‑150 (MQ=255)
                                          tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCCTGATCagg                                                                                                           >  1:4747811/1‑148 (MQ=255)
                                          tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                           >  1:2432852/1‑150 (MQ=255)
                                          tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                           >  1:1492389/1‑150 (MQ=255)
                                          tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                           >  1:3384562/1‑150 (MQ=255)
                                          tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                           >  1:4701364/1‑150 (MQ=255)
                                             gACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAg                                                                                                        <  1:3750431/150‑1 (MQ=255)
                                              aCATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAgg                                                                                                       >  2:4347564/1‑150 (MQ=255)
                                              aCATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGAGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAgg                                                                                                       <  2:4676730/150‑1 (MQ=255)
                                               cATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGc                                                                                                      <  1:1002277/150‑1 (MQ=255)
                                                aTCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCg                                                                                                     >  1:3979877/1‑150 (MQ=255)
                                                 tCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGa                                                                                                    <  1:4621641/150‑1 (MQ=255)
                                                 tCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGa                                                                                                    <  1:4839071/150‑1 (MQ=255)
                                                 tCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGa                                                                                                    <  1:2555670/150‑1 (MQ=255)
                                                   cccACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATg                                                                                                  >  2:1736113/1‑150 (MQ=255)
                                                      accaGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAg                                                                                               >  2:5127201/1‑150 (MQ=255)
                                                        caGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGcc                                                                                             >  1:3944218/1‑150 (MQ=255)
                                                           gATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCa                                                                                          >  2:2641235/1‑150 (MQ=255)
                                                           gATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCa                                                                                          <  2:1155589/150‑1 (MQ=255)
                                                              ataCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGAt                                                                                       >  2:1991413/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)
                                                                aCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATc                                                                                                              >  1:5425199/1‑125 (MQ=255)
                                                                aCAGCACCCAGCAGGTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTTGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcg                                                                                     <  2:2368564/150‑1 (MQ=255)
                                                                 cAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcgc                                                                                    >  2:239164/1‑150 (MQ=255)
                                                                 cAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcgc                                                                                    >  1:272546/1‑150 (MQ=255)
                                                                 cAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcgc                                                                                    <  1:4880728/150‑1 (MQ=255)
                                                                  aGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAGGGCGATGAAGCCGCAGATcgcg                                                                                   >  2:3531182/1‑150 (MQ=255)
                                                                    cACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTa                                                                                 <  1:1823395/150‑1 (MQ=255)
                                                                    cACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTa                                                                                 <  1:2182949/150‑1 (MQ=255)
                                                                       ccAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGcc                                                                              >  2:680937/1‑150 (MQ=255)
                                                                         cggagTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGTCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGccta                                                                            <  2:4066248/147‑3 (MQ=255)
                                                                         agcagTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCaa                                                                            >  2:689807/1‑150 (MQ=255)
                                                                          gcagTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAg                                                                           >  1:811111/1‑150 (MQ=255)
                                                                                ccAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCa                                                                                                                                <  1:361437/91‑1 (MQ=255)
                                                                                ccAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGatcat                                                                     <  2:1533598/150‑1 (MQ=255)
                                                                                    cccGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGAc                                                                 <  1:1049277/150‑1 (MQ=255)
                                                                                    cccGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGAc                                                                 <  1:3132418/150‑1 (MQ=255)
                                                                                            gAAGCTGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCcggc                                                         >  2:1327641/1‑150 (MQ=255)
                                                                                              aGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCcggcgg                                                       <  2:4202115/150‑1 (MQ=255)
                                                                                                 cGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGAtt                                                    <  2:3829675/150‑1 (MQ=255)
                                                                                                    cGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTtca                                                 <  2:1411672/150‑1 (MQ=255)
                                                                                                    cGAACCCACCCCAGTAGGTGGGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTtca                                                 <  1:5137792/150‑1 (MQ=255)
                                                                                                     gAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTtcat                                                >  2:3940104/1‑150 (MQ=255)
                                                                                                      cgcccacccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTtcatc                                               >  1:2491796/3‑150 (MQ=255)
                                                                                                        cccacccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCAc                                             <  1:2945924/150‑1 (MQ=255)
                                                                                                         ccacccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACg                                            >  1:5234232/1‑150 (MQ=255)
                                                                                                         ccacccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACg                                            >  2:4652695/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)
                                                                                                             cccAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTg                                        <  2:2106819/150‑1 (MQ=255)
                                                                                                                 gTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTgccga                                    <  2:2388047/150‑1 (MQ=255)
                                                                                                                      tgtgGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTgccgagccga                               <  1:1736113/150‑1 (MQ=255)
                                                                                                                        tgGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAgg                             <  1:544302/150‑1 (MQ=255)
                                                                                                                        tgGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAgg                             >  2:4565216/1‑150 (MQ=255)
                                                                                                                          ggCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGaa                           >  1:951246/1‑150 (MQ=255)
                                                                                                                              aaaTGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAgg                       <  1:1665967/150‑1 (MQ=255)
                                                                                                                                aTGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGcc                     <  2:310529/150‑1 (MQ=255)
                                                                                                                                   gCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGatcat                                                                     <  2:5141220/99‑1 (MQ=255)
                                                                                                                                   gCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGatcat                                                                     >  1:5141224/1‑99 (MQ=255)
                                                                                                                                   gCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGt                  >  2:4884502/1‑150 (MQ=255)
                                                                                                                                    cTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTa                 <  2:4927076/150‑1 (MQ=255)
                                                                                                                                    cTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTa                 <  2:638023/150‑1 (MQ=255)
                                                                                                                                     tcgaagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGACCAGGCCGGTAc                >  1:4575211/6‑150 (MQ=38)
                                                                                                                                      cgaagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAgg                             >  2:4378286/5‑136 (MQ=255)
                                                                                                                                       gaagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCa              >  2:3623266/4‑150 (MQ=255)
                                                                                                                                       aaCGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCa              >  1:2878254/1‑150 (MQ=255)
                                                                                                                                        aagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCaa             >  2:5122370/3‑150 (MQ=255)
                                                                                                                                          gAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATc           >  2:3267678/1‑150 (MQ=255)
                                                                                                                                          gAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATc           >  2:981553/1‑150 (MQ=255)
                                                                                                                                             ccGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCa        >  2:1920911/1‑150 (MQ=255)
                                                                                                                                                  cACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCAccgcc   >  2:81602/1‑150 (MQ=255)
                                                                                                                                                   aCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCAccgccg  >  2:5326247/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCCGGCCACCAGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCACCGCCG  >  NC_002947/5112305‑5112601

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: