Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,685,607 T→C 5.8% D391G (GAT→GGT)  patD ← medium chain aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,685,6070TC 5.8% 203.0 / 2.1 69D391G (GAT→GGT) patDmedium chain aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base C (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

TGACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGCGG  >  NC_002947/1685461‑1685755
                                                                                                                                                   |                                                                                                                                                    
tGACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCATTGGCCCATTCCAGTGCCTGCGCTTCATCg                                                                                                                                                    >  7:283492/1‑149 (MQ=255)
 gACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCgg                                                                                                                                                   >  3:188728/1‑149 (MQ=255)
  aCCAGCATGACGTGGGTATTCACCCAGGTGCAGCCATCCTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGGCGCGGGGCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCg                                                                                                                                                    <  5:246018/147‑1 (MQ=255)
      gCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCCCGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGaa                                                                                                                                               <  1:266270/149‑1 (MQ=255)
        aTGACGTGTGTATTCACCCAGGTGCAGCCATCCTGCAAGCTCGCCGCCAGGCGGTGGGCGCGGCCGGGGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTGATCGGTGAAGCg                                                                                                                                            <  6:202549/149‑1 (MQ=255)
           aaGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCACCGATGAAGCgcag                                                                                                                                         <  1:226054/149‑4 (MQ=255)
             gTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGa                                                                                                                                       >  5:295762/1‑149 (MQ=255)
                    ttCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAac                                                                                                                                >  6:127329/1‑149 (MQ=255)
                    ttCACCCAGGTGCAGCCATACTGCAAGCGCGACGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGaaa                                                                                                                                 <  8:265186/148‑1 (MQ=255)
                      cACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAacca                                                                                                                              >  3:13845/1‑149 (MQ=255)
                        cccAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAaccacc                                                                                                                            <  5:127329/149‑1 (MQ=255)
                         ccAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCg                                                                                                                           <  1:327215/149‑1 (MQ=255)
                         ccAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCg                                                                                                                           <  5:177111/149‑1 (MQ=255)
                         ccAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCACCGATGAAGCgcagg                                                                                                                                        >  1:30319/1‑132 (MQ=255)
                         ccAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCACCGATGAAGCgcagg                                                                                                                                        <  2:30319/136‑5 (MQ=255)
                          cAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCgg                                                                                                                          <  3:220187/149‑1 (MQ=255)
                             gTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCATTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCAGGTGACCGAAACCACCGGGcc                                                                                                                       <  8:283492/149‑1 (MQ=255)
                                        cTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGt                                                                                                            >  1:342033/1‑149 (MQ=255)
                                         tGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAac                                                                                                                   <  1:132228/141‑1 (MQ=255)
                                         tGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAac                                                                                                                   >  2:132228/1‑141 (MQ=255)
                                         tGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                           >  2:84558/1‑149 (MQ=255)
                                         tGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                           >  7:212325/1‑149 (MQ=255)
                                           cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCgg                                                                                                                                           <  4:158695/115‑1 (MQ=255)
                                           cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCgg                                                                                                                                           >  3:158695/1‑115 (MQ=255)
                                           cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                           <  4:188728/147‑1 (MQ=255)
                                            aaGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTgcg                                                                                                        <  4:13845/149‑1 (MQ=255)
                                                  gccgccAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGAt                                                                                                  >  1:93565/1‑149 (MQ=255)
                                                     gccAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGAtctc                                                                                               <  8:212325/149‑1 (MQ=255)
                                                     gccAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGAtctc                                                                                               <  1:84558/149‑1 (MQ=255)
                                                     gccAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGAtctc                                                                                               <  6:295762/149‑1 (MQ=255)
                                                        aGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCtcgtc                                                                                            >  4:293454/1‑149 (MQ=255)
                                                         ggCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCtcgtcg                                                                                           >  6:228221/1‑149 (MQ=255)
                                                          gCGGTGGGCGGGGCCGGTGTCGCTTGTCCATACCGATGAGGCAAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCtcgtc                                                                                            <  7:64664/147‑1 (MQ=255)
                                                           cGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCgtcgtc                                                                                         <  5:228221/149‑1 (MQ=255)
                                                            ggtggtCGCGGCAGGTGTAGCTAGTACATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCCCCGGGCCAAACCCTTCGTTGCGCTCGATCTCGTCGTcc                                                                                        >  8:142758/1‑149 (MQ=255)
                                                                      gCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGCAACTACCGGGCCAAACACCTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGt                                                                              >  7:16647/1‑149 (MQ=255)
                                                                             tCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTgcgc                                                                                                       >  7:152909/1‑117 (MQ=255)
                                                                             tCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTgcgc                                                                                                       <  8:152909/117‑1 (MQ=255)
                                                                             tCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCa                                                                        <  2:342033/148‑1 (MQ=255)
                                                                             tCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAg                                                                       >  3:237824/1‑149 (MQ=255)
                                                                                  ggTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACgg                                                                  <  4:237824/149‑1 (MQ=255)
                                                                                   gTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGt                                                                 <  2:93565/149‑1 (MQ=255)
                                                                                       aTACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGc                                                             >  8:80156/1‑149 (MQ=255)
                                                                                           cGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCg                                                          >  4:81836/1‑148 (MQ=255)
                                                                                            gATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCgaa                                                        <  1:164819/149‑1 (MQ=255)
                                                                                                  gccaggccGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCTGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGaaaaa                                                  <  5:325754/149‑1 (MQ=255)
                                                                                                  gccaggccGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGaaaa                                                   >  6:325754/1‑148 (MQ=255)
                                                                                                    ccaggccGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCATCGGCCAGCACGGTAGGCTCGAAGAAAAAg                                                 <  3:81836/149‑1 (MQ=255)
                                                                                                        gccGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCAt                                             <  7:142758/149‑1 (MQ=255)
                                                                                                         ccGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATc                                            <  3:293454/149‑1 (MQ=255)
                                                                                                           gTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAg                                                 <  6:329130/142‑1 (MQ=255)
                                                                                                           gTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAg                                                 >  5:329130/1‑142 (MQ=255)
                                                                                                               ttGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGGAACCACCGGGCCAAACACGTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCGCGAAGAAAACGCGCTCGCCCGc                                      >  4:212031/1‑149 (MQ=255)
                                                                                                                    gTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTg                                 >  3:206327/1‑149 (MQ=255)
                                                                                                                     tCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTgg                                <  7:80156/149‑1 (MQ=255)
                                                                                                                       gTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCt                              >  4:123832/1‑149 (MQ=255)
                                                                                                                           gCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACTACCGGGCCAAACACCTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGc                           <  8:16647/148‑1 (MQ=255)
                                                                                                                             ccATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCAc                        <  4:206327/149‑1 (MQ=255)
                                                                                                                             ccATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGCCTTGCCAc                        >  8:225559/1‑149 (MQ=255)
                                                                                                                                   cAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCg                                                                            <  7:85969/91‑1 (MQ=255)
                                                                                                                                   cAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCg                                                                            >  8:85969/1‑91 (MQ=255)
                                                                                                                                        ccTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCgaa                                                        <  2:62719/106‑1 (MQ=255)
                                                                                                                                        ccTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCgaa                                                        >  1:62719/1‑106 (MQ=255)
                                                                                                                                        ccTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGa             >  8:100705/1‑149 (MQ=255)
                                                                                                                                           ccgcTTCACCGATGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCAGCTATTGGCTTGACACCGGGTACCAGACg           >  7:308313/2‑148 (MQ=255)
                                                                                                                                               ttCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAAtgt      >  3:48874/1‑149 (MQ=255)
                                                                                                                                               ttCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAAtgt      >  2:132252/1‑149 (MQ=255)
                                                                                                                                                  aTCGGTGAAGCGGGTTACCGAAACCACCTGGCCAAAGACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCTGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGc    <  3:290680/148‑1 (MQ=255)
                                                                                                                                                  aTCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGc    <  3:123832/148‑1 (MQ=255)
                                                                                                                                                   tCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGCgg  >  8:2478/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGCGG  >  NC_002947/1685461‑1685755

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: