Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,685,607 T→C 7.0% 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 7.0% 203.1 / 4.7 71D391G (GAT→GGT) patDmedium chain aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (21/45);  new base C (2/3);  total (23/48)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GCTGACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGC  >  NC_002947/1685459‑1685753
                                                                                                                                                    |                                                                                                                                                   
gCTGACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCAt                                                                                                                                                     >  4:61519/1‑149 (MQ=255)
gCTGACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCAt                                                                                                                                                     >  3:113074/1‑149 (MQ=255)
gCTGACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCAt                                                                                                                                                     >  5:71041/1‑149 (MQ=255)
    gccAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCgg                                                                                                                                                  <  1:270728/147‑1 (MQ=255)
      cAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCg                                                                                                                                                   <  8:53801/145‑1 (MQ=255)
      cAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCg                                                                                                                                                   >  7:53801/1‑145 (MQ=255)
       aGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGaa                                                                                                                                              <  2:10666/149‑1 (MQ=255)
           tgatctGGGTAATCACCCAGGTGCAGACATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCgg                                                                                                                                          <  8:90425/144‑1 (MQ=255)
           tGAAGTGGGTATTCACCGAGGTGCAGGCATACTGCAAGCGCGCCGCCAGGGGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCGGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCACCGATGAAGCgc                                                                                                                                          <  4:355469/149‑2 (MQ=255)
             aaGTGGGTATTCACCCAGTTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCggg                                                                                                                                         <  5:353160/148‑1 (MQ=255)
             aaGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCggg                                                                                                                                         <  8:149250/148‑1 (MQ=255)
             aaGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGt                                                                                                                                        <  7:277325/149‑1 (MQ=255)
              aGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGt                                                                                                                                        <  2:320890/148‑1 (MQ=255)
              aGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCACCGATGAAGCgcagg                                                                                                                                       <  2:38967/149‑5 (MQ=255)
              aGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCACCGATGAAGCgcagg                                                                                                                                       <  3:134360/149‑5 (MQ=255)
               gTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTTGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGa                                                                                                                                      >  6:363517/1‑149 (MQ=255)
                tGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGAc                                                                                                                                     >  2:181633/1‑149 (MQ=255)
                        cACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAacca                                                                                                                             >  8:184934/1‑149 (MQ=255)
                           ccAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCg                                                                                                                          <  4:358653/149‑1 (MQ=255)
                           ccAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCg                                                                                                                          >  8:107839/1‑149 (MQ=255)
                            cAGGTGCAGCCATACTGCCAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGGCGCGGGGCCATACCGATGCGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCgg                                                                                                                         <  1:363523/149‑1 (MQ=255)
                                   aGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTTCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAac                                                                                                                  <  1:73511/149‑1 (MQ=255)
                                      cATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACt                                                                                                               <  6:121620/149‑1 (MQ=255)
                                         aCTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCg                                                                                                            >  2:334500/1‑149 (MQ=255)
                                            gCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTg                                                                                                         >  7:144178/1‑149 (MQ=255)
                                             cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                          <  1:181633/147‑1 (MQ=255)
                                             cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                          <  5:360239/147‑1 (MQ=255)
                                             cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                          <  7:184934/147‑1 (MQ=255)
                                             cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                          <  7:326328/147‑1 (MQ=255)
                                             cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTgc                                                                                                        <  6:6639/149‑1 (MQ=255)
                                              aaGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTgcg                                                                                                       >  4:177827/1‑149 (MQ=255)
                                                   cgccgcCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCATGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACATCGTTGCGCACGa                                                                                                 >  3:257853/1‑149 (MQ=255)
                                                   cgccgcCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGa                                                                                                 <  3:61519/149‑1 (MQ=255)
                                                     ccgccAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGGTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGAtc                                                                                               <  4:157869/149‑1 (MQ=255)
                                                      cgccAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGAtct                                                                                              >  1:236122/1‑149 (MQ=255)
                                                              ggTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTcc                                                                                      >  3:354774/1‑149 (MQ=255)
                                                                      cgGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCACATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGc                                                                              <  6:71041/149‑1 (MQ=255)
                                                                        gCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGt                                                                            <  3:213527/149‑1 (MQ=255)
                                                                         ccGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCCACGATCTCGTCGTCCTGCAGGGCGt                                                                            <  1:196279/149‑1 (MQ=255)
                                                                           ggTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCgg                                                                         <  2:80350/149‑1 (MQ=255)
                                                                               tCGCGGGTCCATACCGATGCGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCa                                                                      <  3:200814/148‑1 (MQ=255)
                                                                                   gggTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACgg                                                                <  1:136507/149‑1 (MQ=255)
                                                                                        cATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAg                                                             <  1:334500/148‑1 (MQ=255)
                                                                                          tACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCt                                                          <  5:210805/149‑1 (MQ=255)
                                                                                           aCCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTAGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTc                                                         >  6:25025/1‑149 (MQ=255)
                                                                                              gATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCgaa                                                      <  5:176284/149‑1 (MQ=255)
                                                                                                         gccGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCAt                                           <  4:354774/149‑1 (MQ=255)
                                                                                                          ccGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATc                                          <  8:144178/149‑1 (MQ=255)
                                                                                                           cGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCg                                         >  5:130660/1‑149 (MQ=255)
                                                                                                             tAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCgaa                                                      >  5:66658/1‑134 (MQ=255)
                                                                                                             tAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCgaa                                                      <  6:66658/134‑1 (MQ=255)
                                                                                                               gTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGAtctc                                                                                             <  2:74753/93‑1 (MQ=255)
                                                                                                               gTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGAtctc                                                                                             >  1:74753/1‑93 (MQ=255)
                                                                                                               gTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCg                                     <  7:107839/149‑1 (MQ=255)
                                                                                                               gTTGGAGGCGTTGGCCCATTCCAGTGCCTGCTCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGGACGATCTGGGCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCg                                     <  5:377190/149‑1 (MQ=255)
                                                                                                                ttGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGc                                                                        <  8:217889/113‑1 (MQ=255)
                                                                                                                ttGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGc                                                                        >  7:217889/1‑113 (MQ=255)
                                                                                                                  ggAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTa                                  >  2:20402/1‑149 (MQ=255)
                                                                                                                   gAGTCGTTGGCCCATTCCTGTGCGTGGTCTTCATCGGTGAAGCGGGTGACCTAAACCTCGGGGCCAAACACTTCGTTGCGAACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTATGCTCGAAGAAAAAGCCATCGCCCGCTa                                  <  3:39770/148‑1 (MQ=255)
                                                                                                                    aGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTAtt                                <  5:363517/149‑1 (MQ=255)
                                                                                                                     gTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTg                               <  8:29534/149‑1 (MQ=255)
                                                                                                                         ttGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCtt                           <  2:256249/149‑1 (MQ=255)
                                                                                                                           ggCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTg                          <  7:34298/148‑1 (MQ=255)
                                                                                                                            gCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGc                         <  1:20402/148‑1 (MQ=255)
                                                                                                                                aTTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCg                    <  2:236122/149‑1 (MQ=255)
                                                                                                                                    cAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTa                >  8:261795/1‑149 (MQ=255)
                                                                                                                                     aGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTAc               <  3:177827/149‑1 (MQ=255)
                                                                                                                                      gTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTAc               >  3:81771/1‑148 (MQ=255)
                                                                                                                                       cgCCTGCTCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCa             <  4:106091/148‑1 (MQ=255)
                                                                                                                                           tGCGCTTCATCGGTGAAGCGGGTGACCGAAACCCCCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGAc          >  5:119183/1‑148 (MQ=255)
                                                                                                                                           tGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTACGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGAc          >  2:253605/1‑148 (MQ=255)
                                                                                                                                            ccgcTTCACCGATGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACg         >  2:41543/2‑148 (MQ=255)
                                                                                                                                            ccgcTTCACCGATGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACg         >  1:337324/2‑148 (MQ=255)
                                                                                                                                                   aTCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGc  <  4:113074/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTGACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGC  >  NC_002947/1685459‑1685753

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: