Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,262,496 T→A 6.2% T223T (ACT→ACA PP_4638 → methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,262,4960TA 6.2% 276.5 / 3.5 89T223T (ACT→ACAPP_4638methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (38/45);  new base A (4/2);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCTTGCAGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAGGGGTTCA  >  NC_002947/5262348‑5262642
                                                                                                                                                    |                                                                                                                                                  
tCACCTTGCAGCTGGGCTTCGATGCCTCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCCTGAGCCTCACCCACGGGCGCGCGACGAGGCTCGCCCACGTGGCCGGCCTGACt                                                                                                                                                    >  6:336501/1‑149 (MQ=255)
tCACCTTGCAGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACt                                                                                                                                                    >  5:320519/1‑149 (MQ=255)
 cACCTTGCAGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTg                                                                                                                                                   >  1:271498/1‑149 (MQ=255)
 cACCTTGCAGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTg                                                                                                                                                   >  7:63361/1‑149 (MQ=255)
    cTTGCAGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCa                                                                                                                                                >  1:56684/1‑149 (MQ=255)
        cAGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACAGTCACGcc                                                                                                                                            <  4:159698/149‑1 (MQ=255)
           cTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGccatc                                                                                                                                         >  3:187713/1‑149 (MQ=255)
              ggCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGccatccat                                                                                                                                      >  3:6351/1‑149 (MQ=255)
                   cGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATgcgc                                                                                                                                  <  5:336501/148‑1 (MQ=255)
                     aTGCCGCCAAGCATCAGGAAGCCATACACTGGTTGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGAc                                                                                                                               <  7:319745/149‑1 (MQ=255)
                     aTGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGAc                                                                                                                               <  1:80288/149‑1 (MQ=255)
                        ccgccAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATg                                                                                                                            <  4:70473/149‑1 (MQ=255)
                          gccAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCt                                                                                                                          <  2:56684/149‑1 (MQ=255)
                          gccAACCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGATGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCCCGCCATCCATGCGCGACATGCt                                                                                                                          >  6:135398/1‑149 (MQ=255)
                            cAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTc                                                                                                                         >  2:178581/1‑148 (MQ=255)
                              aGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCg                                                                                                                        >  3:221836/1‑147 (MQ=255)
                                    aGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCg                                                                                                                   >  8:273820/1‑146 (MQ=255)
                                    aGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGga                                                                                                                 <  5:107896/148‑1 (MQ=255)
                                    aGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATCCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCTAATCGgag                                                                                                                <  8:242346/149‑1 (MQ=255)
                                     ggAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGgaga                                                                                                               <  2:371982/149‑1 (MQ=255)
                                            aTACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGgcg                                                                                                          >  4:19416/1‑147 (MQ=255)
                                              acacTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGt                                                                                                      >  2:321172/1‑149 (MQ=255)
                                                acTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTc                                                                                                     >  6:265151/1‑148 (MQ=255)
                                                acTGGATGCGCCGTCAGGCCTCACCCAAGACCATGCTGGCCATCCTCATGAGCCTCACCCACGGGCGCGCGACGAAGCTCGACCACGTGGCCGGCCTGACAGTCACGCCATCCATGCGCGACATGCTCTAAGCGGACACGGCGCAGTc                                                                                                     >  8:79572/1‑148 (MQ=255)
                                                 cTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTcc                                                                                                    >  3:57878/1‑148 (MQ=255)
                                                  tGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGCTGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTcc                                                                                                    >  7:297029/1‑147 (MQ=255)
                                                    gATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCGACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCaa                                                                                                  <  2:186551/147‑1 (MQ=255)
                                                    gATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAgg                                                                                                >  1:222840/1‑149 (MQ=255)
                                                    gATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAgg                                                                                                <  4:28746/149‑1 (MQ=255)
                                                    gATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAg                                                                                                 <  4:57878/148‑1 (MQ=255)
                                                    gATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAg                                                                                                 <  2:271498/148‑1 (MQ=255)
                                                         gcCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCttt                                                                                           <  4:80034/149‑1 (MQ=255)
                                                                ggCCTCACCCAAGCCCATGCTGGCCTGCCTCATTACCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAgg                                                                                    <  4:215865/149‑1 (MQ=255)
                                                                 gCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAgg                                                                                    >  5:56135/1‑148 (MQ=255)
                                                                 gCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAgg                                                                                    <  6:56135/148‑1 (MQ=255)
                                                                 gCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGgc                                                                                   <  4:181749/149‑1 (MQ=255)
                                                                 gCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGgc                                                                                   <  1:178581/149‑1 (MQ=255)
                                                                   cTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAg                                                                                     >  3:35173/1‑145 (MQ=255)
                                                                   cTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAg                                                                                     <  4:35173/145‑1 (MQ=255)
                                                                    tCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGgcgcg                                                                                >  5:259168/1‑149 (MQ=255)
                                                                        ccAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAgt                                                                            >  5:228759/1‑149 (MQ=255)
                                                                          caGCCCATGCTGGCGTGCCTACTGAGCCTCTCCCACGTGGGCGCGACGATGCTCGACCACGTGGCCGTCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGAAgtgt                                                                          <  8:63361/148‑1 (MQ=255)
                                                                           aGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGtt                                                                         <  4:221836/149‑1 (MQ=255)
                                                                           aGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGtt                                                                         <  1:321172/149‑1 (MQ=255)
                                                                           aGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGtt                                                                         <  7:273820/149‑1 (MQ=255)
                                                                            gCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGtt                                                                         >  4:95327/1‑148 (MQ=255)
                                                                            gCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTg                                                                        <  6:259168/149‑1 (MQ=255)
                                                                              ccATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGAc                                                                      >  3:244766/1‑149 (MQ=255)
                                                                                     ggCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTgg                                                               <  8:297029/149‑1 (MQ=255)
                                                                                         tgcctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATGCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGATGACCGCCTGGCCtt                                                           >  8:304756/1‑149 (MQ=255)
                                                                                           cctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGTCCAGGCGCGCAGTGTTGACCTCCTGGCCtt                                                           <  4:150322/147‑1 (MQ=255)
                                                                                           cctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTg                                                          >  5:283573/1‑148 (MQ=255)
                                                                                           cctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTg                                                          <  5:139964/148‑1 (MQ=255)
                                                                                           cctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGc                                                         <  3:19416/149‑1 (MQ=255)
                                                                                            ctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGACGCGGAGCCGGCGCAGTCCAAGGCCTTTGCCCGGGCGCGCAGGGTGGAGCGCCTGGCCTTGcc                                                        >  8:248940/1‑148 (MQ=255)
                                                                                            ctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCa                                                        <  5:135398/149‑1 (MQ=255)
                                                                                            ctCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCCCAGTGTTGACCGCCTGGCCTTGCa                                                        <  7:54916/149‑1 (MQ=255)
                                                                                               aTGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATGCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGa                                                      <  7:304756/148‑1 (MQ=255)
                                                                                                 gAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTa                                                   >  2:355910/1‑149 (MQ=255)
                                                                                                  aGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTa                                                   <  5:265151/148‑1 (MQ=255)
                                                                                                  aGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTa                                                   <  1:38375/148‑1 (MQ=255)
                                                                                                   gCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTAt                                                  <  4:6351/148‑1 (MQ=255)
                                                                                                      tCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATc                                                 >  3:155994/1‑146 (MQ=255)
                                                                                                      tCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATc                                                 <  4:155994/146‑1 (MQ=255)
                                                                                                       cACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGt                                              >  2:205785/1‑148 (MQ=255)
                                                                                                        aCCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGAt                                                     <  5:368975/140‑1 (MQ=255)
                                                                                                        aCCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACAAGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGAt                                                     >  6:368975/1‑140 (MQ=255)
                                                                                                          ccACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACa                                           >  2:130655/1‑148 (MQ=255)
                                                                                                            aCGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAg                                        <  3:95327/149‑1 (MQ=255)
                                                                                                              gggCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCt                                      >  4:339240/1‑149 (MQ=255)
                                                                                                                    gcgACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATgt                                 >  1:333397/1‑147 (MQ=255)
                                                                                                                      gacgaTGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGc                               >  5:114623/1‑148 (MQ=255)
                                                                                                                      gacgaTGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGc                               >  8:167296/1‑148 (MQ=255)
                                                                                                                        cgaTGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTAc                            <  6:320519/149‑1 (MQ=255)
                                                                                                                             cTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCg                        <  6:283573/148‑1 (MQ=255)
                                                                                                                              tCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGgt                      >  1:344322/1‑149 (MQ=255)
                                                                                                                              tCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGAGCGCAGAGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATCGGCTACGCCGgt                      >  4:232038/1‑149 (MQ=255)
                                                                                                                                gcgcccGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGgtgt                    <  1:235474/144‑1 (MQ=255)
                                                                                                                                gACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGgt                      <  3:339240/147‑1 (MQ=255)
                                                                                                                                 aCCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGgtgt                    >  1:185284/1‑148 (MQ=255)
                                                                                                                                 aCCACGTGGCCGGCGTGAATGTCACGACATCCATGCGCGACATGATCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGtgtg                   <  2:333397/149‑1 (MQ=255)
                                                                                                                                   cACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCa                 <  4:323804/149‑1 (MQ=255)
                                                                                                                                   cACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCa                 <  2:255833/149‑1 (MQ=255)
                                                                                                                                          atggCGTGACAGTCACGCCGTCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGt            >  2:100367/3‑147 (MQ=255)
                                                                                                                                          atggCGTGACAGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGt            >  5:243128/3‑147 (MQ=255)
                                                                                                                                          atggCGTGACAGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATCGGCTACGCCGGTGTGCACCTGt            >  2:35522/3‑147 (MQ=255)
                                                                                                                                           cGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAg         >  4:59236/1‑149 (MQ=255)
                                                                                                                                            ggCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAgg        >  8:241271/1‑149 (MQ=255)
                                                                                                                                             gCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAggg       <  4:187713/149‑1 (MQ=255)
                                                                                                                                                tGACAGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATCGGCTACGCCGGTGTGCACCTGTCAGGGGtt    <  1:35522/149‑1 (MQ=255)
                                                                                                                                                   cTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAGGGGTTca  <  7:167296/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCACCTTGCAGCTGGGCTTCGATGCCGCCAAGCATCAGGAAGCCATACACTGGATGCGCCGTCAGGCCTCACCCAAGCCCATGCTGGCCTGCCTCATGAGCCTCACCCACGGGCGCGCGACGATGCTCGACCACGTGGCCGGCGTGACTGTCACGCCATCCATGCGCGACATGCTCGAAGCGGAGACGGCGCAGTCCAAGGCCTTTGCCCAGGCGCGCAGTGTTGACCGCCTGGCCTTGCAGATTATCGGTACAAAGCTGATGGGCTACGCCGGTGTGCACCTGTCAGGGGTTCA  >  NC_002947/5262348‑5262642

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: