Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,915,748 C→T 6.6% G114S (GGC→AGC)  ccmC ← protoheme IX reservoir complex subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,915,7480CT 6.6% 264.5 / 4.3 90G114S (GGC→AGC) ccmCprotoheme IX reservoir complex subunit
Reads supporting (aligned to +/- strand):  ref base C (40/44);  new base T (3/3);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

CACAGGCCTTGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCC  >  NC_002947/4915601‑4915893
                                                                                                                                                    |                                                                                                                                                  
cacaGGCCTTGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCg                                                                                                                                                   <  3:34079/149‑1 (MQ=255)
cacaGGCCTTGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCg                                                                                                                                                   <  6:136773/149‑1 (MQ=255)
   aGGCCTTGGCCGTACTATCACGATTACTGATCGCCTGGCCTAGCGCAAATATGCCGAAGTACAGGATCAGCAGGATGAGCTTGTACGTAAGGCGGGCATCGGACACCCACCAGCTGCCCCAGTTCGGCTTGCCCCAGATGGCGCCGGTa                                                                                                                                                <  4:354706/149‑1 (MQ=255)
   aGGCCTTGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGt                                                                                                                                                 >  3:286804/1‑148 (MQ=255)
      ccTTGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACACGAACAGCAGGATGAGCATGTACGTAAGGAGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAAc                                                                                                                                              >  7:304806/1‑148 (MQ=255)
       cTTGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCa                                                                                                                                            <  2:12296/149‑1 (MQ=255)
        ttGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAg                                                                                                                                           >  8:201779/1‑149 (MQ=255)
        ttGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAg                                                                                                                                           >  7:149561/1‑149 (MQ=255)
             cGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCa                                                                                                                                      >  2:360271/1‑149 (MQ=255)
                aCTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCTGGTTACCGGCGCCatct                                                                                                                                   <  3:42245/149‑4 (MQ=255)
                 cTATCACGATTCCTCATCGCCTGGCCCAGCGAAATAATGCCGAAGAACAGCAACAGCAGTATCAGAATGCCCGTACGGCGGCCACAGCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCAcggc                                                                                                                                  <  1:347983/149‑1 (MQ=255)
                  tATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCAcggcg                                                                                                                                >  1:293183/1‑149 (MQ=255)
                            aCTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAg                                                                                                                       <  7:1182/148‑1 (MQ=255)
                              tGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGgcg                                                                                                                    >  3:322695/1‑149 (MQ=255)
                              tGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGgcg                                                                                                                    >  7:254873/1‑149 (MQ=255)
                              tGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGgcg                                                                                                                    >  3:428530/1‑149 (MQ=255)
                               gATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGgcgc                                                                                                                   >  3:136223/1‑149 (MQ=255)
                                aTCGCTTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCCTCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGCTGGCGCCGGTAACCCGCGCCCCGGCGGTCATCCAGGCGcc                                                                                                                  >  2:399165/1‑149 (MQ=255)
                                 tCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCCGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCCCCACGGCGGACATCCAgccgcc                                                                                                                  >  5:402820/1‑148 (MQ=255)
                                    ccTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCTGGTTACCGGCGcc                                                                                                                                       <  1:351683/125‑1 (MQ=255)
                                    ccTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCTGGTTACCGGCGcc                                                                                                                                       >  2:351683/1‑125 (MQ=255)
                                     cTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATgg                                                                                                             >  4:236085/1‑149 (MQ=255)
                                        gCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGgc                                                                                                           <  1:360271/148‑1 (MQ=255)
                                         cccAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGgcgc                                                                                                         <  4:136223/149‑1 (MQ=255)
                                         cccAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGgcgc                                                                                                         >  8:300483/1‑149 (MQ=255)
                                          ccAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGcgcg                                                                                                        >  5:260236/1‑149 (MQ=255)
                                                  ataatGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTg                                                                                                <  2:176194/149‑1 (MQ=255)
                                                     atGCCGCAGTACAGGAACAGCAGGATGAGCCTGGACGGACGGCGGGCATCCCACCCCCACCAGCTGCCCCAGGGCGGCTTGCCCCAGAGGGCGCCGGTAACCAGCGCCCCGGCGGGCATCCAGGCGCCAATGGGCGCGGCACACTGCAg                                                                                             <  1:399165/149‑1 (MQ=255)
                                                      tGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAgg                                                                                            <  4:344635/149‑1 (MQ=255)
                                                        ccGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAggg                                                                                           <  1:119927/149‑1 (MQ=255)
                                                        ccGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGc                                                                                          >  4:141625/1‑149 (MQ=255)
                                                         cGAAGTACAGGAACAGCACGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATTGGCGCGGCACACTGCGGGGc                                                                                          <  6:444553/148‑1 (MQ=255)
                                                          gAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGa                                                                                        <  3:141625/149‑1 (MQ=255)
                                                                  ggAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGcc                                                                                <  4:428530/149‑1 (MQ=255)
                                                                   gAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGc                                                                                 >  6:60827/1‑147 (MQ=255)
                                                                   gAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGc                                                                                 >  4:50656/1‑147 (MQ=255)
                                                                   gAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGc                                                                                 >  8:409390/1‑147 (MQ=255)
                                                                      cagcagGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGtt                                                                            >  4:361139/1‑149 (MQ=255)
                                                                       agcagGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGttt                                                                           >  1:388797/1‑149 (MQ=255)
                                                                         cagGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCa                                                                         >  5:174696/1‑149 (MQ=255)
                                                                                aGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCca                                                                  <  4:322695/149‑1 (MQ=255)
                                                                                     ggACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCacac                                                               >  5:440915/1‑147 (MQ=255)
                                                                                         gTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTc                                                                    <  5:15362/138‑1 (MQ=255)
                                                                                         gTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTc                                                                    >  6:15362/1‑138 (MQ=255)
                                                                                            aGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCccac                                                      <  2:388797/149‑1 (MQ=255)
                                                                                            aGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCccac                                                      >  2:289733/1‑149 (MQ=255)
                                                                                              gCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCccacca                                                    >  4:195567/1‑149 (MQ=255)
                                                                                                   cATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCg                                               <  3:50656/149‑1 (MQ=255)
                                                                                                      cccACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGc                                             >  5:344282/1‑148 (MQ=255)
                                                                                                        cacaACCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATGAAGTCGCCAATGGGCTCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCAAACCAGCCACACCACCCCGGCCAc                                          >  7:62950/1‑149 (MQ=255)
                                                                                                          cacccaccAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCg                                        <  6:174696/149‑1 (MQ=255)
                                                                                                            cccaccAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGc                                       >  2:5182/1‑148 (MQ=255)
                                                                                                             ccaccaGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCca                                     <  1:289733/149‑1 (MQ=255)
                                                                                                             ccaccaGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCca                                     <  2:293183/149‑1 (MQ=255)
                                                                                                              caccaGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCCCCCCGGCCACCGCcag                                    >  4:283333/1‑149 (MQ=255)
                                                                                                               accaGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCcagc                                   <  1:145345/149‑1 (MQ=255)
                                                                                                                 caGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACCCCAGCCCCACCACCCCGGCCACCGCcag                                    >  1:332042/1‑147 (MQ=255)
                                                                                                                  aGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATACAGGCGCCAATGGGCGCGGCACACTGCAGGGCGAAGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCCGCCACCGCCagcagc                                >  6:310004/1‑149 (MQ=255)
                                                                                                                   gCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGTCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCGCGGCCACCGCCAgcagca                               <  8:304806/149‑1 (MQ=255)
                                                                                                                      gCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGcaca                             >  3:264929/1‑148 (MQ=255)
                                                                                                                       ccccAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCacac                                                               >  8:138526/1‑113 (MQ=255)
                                                                                                                       ccccAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCacac                                                               <  7:138526/113‑1 (MQ=255)
                                                                                                                       ccccAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGcaca                             <  8:149561/147‑1 (MQ=255)
                                                                                                                       ccccAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACAt                            <  3:236085/148‑1 (MQ=255)
                                                                                                                        cccAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAg                          <  1:40372/149‑1 (MQ=255)
                                                                                                                           aGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAg                       <  4:286804/149‑1 (MQ=255)
                                                                                                                           aGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAg                       >  8:346640/1‑149 (MQ=255)
                                                                                                                           aGGGCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAg                       <  4:149372/149‑1 (MQ=255)
                                                                                                                             gTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGa                     <  6:260236/149‑1 (MQ=255)
                                                                                                                             gTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGa                     <  6:440915/149‑1 (MQ=255)
                                                                                                                                ggCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACt                   <  3:195567/148‑1 (MQ=255)
                                                                                                                                 gCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTgc                 <  3:283333/149‑1 (MQ=255)
                                                                                                                                   ttGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACAt                            <  6:85414/136‑1 (MQ=255)
                                                                                                                                   ttGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACAt                            >  5:85414/1‑136 (MQ=255)
                                                                                                                                     gCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTgcgc               >  5:188075/1‑147 (MQ=255)
                                                                                                                                     gCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTgcgc               <  7:300483/147‑1 (MQ=255)
                                                                                                                                     gCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGcc              <  7:409390/148‑1 (MQ=255)
                                                                                                                                     cccccATATGCGGCCGGTAACCAGCGCCACGGCGGTCATCCACGCGCCAATGGGCGCGGCACACTGCACGGCGACGGCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGcc              <  8:191235/147‑1 (MQ=255)
                                                                                                                                       cccAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAg            >  2:314777/1‑148 (MQ=255)
                                                                                                                                       cccAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGc           <  6:188075/149‑1 (MQ=255)
                                                                                                                                        cgccgtGGCGCTGGTTACCGGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCa          >  1:124133/6‑149 (MQ=255)
                                                                                                                                          ccgtGGCGCTGGTTACCGGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATc        >  8:377562/4‑149 (MQ=255)
                                                                                                                                          ccgtGGCGCTGGTTACCGGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATc        <  2:124133/146‑1 (MQ=255)
                                                                                                                                          aGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATc        <  6:344282/149‑1 (MQ=255)
                                                                                                                                               gcgcCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATcgccgc   <  2:332042/149‑1 (MQ=255)
                                                                                                                                                cgcCTGTAACCAGCGCCACGTCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCACACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCgccgcc  <  8:62950/149‑1 (MQ=255)
                                                                                                                                                cgcCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAg                          <  4:40259/125‑1 (MQ=255)
                                                                                                                                                cgcCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAg                          >  3:40259/1‑125 (MQ=255)
                                                                                                                                                 gcCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTgcg                >  6:350591/1‑134 (MQ=255)
                                                                                                                                                 gcCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTgcg                <  5:350591/134‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CACAGGCCTTGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCC  >  NC_002947/4915601‑4915893

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