Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,457,059 T→G 8.2% E77A (GAG→GCG)  hbdH ← 3‑hydroxybutyrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,457,0590TG 8.2% 247.8 / 9.2 83E77A (GAG→GCG) hbdH3‑hydroxybutyrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base G (5/2);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGC  >  NC_002947/3456915‑3457207
                                                                                                                                                |                                                                                                                                                    
cccGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTcgcg                                                                                                                                                  <  7:140406/149‑1 (MQ=255)
cccGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTcgc                                                                                                                                                   <  3:166064/148‑1 (MQ=255)
 ccGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGt                                                                                                                                                 >  6:105250/1‑149 (MQ=255)
   ggCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTc                                                                                                                                               <  4:53434/149‑1 (MQ=255)
    gCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCg                                                                                                                                              <  4:345888/149‑1 (MQ=255)
          ccAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATAATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCAACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAAAACGCGTTCGGCCAgt                                                                                                                                        >  8:56754/1‑148 (MQ=255)
           cAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGt                                                                                                                                       <  5:152526/149‑1 (MQ=255)
           cAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCAAGCTTTATGGAGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACACGCGTAAGGCCAgct                                                                                                                                       >  1:333731/1‑147 (MQ=255)
             ggCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATTATCTTGTCCCAGCTTTCTGGCGGTAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGTATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTg                                                                                                                                     <  7:56754/149‑1 (MQ=255)
              gcgcGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTg                                                                                                                                     >  6:294007/1‑148 (MQ=255)
                gcgTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGaa                                                                                                                                   >  2:176249/1‑148 (MQ=255)
                      cATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGcc                                                                                                                            >  3:60645/1‑149 (MQ=255)
                        tGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTc                                                                                                                          <  6:327067/149‑1 (MQ=255)
                        tGGAATACGGCCGACAGGTTCAGGGCGATGATCTTCTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACAGGCGTTGTTCACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACACAGCCTc                                                                                                                          >  7:188235/1‑149 (MQ=255)
                           aaTACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGATCGGCCAGGTTGAACAAAGCCTCGa                                                                                                                        >  7:198337/1‑148 (MQ=255)
                               cGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTgg                                                                                                                   <  2:179036/149‑1 (MQ=255)
                               cGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTgg                                                                                                                   <  3:30631/149‑1 (MQ=255)
                                ggCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGa                                                                                                                  <  7:302717/149‑1 (MQ=255)
                                ggCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGa                                                                                                                  >  2:14800/1‑149 (MQ=255)
                                ggCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACGCGAGTTCGGcg                                                                                                                                           >  4:28904/1‑123 (MQ=255)
                                ggCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACGCGAGTTCGGcg                                                                                                                                           <  3:28904/124‑2 (MQ=255)
                                ggCCGACAGGTTCAGGGCGAAGATCTTGTCCCAGCTGTCTGGCGGGAACTGCTCCACCGGTGCCACATGCTGGATACCGGCGGGGTTGACCAGGAAGTCGACGCCGCCGAACTCGCGTACGGCCAGGAGGAACAAAGCCTCGCTCTggc                                                                                                                  >  5:226417/1‑148 (MQ=255)
                                  ccGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGAc                                                                                                                 >  1:22757/1‑148 (MQ=255)
                                       tctttCATGTCGATGATCTTGTCCCAGCTTTCTGGCGGTAACTGCTCAACGGGTGCCACTTGCTGTATACCGGCGTTGTTGACCAGGATGTCGACGGCGCAGAACTCGCGTTCGGCCTGGTTGAACAAAGCCTCGATCTGGACCACAt                                                                                                            <  7:68776/145‑1 (MQ=255)
                                        ggTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATcc                                                                                                          <  5:39182/149‑1 (MQ=255)
                                           taaGGGAGATGATATTGTCCCAGCTTTCTGTCGGGAACTGATCAACAGTTGCCACATTCTGCATACAGTCGTTGTTGACAAGGATGTAGACGCCGCCGAACTCGCGTTCTGCAAGGTTGAACAAAGACTCGATCTGCACCACAt                                                                                                            >  8:68776/3‑144 (MQ=255)
                                              gggCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTg                                                                                                                                     >  6:187678/1‑116 (MQ=255)
                                              gggCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTg                                                                                                                                     <  5:187678/116‑1 (MQ=255)
                                              gggCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAgg                                                                                                    >  1:160686/1‑149 (MQ=255)
                                               ggCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACAt                                                                                                            >  6:195518/1‑140 (MQ=255)
                                               ggCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACAt                                                                                                            <  5:195518/140‑1 (MQ=255)
                                                gCGATGTTCTTGTCCCAGGTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTc                                                                                                  <  3:157712/149‑1 (MQ=255)
                                                       tCTTGTCCCAGCTTTCGGGCTGGAACTGCTCACCCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGCTGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAACGCCTCGATCTGGACCACATCCGACAGGTCggactg                                                                                            <  6:226417/148‑5 (MQ=255)
                                                           gTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGggtggt                                                                                       >  3:175482/1‑149 (MQ=255)
                                                            tCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGggtggt                                                                                       <  1:14800/148‑1 (MQ=255)
                                                            tCCCAGCTTTCAGGCGGGAACGGCTCAACCGGTGCCACATGCGGGATACCGGCGTTGTTGACCAGCATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGggtggt                                                                                       <  8:143265/148‑1 (MQ=255)
                                                              ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATcc                                                                                                          >  7:244352/1‑127 (MQ=255)
                                                              ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATcc                                                                                                          <  8:244352/127‑1 (MQ=255)
                                                              ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGgtggtg                                                                                      >  3:226378/1‑147 (MQ=255)
                                                              ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGgtggtg                                                                                      <  4:226378/147‑1 (MQ=255)
                                                              ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGca                                                                                    >  5:106978/1‑149 (MQ=255)
                                                              ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGc                                                                                     <  5:294007/148‑1 (MQ=255)
                                                               cAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGca                                                                                    >  1:47136/1‑148 (MQ=255)
                                                                         gCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCTTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCAc                                                                          <  5:70643/148‑1 (MQ=255)
                                                                         gCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCAc                                                                          <  1:149544/148‑1 (MQ=255)
                                                                         gCGGGAACTGCTCAACCGGTGCCACATGCTGGATACAGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCAcc                                                                         >  5:179558/1‑149 (MQ=255)
                                                                           gggAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCAcccc                                                                       >  6:133098/1‑149 (MQ=255)
                                                                               aCTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGc                                                                   >  7:262731/1‑149 (MQ=255)
                                                                                   cTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCgg                                                                <  6:106978/148‑1 (MQ=255)
                                                                                              ccACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGc                                                     >  3:203859/1‑148 (MQ=255)
                                                                                                 caTGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATg                                                 >  5:350132/1‑149 (MQ=255)
                                                                                                   tGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGc                                               >  4:327360/1‑149 (MQ=255)
                                                                                                   tGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGc                                               >  1:219259/1‑149 (MQ=255)
                                                                                                   tGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGc                                               >  7:181034/1‑149 (MQ=255)
                                                                                                      tGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCggg                                            >  7:176004/1‑149 (MQ=255)
                                                                                                          tACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGccc                                        <  2:22757/149‑1 (MQ=255)
                                                                                                          tACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGccc                                        >  5:190113/1‑149 (MQ=255)
                                                                                                              ggCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAAc                                                                                <  4:174724/105‑1 (MQ=255)
                                                                                                              ggCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAAc                                                                                >  3:174724/1‑105 (MQ=255)
                                                                                                                  ttgttgACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCAcc                                <  2:219259/149‑1 (MQ=255)
                                                                                                                    gttgACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACCACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCCCCa                               >  4:3436/1‑148 (MQ=255)
                                                                                                                         ccAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCg                         >  2:27/1‑149 (MQ=255)
                                                                                                                          cAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGt                        >  3:241605/1‑149 (MQ=255)
                                                                                                                           agtttGTCGAGGCCTCCGAACTCGCGTTGGGCCAGGTTGAACAAAGCCTCTATCTGGACCACATCCGACAGGTCGGCCGGGTGTTGCACAACCTTCTCCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAATCCGtt                       <  2:333731/145‑1 (MQ=255)
                                                                                                                           aGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGtt                       >  6:52419/1‑149 (MQ=255)
                                                                                                                            ggATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTg                      <  6:179558/149‑1 (MQ=255)
                                                                                                                             gATGTCGACGGCGCGGAACTCGGGGTCGGCCAGGTTGAACAATGCCTCGATCTGGAGCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGa                     <  8:198337/149‑1 (MQ=255)
                                                                                                                                gTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCa                  >  1:285457/1‑149 (MQ=255)
                                                                                                                                gTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCa                  >  3:3260/1‑149 (MQ=255)
                                                                                                                                    aCGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGAt              <  2:47136/149‑1 (MQ=255)
                                                                                                                                    aCGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGAt              <  8:181034/149‑1 (MQ=255)
                                                                                                                                     cctggccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATg             >  7:212802/5‑149 (MQ=255)
                                                                                                                                      gccgccGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGt            <  8:262731/149‑1 (MQ=255)
                                                                                                                                       tggccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGtt           >  4:101246/3‑149 (MQ=255)
                                                                                                                                       tggccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGtt           >  1:312997/3‑149 (MQ=255)
                                                                                                                                       ccgccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGtt           <  6:350132/149‑1 (MQ=255)
                                                                                                                                        ggccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGACATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTg          >  3:12536/2‑149 (MQ=255)
                                                                                                                                        cgccgAACTCGCGTTCGGCCAGGTTGAACACAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAACGCCGTTGAGCACGATGTTg          >  3:295804/1‑149 (MQ=255)
                                                                                                                                        cgccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTg          <  5:52419/149‑1 (MQ=255)
                                                                                                                                         gccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTgg         >  2:32411/1‑149 (MQ=255)
                                                                                                                                         gccgAACTCGCGTTCGGCCAGGTAGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCTCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTgg         >  6:219012/1‑149 (MQ=255)
                                                                                                                                          ccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTTgc        <  2:312997/149‑1 (MQ=255)
                                                                                                                                           cgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGgcg       <  4:203859/149‑1 (MQ=255)
                                                                                                                                                tCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGc  >  3:28417/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGC  >  NC_002947/3456915‑3457207

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