Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,331,401 C→A 6.2% G54C (GGC→TGC) ‡ 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,4010CA 6.2% 257.1 / 4.5 80G54C (GGC→TGC) ‡PP_4692aspartate/tyrosine/aromatic aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (37/38);  new base A (3/2);  total (40/40)
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.46e-01

CGGAGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCACTGTA  >  NC_002947/5331254‑5331548
                                                                                                                                                      |                                                                                                                                                   
cGGAGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTggccg                                                                                                                                                    <  1:275781/149‑1 (MQ=255)
   aGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTggccggcc                                                                                                                                                 <  8:368706/149‑1 (MQ=255)
   aGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTggccggcc                                                                                                                                                 <  2:181338/149‑1 (MQ=255)
   aGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTggccggcc                                                                                                                                                 >  5:157246/1‑149 (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                                                                                                                                     <  4:286359/149‑1 (MQ=255)
               taataaGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGt                                                                                                                                     <  6:265223/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)
                     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                                                                                                                         >  7:212194/1‑148 (MQ=255)
                            ttCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCgtggtg                                                                                                                         <  8:212194/148‑1 (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                                                                                >  8:308481/1‑136 (MQ=255)
                                                                                  cACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGccc                                                                                <  7:308481/136‑1 (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                                                        >  7:212827/1‑149 (MQ=255)
                                                                                               ggCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCg                                                        <  5:396243/149‑1 (MQ=255)
                                                                                               ggCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCg                                                        >  6:396243/1‑149 (MQ=255)
                                                                                                gCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCgg                                                       <  6:352405/149‑1 (MQ=255)
                                                                                                gCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCgg                                                       >  2:193780/1‑149 (MQ=255)
                                                                                                                    tAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCgg                                                       <  8:381626/129‑1 (MQ=255)
                                                                                                                    tAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCgg                                                       >  7:381626/1‑129 (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                          >  4:396145/1‑149 (MQ=255)
                                                                                                                             tgCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTc                          <  6:115686/149‑1 (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:255033/139‑1 (MQ=255)
                                                                                                                                   gCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAg                              <  5:260838/139‑1 (MQ=255)
                                                                                                                                   gCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAg                              >  6:260838/1‑139 (MQ=255)
                                                                                                                                   gCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAg                              >  6:255033/1‑139 (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      <  3:346082/149‑1 (MQ=255)
                                                                                                                                                 cTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCAc      >  7:181359/1‑149 (MQ=255)
                                                                                                                                                 cTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCAc      >  4:346082/1‑149 (MQ=255)
                                                                                                                                                  tGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCACt     >  6:361304/1‑149 (MQ=255)
                                                                                                                                                    gccggccgCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCACTg    >  1:211025/1‑148 (MQ=255)
                                                                                                                                                     ccggccgCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCACTGTa  <  8:190483/149‑1 (MQ=255)
                                                                                                                                                      cggccgCAACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGtt                           <  7:318995/123‑1 (MQ=255)
                                                                                                                                                      cggccgCAACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGtt                           >  8:318995/1‑123 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CGGAGCCACCGGGGGTAATAAGGATGCGTTCCGGGTCGATGTTCAGGCCATAGCGCTGCCCATAAAAACCGGCAATTGCCTCACGCAGTGCCGGCAGACCGCGTGCGGCGGTGTAGCGAGTGTGCCCTGCGGCCAGCGCGGCCTGGCCGGCCGCGACGATGGGTGCGGCCGTGGTGAAGTCCGGCTCGCCGATTTCCAGGTGGATCACGTCGTGCCCGGCCGCCTGCAGCTCGTTGGCCCGGGCCAGCAGCGCCATGACGTGGAAAGGTTCGATGGCGCGGCTGCGCGCACTGTA  >  NC_002947/5331254‑5331548

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