Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,185,401 G→T 6.4% R353R (CGC→CGA)  aceK ← isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,185,4010GT 6.4% 231.6 / 6.0 78R353R (CGC→CGA) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base T (2/3);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01

CGGCGCCACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGGTTGATCA  >  NC_002947/5185257‑5185549
                                                                                                                                                |                                                                                                                                                    
cGGCGCCACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGAt                                                                                                                                                  <  2:232826/149‑1 (MQ=255)
     ccACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACg                                                                                                                                             >  1:143947/1‑149 (MQ=255)
       aCCTCCAGCAACTCGGCCATGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACGGt                                                                                                                                           <  7:85161/149‑1 (MQ=255)
       aCCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGt                                                                                                                                           >  2:141571/1‑149 (MQ=255)
       aCCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACGGt                                                                                                                                           <  6:298067/149‑1 (MQ=255)
        ccTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGtt                                                                                                                                          >  4:367845/1‑149 (MQ=255)
          tCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGtttt                                                                                                                                        >  8:286930/1‑149 (MQ=255)
           ccAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTc                                                                                                                                       >  5:58039/1‑149 (MQ=255)
            cAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGATAATCGGCATACTCCTGGGTGTCGGCAATACGCCCCACCCGGTCGACAATCTTCACCAGCCGATACTTGTCTATAAACGTGGCGCGATCGACGGTTTTc                                                                                                                                       >  5:319379/1‑148 (MQ=255)
            cAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTc                                                                                                                                       <  1:141571/148‑1 (MQ=255)
            cAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGGCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTc                                                                                                                                       <  5:162111/148‑1 (MQ=255)
               cAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATg                                                                                                                                   <  6:343506/149‑1 (MQ=255)
               cAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATg                                                                                                                                   <  6:137557/149‑1 (MQ=255)
                aaCTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACgg                                                                                                                                            >  7:298303/1‑139 (MQ=255)
                aaCTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACgg                                                                                                                                            <  8:298303/139‑1 (MQ=255)
                    cGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGc                                                                                                                                 >  2:58916/1‑146 (MQ=255)
                        cAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGc                                                                                                                          <  3:44818/149‑1 (MQ=255)
                        cAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGc                                                                                                                          >  8:267177/1‑149 (MQ=255)
                         aGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCg                                                                                                                         <  3:367845/149‑1 (MQ=255)
                                ccGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCtt                                                                                                                   >  5:362684/1‑148 (MQ=255)
                                 cGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTtga                                                                                                                 <  3:287011/149‑1 (MQ=255)
                                   gCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTtgatg                                                                                                               <  2:85320/149‑1 (MQ=255)
                                       gAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATc                                                                                                            <  6:358227/148‑1 (MQ=255)
                                        aaCTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATATTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCtt                                                                                                          >  1:39420/1‑149 (MQ=255)
                                         aCTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGaa                                                                                                                            <  6:334920/130‑1 (MQ=255)
                                         aCTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGaa                                                                                                                            >  5:334920/1‑130 (MQ=255)
                                         aCTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTg                                                                                                         <  6:58039/149‑1 (MQ=255)
                                             gctgcGCGGAAAGCGAAAATCGCCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGaa                                                                                                       <  2:300446/147‑1 (MQ=255)
                                                      aaaGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGa                                                                                            <  4:354808/149‑1 (MQ=255)
                                                       aaGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGaa                                                                                           <  6:362684/149‑1 (MQ=255)
                                                           gcaactCGCCAACCTCCGGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCCCCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGACCACCGTATTGAAGCCg                                                                                       <  4:240044/144‑1 (MQ=255)
                                                           gAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                       <  1:181901/149‑1 (MQ=255)
                                                           gAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                       >  3:135626/1‑149 (MQ=255)
                                                           gAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                       >  2:181901/1‑149 (MQ=255)
                                                           gAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                       >  5:271280/1‑149 (MQ=255)
                                                           gAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                       >  8:346628/1‑149 (MQ=255)
                                                            aaaaTCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCgg                                                                                      <  4:327360/149‑1 (MQ=255)
                                                            aaaaTCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCgg                                                                                      <  7:267177/149‑1 (MQ=255)
                                                                tCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCaa                                                                                  <  1:100931/149‑1 (MQ=255)
                                                                  ggCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAgg                                                                                <  5:250875/149‑1 (MQ=255)
                                                                   gCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCTATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGt                                                                               <  7:346628/149‑1 (MQ=255)
                                                                   gCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGt                                                                               <  6:187300/149‑1 (MQ=255)
                                                                   gCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGt                                                                               <  7:286930/149‑1 (MQ=255)
                                                                   gCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGt                                                                               >  7:289257/1‑149 (MQ=255)
                                                                   gCAAACTCCTGGGTGTCGGCCATACGCCACACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGt                                                                               >  4:19231/1‑149 (MQ=255)
                                                                           gtGGGTGTCGGCCAAACGCACCACCCGGTCGACACTCTTCACCATTCGATACTTGTCGATCACCGTGTAGCGATCCACTGTTTTCGATGGCTAGAAGCGGTCATTTATGATCTTTAACACCTTATTTAAGCCGGGCAAGTTGAACACGc                                                                       <  6:319379/148‑1 (MQ=255)
                                                                           cTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGc                                                                       >  3:335763/1‑149 (MQ=255)
                                                                            tGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCAACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCt                                                                      >  4:178310/1‑149 (MQ=255)
                                                                                     gCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCAt                                                             >  7:97951/1‑149 (MQ=255)
                                                                                          aCGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATg                                                                  >  7:262692/1‑139 (MQ=255)
                                                                                          aCGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATg                                                                  <  8:262692/139‑1 (MQ=255)
                                                                                            gCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCg                                                       >  4:13341/1‑148 (MQ=255)
                                                                                            gCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCg                                                       >  7:114599/1‑148 (MQ=255)
                                                                                              cccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCAACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCAc                                                    >  4:91763/1‑149 (MQ=255)
                                                                                              ccaaCCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCCTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCAc                                                    >  8:100547/4‑149 (MQ=255)
                                                                                               ccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCAcc                                                   >  2:328434/1‑149 (MQ=255)
                                                                                                 aCCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCAcccc                                                 <  2:143947/149‑1 (MQ=255)
                                                                                                          acacTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                        <  3:13341/149‑1 (MQ=255)
                                                                                                          acacTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                        <  8:114599/149‑1 (MQ=255)
                                                                                                          acacTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                        <  1:363857/149‑1 (MQ=255)
                                                                                                              tctTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                        <  8:169941/145‑1 (MQ=255)
                                                                                                              tctTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                        >  7:169941/1‑145 (MQ=255)
                                                                                                              tctTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGAc                                     >  6:168416/1‑148 (MQ=255)
                                                                                                                  cACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACaaa                                  <  8:97951/147‑1 (MQ=255)
                                                                                                                   aCCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGc                                <  3:178310/148‑1 (MQ=255)
                                                                                                                     cAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGt                             >  7:21779/1‑149 (MQ=255)
                                                                                                                         cGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTACTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGCTGAACACGATCAt                                                                   >  3:251617/1‑107 (MQ=255)
                                                                                                                               ttGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCActgct                   <  4:135626/149‑1 (MQ=255)
                                                                                                                                 gTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCTCACCCCTGGTGCCATGACAAAGCGGTCATCActgct                   <  3:19231/147‑1 (MQ=255)
                                                                                                                                 gTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCActgct                   <  6:271280/147‑1 (MQ=255)
                                                                                                                                 gTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTgg                 >  8:270640/1‑149 (MQ=255)
                                                                                                                                 gTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTgg                 >  6:106098/1‑149 (MQ=255)
                                                                                                                                 gTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGACCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTgg                 >  6:359093/1‑149 (MQ=255)
                                                                                                                                  tCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGc                <  1:328434/149‑1 (MQ=255)
                                                                                                                                   cgaaaaCCGTGGATCGCGCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGc                >  3:365322/6‑148 (MQ=255)
                                                                                                                                      tCACCGTGGCGCGATCCAGGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAgg            <  4:335763/149‑1 (MQ=255)
                                                                                                                                        aCCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCCTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAggt           <  7:100547/148‑1 (MQ=255)
                                                                                                                                         ccGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAggtgg         <  3:91763/149‑1 (MQ=255)
                                                                                                                                              gcgcGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGTGGTCATCACTGCTGGCCAGGTGGTTGAt    >  6:41932/1‑149 (MQ=255)
                                                                                                                                                gcgATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGtca                           >  3:358831/1‑124 (MQ=255)
                                                                                                                                                gcgATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAggtggt        >  8:121980/1‑143 (MQ=255)
                                                                                                                                                gcgATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAggtggt        <  7:121980/143‑1 (MQ=255)
                                                                                                                                                gcgATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGGTTGAt    >  8:47776/1‑147 (MQ=255)
                                                                                                                                                gcgATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGGTTGATCa  <  1:58916/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGGCGCCACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGGTTGATCA  >  NC_002947/5185257‑5185549

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: ‑