Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,915,756 A→G 6.9% L111P (CTG→CCG)  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,7560AG 6.9% 247.6 / 6.7 88L111P (CTG→CCG) ccmCprotoheme IX reservoir complex subunit
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base G (3/3);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

CTTGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCACATGGA  >  NC_002947/4915608‑4915903
                                                                                                                                                     |                                                                                                                                                    
cTTGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCa                                                                                                                                                      <  2:12296/149‑1 (MQ=255)
 ttGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAg                                                                                                                                                     >  7:149561/1‑149 (MQ=255)
 ttGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAg                                                                                                                                                     >  8:201779/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                                                                                                                              >  7:254873/1‑149 (MQ=255)
                       tGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGgcg                                                                                                                              >  3:322695/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                                                                              >  6:15362/1‑138 (MQ=255)
                                                                                  gTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTc                                                                              <  5:15362/138‑1 (MQ=255)
                                                                                     aGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCccac                                                                >  2:289733/1‑149 (MQ=255)
                                                                                     aGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCccac                                                                <  2:388797/149‑1 (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                                 >  8:346640/1‑149 (MQ=255)
                                                                                                                    aGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAg                                 <  4:286804/149‑1 (MQ=255)
                                                                                                                    aGGGCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAg                                 <  4:149372/149‑1 (MQ=255)
                                                                                                                      gTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGa                               <  6:440915/149‑1 (MQ=255)
                                                                                                                      gTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGa                               <  6:260236/149‑1 (MQ=255)
                                                                                                                         ggCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACt                             <  3:195567/148‑1 (MQ=255)
                                                                                                                          gCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTgc                           <  3:283333/149‑1 (MQ=255)
                                                                                                                            ttGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACAt                                      >  5:85414/1‑136 (MQ=255)
                                                                                                                            ttGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACAt                                      <  6:85414/136‑1 (MQ=255)
                                                                                                                              gCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTgcgc                         <  7:300483/147‑1 (MQ=255)
                                                                                                                              gCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTgcgc                         >  5:188075/1‑147 (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)
                                                                                                                                              ggTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGca        <  4:264929/149‑1 (MQ=255)
                                                                                                                                                   ccAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCACATgg   >  5:366206/1‑149 (MQ=255)
                                                                                                                                                    cAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCACATGGa  >  3:284409/1‑149 (MQ=255)
                                                                                                                                                    cAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCACATGGa  >  3:201092/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTGGCCGCACTATCACGATTACTGATCGCCTGGCCCAGCGCAATAATGCCGAAGTACAGGAACAGCAGGATGAGCATGGACGTAAGGCGGGCATCCCACACCCACCAGCTGCCCCAGGTCGGCTTGCCCCAGATGGCGCCGGTAACCAGCGCCACGGCGGTCATCCAGGCGCCAATGGGCGCGGCACACTGCAGGGCGACGTCGGCCAGTTTCATCTTCCACACCAGCCCCACCACCCCGGCCACCGCCAGCAGCACATAGCAGGACTGCGCCAGCATCGCCGCCGGCACATGGA  >  NC_002947/4915608‑4915903

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: