Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,112,451 C→G 5.3% V75V (GTG→GTC) ‡ 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,4510CG 5.3% 562.9 / 15.4 169V75V (GTG→GTC) ‡aroP‑IIaromatic amino acid transport protein
Reads supporting (aligned to +/- strand):  ref base C (89/71);  new base G (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.94e-01

Only 100 of 172 total aligned reads displayed.

TCCGGCCACCAGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCACCGCC  >  NC_002947/5112305‑5112600
                                                                                                                                                  |                                                                                                                                                     
tcCGGCCACCAGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCAcc                                                                                                                                                    >  1:2388047/1‑150 (MQ=255)
       accaGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGTGGCCGGctcgt                                                                                                                                             <  2:2056443/150‑3 (MQ=255)
         caGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGACCACCGGCTCTTcg                                                                                                                                           >  1:5049547/1‑150 (MQ=255)
          agtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcga                                                                                                                                          <  2:2664590/150‑1 (MQ=255)
          agtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAAGAGCCGGTGGCCGGctcgttcag                                                                                                                                         <  1:1450446/150‑7 (MQ=255)
           gtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcgac                                                                                                                                         <  1:2353175/150‑1 (MQ=255)
           gtagtGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcgac                                                                                                                                         <  1:3847638/150‑1 (MQ=255)
                gAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCa                                                                                                                                    >  2:3756901/1‑150 (MQ=255)
                gAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCa                                                                                                                                    <  1:4763645/150‑1 (MQ=255)
                 aaCATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCAt                                                                                                                                   >  2:4497072/1‑150 (MQ=255)
                 aaCATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCAt                                                                                                                                   >  2:1220677/1‑150 (MQ=255)
                     tACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCa                                                                                                                               <  1:1368807/150‑1 (MQ=255)
                     tACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCa                                                                                                                               <  1:1285361/150‑1 (MQ=255)
                       cTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCAcc                                                                                                                             >  1:4568279/1‑150 (MQ=255)
                             gACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTg                                                                                                                       >  2:3480904/1‑150 (MQ=255)
                              aCCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTgg                                                                                                                      >  2:1665967/1‑150 (MQ=255)
                              aCCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTgg                                                                                                                      <  1:4356627/150‑1 (MQ=255)
                               ccgccgAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGgc                                                                                                                     <  2:4068934/150‑1 (MQ=255)
                                 gccgAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGgcgc                                                                                                                   <  2:1008298/150‑1 (MQ=255)
                                    gAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATg                                                                                                                <  1:241792/150‑1 (MQ=255)
                                          tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCTGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                          >  1:1149231/1‑150 (MQ=255)
                                          tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                          >  1:4679590/1‑150 (MQ=255)
                                          tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                          >  1:1492389/1‑150 (MQ=255)
                                          tCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATg                                                                                                          >  1:3384562/1‑150 (MQ=255)
                                             gACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAg                                                                                                       <  1:3750431/150‑1 (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:2555670/150‑1 (MQ=255)
                                                  ccccACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCCTGAAGGCGAt                                                                                                  >  2:4110079/1‑150 (MQ=255)
                                                  ccccACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGAt                                                                                                  <  1:261894/150‑1 (MQ=255)
                                                   cccACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATg                                                                                                 >  2:3915751/1‑150 (MQ=255)
                                                   cccACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATg                                                                                                 >  2:1736113/1‑150 (MQ=255)
                                                      accaGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAg                                                                                              >  2:5127193/1‑150 (MQ=255)
                                                      accaGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAg                                                                                              >  2:5128064/1‑150 (MQ=255)
                                                         aGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTcgac                                                                                                                                         >  2:2363729/1‑104 (MQ=255)
                                                           gATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCa                                                                                         >  1:3432061/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)
                                                                aCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAACTGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGACGCCGCAGATcg                                                                                    >  2:2982557/1‑150 (MQ=255)
                                                                aCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCAc                                                                                                                                                     <  1:1370857/85‑1 (MQ=255)
                                                                aCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATc                                                                                                             >  1:5425199/1‑125 (MQ=255)
                                                                aCAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcg                                                                                    >  1:2024595/1‑150 (MQ=255)
                                                                 cAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcgc                                                                                   >  1:272546/1‑150 (MQ=255)
                                                                 cAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATcgc                                                                                   >  2:239164/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)
                                                                       ccAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGcc                                                                             >  2:532357/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                                                                                                                               >  2:361437/1‑91 (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)
                                                                                      cGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACgg                                                              >  2:4967256/1‑150 (MQ=255)
                                                                                       gACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACggg                                                             >  2:4913638/1‑150 (MQ=255)
                                                                                            gAAGCTGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCcggc                                                        >  2:1327641/1‑150 (MQ=255)
                                                                                            gAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCcggc                                                        >  2:4795455/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)
                                                                                                    cGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTtca                                                <  2:1411672/150‑1 (MQ=255)
                                                                                                    cGAACCCACCCCAGTAGGTGGGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTtca                                                <  1:5137792/150‑1 (MQ=255)
                                                                                                      cgcccacccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTtcatc                                              >  1:2491796/3‑150 (MQ=255)
                                                                                                         ccacccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACg                                           >  1:5234232/1‑150 (MQ=255)
                                                                                                         ccacccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACg                                           >  2:3237843/1‑150 (MQ=255)
                                                                                                           acccCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGcc                                         >  2:5156113/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)
                                                                                                                     gtgtGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTgccgagccg                               >  1:2985729/1‑150 (MQ=255)
                                                                                                                      tgtgGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTgccgagccga                              <  1:1736113/150‑1 (MQ=255)
                                                                                                                        tgGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAgg                            <  1:544302/150‑1 (MQ=255)
                                                                                                                          ggCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGaa                          >  2:4343890/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                                                                    >  1:5141224/1‑99 (MQ=255)
                                                                                                                                   gCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGatcat                                                                    <  2:5141224/99‑1 (MQ=255)
                                                                                                                                    cTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTa                <  2:638023/150‑1 (MQ=255)
                                                                                                                                      cgaagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAgg                            <  1:4378286/132‑1 (MQ=255)
                                                                                                                                      cgaagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAgg                            >  2:4378286/5‑136 (MQ=255)
                                                                                                                                       aaCGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCa             >  2:4607737/1‑150 (MQ=255)
                                                                                                                                       aaCGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCa             >  1:2878439/1‑150 (MQ=255)
                                                                                                                                       aaCGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCa             >  1:2878254/1‑150 (MQ=255)
                                                                                                                                       aaCGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGACCGGTACCa             >  2:4517984/1‑150 (MQ=255)
                                                                                                                                        aagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCaa            >  1:255970/3‑150 (MQ=255)
                                                                                                                                        aagAGCCGGTGGCCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCaa            >  2:5122370/3‑150 (MQ=255)
                                                                                                                                          gAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATc          >  2:981553/1‑150 (MQ=255)
                                                                                                                                             ccGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCa       >  2:1920911/1‑150 (MQ=255)
                                                                                                                                               ggCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCAcc     >  1:4730599/1‑150 (MQ=255)
                                                                                                                                                  cACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCAccgcc  >  2:81602/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCCGGCCACCAGTAGTGAACATACTTGCCGACCGCCGAAAGCTCCGACATCCCCACCAGGATATACAGCACCCAGCAGTTCCAGCCCGACAGGAAGCCGGCGAACCCACCCCAGTAGGTGTGGGCGAAATGGCTGAACGAGCCGGCCACCGGCTCTTCGACGATCATTTCACCCAGCTGGCGCATGATCATGAAGGCGATGAAGCCGCAGATCGCGTAGCCAAGGATCATCGACGGGCCGGCGGATTTCATCACGCCTGCCGAGCCGAGGAACAGGCCGGTACCAATCGCACCGCC  >  NC_002947/5112305‑5112600

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: