Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,248,056 G→T 5.4% A657A (GCC→GCA PP_3724 ← acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,248,0560GT 5.4% 243.2 / 3.4 74A657A (GCC→GCAPP_3724acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base G (36/34);  new base T (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

GTCATTCATGTCGTTCTCCTTGCGCCTGGCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTCCTTGAGCAGTTCGAC  >  NC_002947/4247908‑4248202
                                                                                                                                                    |                                                                                                                                                    
gTCATTCATGTCGTTCTCCTTGCGCCTGGCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATg                                                                                                                                                      <  4:40979/149‑1 (MQ=255)
gTCATTCATGTCGTTCTCCTTGCGCCTGGCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATg                                                                                                                                                      >  2:28035/1‑149 (MQ=255)
    ttcatGTCGTTCTCCTTGCGCCTGGCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGcc                                                                                                                                                   >  3:216262/1‑148 (MQ=255)
     tcatGTCGTTCTCCTTGCGCCTGGCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGcc                                                                                                                                                   >  5:180048/1‑147 (MQ=255)
         gTCGTTCTCCTTGCGCCTGGCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGAc                                                                                                                                              >  3:213967/1‑148 (MQ=255)
                  cTTGCGCCTGGCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATc                                                                                                                                     <  2:287895/148‑1 (MQ=255)
                       gcCTGGCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGtt                                                                                                                               >  8:201019/1‑149 (MQ=255)
                           ggCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTggcg                                                                                                                          >  7:63207/1‑149 (MQ=255)
                            gCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCACGGATGGCCGAGACCAAGGCATCGAGGTTggcg                                                                                                                          <  6:175217/148‑1 (MQ=255)
                                aGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAg                                                                                                                                  <  2:64224/137‑1 (MQ=255)
                                   cGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACg                                                                                                                  <  4:213967/149‑1 (MQ=255)
                                    ggATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTg                                                                                                                             <  5:46045/137‑1 (MQ=255)
                                    ggATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTg                                                                                                                             >  6:46045/1‑137 (MQ=255)
                                          gcgcGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATTGCCGAGACCAAGGCATCGAGGTTGGTGGCCGGACGGCCACg                                                                                                           >  4:291993/1‑148 (MQ=255)
                                          gcgcGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATTGCCGAGACCAAGGCATCGAGGTTGGTGGCCGGACGGCCACg                                                                                                           >  4:287530/1‑148 (MQ=255)
                                          gcgcGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACg                                                                                                           >  8:235069/1‑148 (MQ=255)
                                            gcgTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGaa                                                                                                        <  4:127980/149‑1 (MQ=255)
                                                   cTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATTGCCGAGACCAAGGCATCGAGGTTGGTGGCCGGACGGCCACGGAAGCCTTgc                                                                                                 <  3:287530/149‑1 (MQ=255)
                                                   cTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATTGCCGAGACCAAGGCATCGAGGTTGGTGGCCGGACGGCCACGGAAGCCTTgc                                                                                                 <  3:291993/149‑1 (MQ=255)
                                                   cTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTgc                                                                                                 <  7:201019/149‑1 (MQ=255)
                                                    tGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATAAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTg                                                                                                  >  8:153997/1‑147 (MQ=255)
                                                      cGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTgcag                                                                                               >  2:10544/1‑148 (MQ=255)
                                                      cGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTgca                                                                                                <  7:235069/147‑1 (MQ=255)
                                                      cGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTgca                                                                                                <  7:153997/147‑1 (MQ=255)
                                                       ggtggtACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTgcag                                                                                               >  2:22789/1‑147 (MQ=255)
                                                                 ttGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGCGCCACGGAAGCCTTGCAGCAGGCTGGCAc                                                                                    <  1:22789/149‑1 (MQ=255)
                                                                    gCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGTCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTgc                                                                                 <  8:146755/148‑1 (MQ=255)
                                                                        cGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAgg                                                                            >  4:134098/1‑149 (MQ=255)
                                                                          aCCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTgcg                                                                                <  5:116470/143‑1 (MQ=255)
                                                                          aCCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTgcg                                                                                >  6:116470/1‑143 (MQ=255)
                                                                            cATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGtt                                                                          >  3:143601/1‑147 (MQ=255)
                                                                             aTCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTcag                                                                       <  3:148180/149‑1 (MQ=255)
                                                                             aTCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTcag                                                                       <  3:134098/149‑1 (MQ=255)
                                                                                   gggTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTc                                                                         >  8:107127/1‑141 (MQ=255)
                                                                                   gggTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTc                                                                         <  7:107127/141‑1 (MQ=255)
                                                                                    ggTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTgcg                                                                <  7:97363/149‑1 (MQ=255)
                                                                                     gTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTgcgg                                                               >  5:178733/1‑149 (MQ=255)
                                                                                     gTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTgcgg                                                               >  5:124731/1‑149 (MQ=255)
                                                                                     gTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTgcgg                                                               >  1:14976/1‑149 (MQ=255)
                                                                                      ttCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTgcggc                                                              >  8:102837/1‑149 (MQ=255)
                                                                                         acaTCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGcgcg                                                           <  4:143601/149‑1 (MQ=255)
                                                                                              cAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATgg                                                      <  8:63207/149‑1 (MQ=255)
                                                                                              cAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCAGCGCGGATgg                                                      <  1:28035/149‑1 (MQ=255)
                                                                                               aGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGaa                                                                                                        >  4:139523/1‑98 (MQ=255)
                                                                                                     aTCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTc                                               <  1:10544/149‑1 (MQ=255)
                                                                                                      tCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCg                                              >  4:252453/1‑149 (MQ=255)
                                                                                                      tCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCg                                              >  3:47996/1‑149 (MQ=255)
                                                                                                        aaCTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCgg                                             >  5:78553/1‑148 (MQ=255)
                                                                                                         aCTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCgg                                             >  8:104659/1‑147 (MQ=255)
                                                                                                         aCTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTa                                           <  4:47996/149‑1 (MQ=255)
                                                                                                          cTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAg                                          >  2:290040/1‑149 (MQ=255)
                                                                                                             cccGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTg                                        >  1:245836/1‑149 (MQ=255)
                                                                                                             cccGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTg                                        >  1:220696/1‑148 (MQ=255)
                                                                                                             cccGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTg                                        <  2:220696/148‑1 (MQ=255)
                                                                                                               cGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCgg                                             >  4:215549/1‑141 (MQ=255)
                                                                                                               cGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCgg                                             <  3:215549/141‑1 (MQ=255)
                                                                                                                ggCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTggg                                      >  4:287229/1‑147 (MQ=255)
                                                                                                                ggCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGgca                                    <  3:287229/149‑1 (MQ=255)
                                                                                                                 gCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGgcag                                   >  3:30092/1‑149 (MQ=255)
                                                                                                                        cGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAg                             <  5:248340/148‑1 (MQ=255)
                                                                                                                        cGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAg                             >  6:248340/1‑148 (MQ=255)
                                                                                                                         gCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGc                           <  7:102837/149‑1 (MQ=255)
                                                                                                                           aGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTg                         <  6:180048/149‑1 (MQ=255)
                                                                                                                            ggCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGa                        <  6:78553/149‑1 (MQ=255)
                                                                                                                                 aaTCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCCGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTg                   >  8:161828/1‑149 (MQ=255)
                                                                                                                                   tCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGCCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTc                 >  6:305272/1‑149 (MQ=255)
                                                                                                                                        gACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAg                             <  3:280222/132‑1 (MQ=255)
                                                                                                                                        gACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAg                             >  4:280222/1‑132 (MQ=255)
                                                                                                                                        gACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTCCTTGa            <  6:178733/149‑1 (MQ=255)
                                                                                                                                        gACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTCCTTGa            <  4:216262/149‑1 (MQ=255)
                                                                                                                                        gACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTCCTTGa            <  1:290040/149‑1 (MQ=255)
                                                                                                                                           tGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTCCTTGAg           >  5:263994/1‑147 (MQ=255)
                                                                                                                                              gcgGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTCCTTGAGCAGtt      <  3:91536/149‑1 (MQ=255)
                                                                                                                                              gcgGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTCCTTGAGCAGt       >  4:91536/1‑148 (MQ=255)
                                                                                                                                                 gATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGTCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTCCTTGAGCAGTTCGa   >  5:230875/1‑149 (MQ=255)
                                                                                                                                                  aTGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTCCTTGAGCAGTTCGAc  <  6:124731/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCATTCATGTCGTTCTCCTTGCGCCTGGCCAAGGCGGATCAGCGCGTCGACTGCGGTGGTACCGTTGGCACCGACCATCAATGGGTTCACATCCAGCTCCATCAACTGCCCGGCGTTTTCGCAGGCGTAATCGGCGACTGCGCGGATGGCCGAGACCAAGGCATCGAGGTTGGCGGCCGGACGGCCACGGAAGCCTTGCAGCAGGCTGGCACTGCGCAGGTTCAGCAGTGCGGCGCGGATGGCGCTGTCGGTAGTGGGCAGCAACAGGCTGACACTGTCCTTGAGCAGTTCGAC  >  NC_002947/4247908‑4248202

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: