Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,828,109 A→T 5.2% S14C (AGC→TGC)  yggW → coproporphyrinogen/heterocyclic compound oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,828,1090AT 5.2% 557.2 / 3.8 169S14C (AGC→TGC) yggWcoproporphyrinogen/heterocyclic compound oxidase
Reads supporting (aligned to +/- strand):  ref base A (30/128);  new base T (7/2);  total (39/130)
Fisher's exact test for biased strand distribution p-value = 4.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

TGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAGCTG  >  NC_002947/5827960‑5828258
                                                                                                                                                      |                                                                                                                                                     
tGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCCGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTa                                                                                                                                                       <  1:4302620/150‑1 (MQ=255)
tGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTa                                                                                                                                                       <  2:3351808/150‑1 (MQ=255)
tGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTa                                                                                                                                                       <  1:2997508/150‑1 (MQ=255)
 gggTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTag                                                                                                                                                      <  1:5792754/150‑1 (MQ=255)
 gggTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTag                                                                                                                                                      <  1:2425678/150‑1 (MQ=255)
  ggTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTagc                                                                                                                                                     <  2:6612017/150‑1 (MQ=255)
    tACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCGTTagcag                                                                                                                                                   >  2:6683491/1‑150 (MQ=255)
       cGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCCGCAGACCAGACCCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCt                                                                                                                                                <  2:6451890/150‑1 (MQ=255)
       cGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCt                                                                                                                                                <  1:5558729/150‑1 (MQ=255)
       cGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCt                                                                                                                                                <  1:5045992/150‑1 (MQ=255)
        ggAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGCGGCGTGTGCAGccag                                                                                                                                               >  2:3309425/1‑148 (MQ=255)
        ggAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGAGGCTTTAGCAGCCTg                                                                                                                                               <  2:4914774/150‑1 (MQ=255)
        ggAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCCGGCGGTTGCAGccgg                                                                                                                                               >  2:2744572/1‑148 (MQ=255)
        ggAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTg                                                                                                                                               <  1:2104140/150‑1 (MQ=255)
        ggAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTg                                                                                                                                               <  2:4148373/150‑1 (MQ=255)
         gAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCCGGCGGTAGCAGCCGgc                                                                                                                                              >  2:3183426/1‑150 (MQ=255)
          aGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTgcc                                                                                                                                             <  2:2181752/150‑1 (MQ=255)
           gcgcAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGTGGCAGGCTTTAGCAGCCTgcca                                                                                                                                            <  2:4981351/150‑2 (MQ=255)
               aaCTGCTCCGGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCcgccgc                                                                                                                                        <  2:2389919/150‑1 (MQ=255)
               aaCTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGCGGCGTTTGCAGCCCGCccccgc                                                                                                                                        >  2:4399391/1‑150 (MQ=255)
               aaCTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCATTAGCAGCCTgccgcccc                                                                                                                                        >  2:3404753/1‑148 (MQ=255)
                 cTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTg                                                                                                                                      <  2:6817410/150‑1 (MQ=255)
                 cTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTg                                                                                                                                      <  2:2439024/150‑1 (MQ=255)
                   gCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCGCCAGCGGCGGCAGGCTGCTAAAGCCTgccgcgccggg                                                                                                                                     <  2:773647/150‑7 (MQ=37)
                    cTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCCGGCGGTTGCCGCCAGCCGCCCCAGgcg                                                                                                                                   >  1:4505640/1‑150 (MQ=255)
                    cTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGgcg                                                                                                                                   <  2:4777302/150‑1 (MQ=255)
                     tCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCCGGGGTTGGCCGCCGGCCCCCCCCGgcgc                                                                                                                                  >  2:1832671/1‑150 (MQ=255)
                      ccAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCt                                                                                                                                 <  1:6742090/150‑1 (MQ=255)
                        aGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGt                                                                                                                               <  1:1080094/150‑1 (MQ=255)
                         gcgcCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTa                                                                                                                              <  1:4146730/150‑1 (MQ=255)
                          cgcCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTAc                                                                                                                             <  2:1351596/150‑1 (MQ=255)
                            ccGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACAt                                                                                                                           <  1:2166462/150‑1 (MQ=255)
                             cGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATc                                                                                                                          <  2:6141010/150‑1 (MQ=255)
                              gATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCCGGGGTTTTCAGCCAGCCCCCCCAGGCGCGGTaaaaac                                                                                                                         >  2:1635545/1‑145 (MQ=255)
                               ctctcGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCAGGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCca                                                                                                                        <  2:4357991/149‑1 (MQ=255)
                                tCTAGCCCCAGCAGACAAGACCCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCcac                                                                                                                       <  1:5971841/150‑1 (MQ=255)
                                     ccccAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTcc                                                                                                                  <  2:5482149/150‑1 (MQ=255)
                                     ccccAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTcc                                                                                                                  <  2:5483033/150‑1 (MQ=255)
                                        cagcagACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCgtg                                                                                                               <  1:6261591/150‑1 (MQ=255)
                                        cagcagACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCgtg                                                                                                               <  2:2541264/150‑1 (MQ=255)
                                        cagcagACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCgtg                                                                                                               <  2:2110256/150‑1 (MQ=255)
                                        cagcagACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCgtg                                                                                                               <  1:986201/150‑1 (MQ=255)
                                          gcagACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCgtggt                                                                                                             <  2:394066/150‑1 (MQ=255)
                                           cagACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCGTTAGCACCCGGCCGCCCCAGGCGCTGTACAACCCCCTACCggggtg                                                                                                            >  2:6205950/1‑150 (MQ=255)
                                              aCAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGt                                                                                                         <  1:5913598/150‑1 (MQ=255)
                                                  gAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCCGGCGGTGGCAGCCCGCCCCCGCCGGCGCGGGACCACCACAATGCGGGGGGCGGACGc                                                                                                     >  2:5078757/1‑150 (MQ=255)
                                                         cTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGcc                                                                                              <  1:6845806/150‑1 (MQ=255)
                                                         cTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGcc                                                                                              <  1:5131497/150‑1 (MQ=255)
                                                           cAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCt                                                                                            <  1:2730739/150‑1 (MQ=255)
                                                              ccACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTAt                                                                                         <  2:884378/150‑1 (MQ=255)
                                                                   gcgcCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGa                                                                                    <  1:1855362/150‑1 (MQ=255)
                                                                       ccGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTc                                                                                <  2:5700689/150‑1 (MQ=255)
                                                                       ccGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTc                                                                                <  1:670556/150‑1 (MQ=255)
                                                                       ccGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTc                                                                                <  2:1759922/150‑1 (MQ=255)
                                                                        cGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCa                                                                               <  1:1247540/150‑1 (MQ=255)
                                                                         gTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCaa                                                                              <  1:5383167/150‑1 (MQ=255)
                                                                               tGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCa                                                                        <  1:5224720/150‑1 (MQ=255)
                                                                                   ccTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCa                                                                    <  2:5418780/150‑1 (MQ=255)
                                                                                              aaCGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGa                                                         <  2:2516925/150‑1 (MQ=255)
                                                                                               aCGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGaa                                                        <  2:2777987/150‑1 (MQ=255)
                                                                                                   cTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTAc                                                    <  2:1326495/150‑1 (MQ=255)
                                                                                                      ggCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCgg                                                 <  2:2625637/150‑1 (MQ=255)
                                                                                                       gCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGa                                                <  2:1045590/150‑1 (MQ=255)
                                                                                                        ccTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGaa                                               <  1:1447512/150‑1 (MQ=255)
                                                                                                          tctcATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGa                                             <  1:6833929/147‑1 (MQ=255)
                                                                                                           ggCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGAt                                            <  2:3558771/150‑1 (MQ=255)
                                                                                                              aTGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTACCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGcc                                         >  1:1127752/1‑150 (MQ=255)
                                                                                                                gATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTa                                       <  2:1245447/150‑1 (MQ=255)
                                                                                                                   cGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGt                                    <  2:3877838/150‑1 (MQ=255)
                                                                                                                     aaaTGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTca                                  <  1:6573242/150‑1 (MQ=255)
                                                                                                                       aTGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTcaca                                <  1:3145906/150‑1 (MQ=255)
                                                                                                                       aTGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTcaca                                <  2:5291885/150‑1 (MQ=255)
                                                                                                                              cGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGc                         <  1:5769008/150‑1 (MQ=255)
                                                                                                                                   cccGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTcgac                    <  1:857459/150‑1 (MQ=255)
                                                                                                                                      ggCGCGGCCGGCTTTAGCAGCCTGCCCCCGCTGGCGCTGTACATCCACAATCCGTGGTGCGTACGCAAATGCCCTTATTTCGACTTCAACTACCACCCAGCCGGGCCAGAAAAACCCGCACATGCCAACGTCACAGCCCAGCTCgacgac                 >  2:5891160/1‑150 (MQ=255)
                                                                                                                                       gcgcGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACCTTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGATGATGCCTACGTCACAGCCCTGCTCGACGAcc                <  1:2254987/150‑1 (MQ=255)
                                                                                                                                       gcgcGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGAcc                <  1:324704/150‑1 (MQ=255)
                                                                                                                                       gcgcGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGAcc                <  1:4436050/150‑1 (MQ=255)
                                                                                                                                          gcgGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACAGTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTc              <  1:3560987/150‑1 (MQ=255)
                                                                                                                                            ggCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGa            <  1:1326465/150‑1 (MQ=255)
                                                                                                                                            ggCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGa            <  1:3183702/150‑1 (MQ=255)
                                                                                                                                            ggCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGa            >  1:6789542/1‑150 (MQ=255)
                                                                                                                                             gCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGAc           <  1:6397826/150‑1 (MQ=255)
                                                                                                                                             gCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGAc           <  2:6313080/150‑1 (MQ=255)
                                                                                                                                              cAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGAcc          >  2:297819/1‑150 (MQ=255)
                                                                                                                                              cAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGAcc          <  2:1109519/150‑1 (MQ=255)
                                                                                                                                              cAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGAcc          >  1:3296203/1‑150 (MQ=255)
                                                                                                                                              cAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGAcc          >  2:5676801/1‑150 (MQ=255)
                                                                                                                                               aGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCCCAGTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCa         <  2:692440/150‑1 (MQ=255)
                                                                                                                                               aGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCa         <  2:2958789/150‑1 (MQ=255)
                                                                                                                                               aGGCTTTAGCAGCCTGCCGCCGCGGGCGCTGTACATCCACAGTCCGCGGTGCGTACGCAAATGCCCTTAGTGCGACTTCAACTCCCCCGCAGCCGGGCCTGAACGACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCa         <  2:6537945/150‑1 (MQ=255)
                                                                                                                                                ggCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAg        <  2:4931349/150‑1 (MQ=255)
                                                                                                                                                 gCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAgg       >  1:6421442/1‑150 (MQ=255)
                                                                                                                                                  cTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGa      <  2:3990090/150‑1 (MQ=255)
                                                                                                                                                   tttAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAg     >  1:6756583/1‑150 (MQ=255)
                                                                                                                                                   tttAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAg     >  1:3999369/1‑150 (MQ=255)
                                                                                                                                                   tttAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAg     >  1:4911742/1‑150 (MQ=255)
                                                                                                                                                   tatAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCACTGCTCGACTACCTCGACCAGGAg     >  2:4634484/3‑150 (MQ=255)
                                                                                                                                                      agcagcCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgctg  >  2:5543236/1‑150 (MQ=255)
                                                                                                                                                      agcagcCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgctg  >  2:171470/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAGCTG  >  NC_002947/5827960‑5828258

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