Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,185,400 C→A 7.4% A354S (GCC→TCC)  aceK ← isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,185,4000CA 7.4% 241.4 / 5.0 78A354S (GCC→TCC) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base A (2/4);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GTCGACGGCGCCACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGGTTGAT  >  NC_002947/5185252‑5185547
                                                                                                                                                    |                                                                                                                                                   
gTCGACGGCGCCACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGgc                                                                                                                                                     >  5:137557/1‑149 (MQ=255)
     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                                                                                                                                          <  8:298303/139‑1 (MQ=255)
                     aaCTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACgg                                                                                                                                          >  7:298303/1‑139 (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                                                                                     >  2:181901/1‑149 (MQ=255)
                                                                gAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                     >  3:135626/1‑149 (MQ=255)
                                                                gAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                     >  8:346628/1‑149 (MQ=255)
                                                                gAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                     >  5:271280/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                                                                             >  7:289257/1‑149 (MQ=255)
                                                                        gCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGt                                                                             <  6:187300/149‑1 (MQ=255)
                                                                        gCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGt                                                                             <  7:286930/149‑1 (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                                                                <  8:262692/139‑1 (MQ=255)
                                                                                               aCGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATg                                                                >  7:262692/1‑139 (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                                      <  8:114599/149‑1 (MQ=255)
                                                                                                               acacTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                      <  3:13341/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               >  6:106098/1‑149 (MQ=255)
                                                                                                                                      gTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTgg               >  8:270640/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)
                                                                                                                                                    |                                                                                                                                                   
GTCGACGGCGCCACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGGTTGAT  >  NC_002947/5185252‑5185547

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