Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,331,399 G→C 7.2% G54G (GGC→GGG) ‡ PP_4692 ← aspartate/tyrosine/aromatic aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,331,3990GC 7.2% 256.9 / 6.8 80G54G (GGC→GGG) ‡PP_4692aspartate/tyrosine/aromatic aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base C (4/2);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CGGAGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCACTG  >  NC_002947/5331254‑5331546
                                                                                                                                                    |                                                                                                                                                   
cGGAGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTggccg                                                                                                                                                  <  1:275781/149‑1 (MQ=255)
   aGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTggccggcc                                                                                                                                               >  5:157246/1‑149 (MQ=255)
   aGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTggccggcc                                                                                                                                               <  8:368706/149‑1 (MQ=255)
   aGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTggccggcc                                                                                                                                               <  2:181338/149‑1 (MQ=255)
    gCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGgccggccg                                                                                                                                              >  5:107451/1‑149 (MQ=255)
     ccACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCcgc                                                                                                                                             >  4:323434/1‑149 (MQ=255)
       aCCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGcga                                                                                                                                           >  5:8397/1‑149 (MQ=255)
       aCCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGcga                                                                                                                                           <  8:235972/149‑1 (MQ=255)
        ccGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCGGCCAGGCcgcgct                                                                                                                                          <  5:280854/149‑3 (MQ=255)
         cGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGcgac                                                                                                                                          >  5:396650/1‑148 (MQ=255)
           ggggTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGAt                                                                                                                                       <  6:4227/149‑1 (MQ=255)
            gggTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATg                                                                                                                                      <  1:212167/149‑1 (MQ=255)
             ggTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATgg                                                                                                                                     <  7:368307/149‑1 (MQ=255)
             ggTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATgg                                                                                                                                     >  1:213302/1‑149 (MQ=255)
             ggTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCAGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATgg                                                                                                                                     <  5:254948/149‑1 (MQ=255)
               taataaGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGt                                                                                                                                   <  6:265223/149‑1 (MQ=255)
               taataaGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGt                                                                                                                                   <  4:286359/149‑1 (MQ=255)
                aataaGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCcgc                                                                                                                                             >  4:178705/1‑138 (MQ=255)
                aataaGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCcgc                                                                                                                                             <  3:178705/138‑1 (MQ=255)
                  taaGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCAAAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGAGGGGTGCg                                                                                                                                >  5:276787/1‑149 (MQ=255)
                   aaGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTgg                                                                                                                                                     <  2:104671/127‑1 (MQ=255)
                   aaGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTgg                                                                                                                                                     >  1:104671/1‑127 (MQ=255)
                     ggATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGcc                                                                                                                             <  8:114753/149‑1 (MQ=255)
                      gATGCGTTACGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGAAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGCCGATGGGTGCGGCCg                                                                                                                            >  6:2321/1‑149 (MQ=255)
                        tGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCgtg                                                                                                                          <  6:107451/149‑1 (MQ=255)
                         gCGTTCGGGGTCGATGTTCAGCCAATAGGGCTGCCAAAAAAACCGGTAAATTCCCCCTACGACCCCCCCCATAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCgt                                                                                                                           <  1:335084/149‑1 (MQ=255)
                          cGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCgtggt                                                                                                                        >  6:210205/1‑149 (MQ=255)
                           gTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGGCCGCGCTggccgcagggc                                                                                                                                            <  8:265127/128‑7 (MQ=255)
                           gTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCgtggt                                                                                                                        >  2:300033/1‑148 (MQ=255)
                           gTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCgtg                                                                                                                          >  1:265361/1‑146 (MQ=255)
                           gTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCgtg                                                                                                                          <  2:265361/146‑1 (MQ=255)
                            ttCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCgtggtg                                                                                                                       <  8:212194/148‑1 (MQ=255)
                            ttCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCgtggtg                                                                                                                       >  7:212194/1‑148 (MQ=255)
                                 ggTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTcc                                                                                                                 <  6:271428/149‑1 (MQ=255)
                                     gATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCt                                                                                                             <  3:323434/149‑1 (MQ=255)
                                         ttCAGGCCATAGCGCTGCCCATAAAAACCGGCAACTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGcc                                                                                                         <  5:39298/149‑1 (MQ=255)
                                               ccATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTc                                                                                                   >  1:246033/1‑149 (MQ=255)
                                                 aTAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATTGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCa                                                                                                 <  2:213302/149‑1 (MQ=255)
                                                      gcTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTgg                                                                                            <  6:157246/149‑1 (MQ=255)
                                                          cccATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATc                                                                                         <  5:328298/148‑1 (MQ=255)
                                                          cccATAAAAACCGGCAATTGACTCACTCACTGCCGGCAGACCGCGTGCGGCGGTGTAGCGACTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGAAGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCa                                                                                        >  4:173721/1‑149 (MQ=255)
                                                                 aaaCCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTg                                                                                 <  7:207842/149‑1 (MQ=255)
                                                                  aaCCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGc                                                                                >  5:352405/1‑149 (MQ=255)
                                                                     cGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCg                                                                             >  4:384312/1‑149 (MQ=255)
                                                                                  cACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGccc                                                                              <  7:308481/136‑1 (MQ=255)
                                                                                  cACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGccc                                                                              >  8:308481/1‑136 (MQ=255)
                                                                                      gCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGGGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCg                                                              <  5:2321/149‑1 (MQ=255)
                                                                                        cAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGt                                                             <  6:396650/149‑1 (MQ=255)
                                                                                         aGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGtt                                                            <  3:384312/149‑1 (MQ=255)
                                                                                          gTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTg                                                           <  3:86393/149‑1 (MQ=255)
                                                                                             ccGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGccc                                                       >  6:197862/1‑149 (MQ=255)
                                                                                               ggCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCg                                                      <  5:396243/149‑1 (MQ=255)
                                                                                               ggCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCg                                                      >  7:212827/1‑149 (MQ=255)
                                                                                               ggCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCg                                                      >  6:396243/1‑149 (MQ=255)
                                                                                                gCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCgg                                                     <  6:352405/149‑1 (MQ=255)
                                                                                                gCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCgg                                                     >  2:193780/1‑149 (MQ=255)
                                                                                                                    tAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCgg                                                     >  7:381626/1‑129 (MQ=255)
                                                                                                                    tAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCgg                                                     <  8:381626/129‑1 (MQ=255)
                                                                                                                      gCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGa                               <  5:129526/149‑1 (MQ=255)
                                                                                                                       cGAGTGTGCCCTGCGGCCAGTGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAg                                                                                                >  4:171834/1‑83 (MQ=255)
                                                                                                                       cGAGTGTGCCCTGCGGCCAGTGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAg                                                                                                <  3:171834/83‑1 (MQ=255)
                                                                                                                       cGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGa                               >  7:35473/1‑148 (MQ=255)
                                                                                                                         aGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCGGGGCCAGCAGCGCCATGACGTGGAAAg                            <  3:173721/149‑1 (MQ=255)
                                                                                                                             tgCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTc                        <  6:115686/149‑1 (MQ=255)
                                                                                                                             tgCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTc                        >  4:396145/1‑149 (MQ=255)
                                                                                                                             tgCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTc                        >  4:190419/1‑149 (MQ=255)
                                                                                                                               cccTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGa                      >  4:243195/1‑149 (MQ=255)
                                                                                                                                ccTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGAt                     >  1:247327/1‑149 (MQ=255)
                                                                                                                                 cTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATg                    <  8:212827/149‑1 (MQ=255)
                                                                                                                                 cTGCGGCCAGCGCGGCCTGCCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATg                    >  6:339245/1‑149 (MQ=255)
                                                                                                                                   gCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAg                            <  5:260838/139‑1 (MQ=255)
                                                                                                                                   gCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAg                            >  6:260838/1‑139 (MQ=255)
                                                                                                                                   gCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAg                            >  6:255033/1‑139 (MQ=255)
                                                                                                                                   gCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAg                            <  5:255033/139‑1 (MQ=255)
                                                                                                                                    cGGCCAGGCCGCGCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGgcg                 >  3:271176/1‑149 (MQ=255)
                                                                                                                                      gCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCAGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGcgcg               >  3:133455/1‑149 (MQ=255)
                                                                                                                                       ccAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCAGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCgg              <  4:133455/149‑1 (MQ=255)
                                                                                                                                       ccAGCGCGGCCTGGCCGGCCGCGACGATGGGTGAGGCCGGGGTGAAGTCCGGCTCGCCGATTTCCAGGGGGATCACGTCGTGCCCGGCCGCCGGCAACTCGGTGGCCCCGGCCAGCAGCGCCATGACGGGGAAAGGTACGAAGGAGCgg              >  3:244529/1‑149 (MQ=255)
                                                                                                                                        cAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGc             <  1:193780/149‑1 (MQ=255)
                                                                                                                                          tcgcgGCCGGCCAGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTgg                                                          <  8:121589/101‑1 (MQ=255)
                                                                                                                                          tcgcgGCCGGCCAGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTgg                                                          >  7:121589/2‑102 (MQ=255)
                                                                                                                                          tcgcgGCCGGCCAGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTg           >  1:322646/2‑149 (MQ=255)
                                                                                                                                          gcgcGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTg           >  7:190483/1‑149 (MQ=255)
                                                                                                                                           cgcgGCCGGCCAGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTgc          >  6:332532/1‑149 (MQ=255)
                                                                                                                                                 cTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCAc    >  4:346082/1‑149 (MQ=255)
                                                                                                                                                 cTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCAc    <  3:346082/149‑1 (MQ=255)
                                                                                                                                                 cTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCAc    >  7:181359/1‑149 (MQ=255)
                                                                                                                                                  tGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCACt   >  6:361304/1‑149 (MQ=255)
                                                                                                                                                    gccggccgCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCACTg  >  1:211025/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGAGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCACTG  >  NC_002947/5331254‑5331546

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: