Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,003,980 A→G 6.9% E66E (GAA→GAG allE → S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,003,9800AG 6.9% 585.3 / 20.9 189E66E (GAA→GAGallES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (76/100);  new base G (8/5);  total (84/105)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 197 total aligned reads displayed.

ACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAACACCACACCCTGGTACCTGGTGG  >  NC_002947/4003831‑4004129
                                                                                                                                                     |                                                                                                                                                     
aCCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGaa                                                                                                                                                       >  2:976250/1‑150 (MQ=255)
  cGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAAc                                                                                                                                                     >  1:1060396/1‑150 (MQ=255)
     gtcgatgactgacCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGtt                                                                                                                                                  <  1:5079356/145‑1 (MQ=255)
     gctgctGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGtt                                                                                                                                                  <  1:590523/150‑1 (MQ=255)
       tcctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtct                                                                                                                                                >  1:861330/3‑150 (MQ=255)
        gctgACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTtctc                                                                                                                                               >  1:2726046/1‑150 (MQ=255)
         ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                                                                                                              <  1:4612988/150‑1 (MQ=255)
         ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                                                                                                              <  1:3341939/150‑1 (MQ=255)
         ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                                                                                                              <  1:2827515/150‑1 (MQ=255)
         ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                                                                                                              <  1:4612866/150‑1 (MQ=255)
         ctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtctcg                                                                                                                                              >  2:2602222/1‑149 (MQ=255)
         ctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtctcg                                                                                                                                              >  2:3070560/1‑149 (MQ=255)
          tgaccgatCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtctcgc                                                                                                                                             >  2:174779/1‑148 (MQ=255)
           gactgacCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCa                                                                                                                                            <  2:3153213/150‑1 (MQ=255)
              cgatcgcgCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAATa                                                                                                                                         >  1:3127040/8‑150 (MQ=255)
                atcgcgCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAATata                                                                                                                                       >  1:1590975/6‑148 (MQ=255)
                 ccGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAt                                                                                                                                      <  2:3115458/150‑1 (MQ=255)
                  cgcgCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatatt                                                                                                                                     <  1:1887276/147‑2 (MQ=255)
                  cGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtc                                                                                                                                     <  2:3612626/150‑1 (MQ=255)
                  cGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCACAGTACAtc                                                                                                                                     >  2:3831942/1‑150 (MQ=255)
                     gCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatg                                                                                                                                  <  1:739164/150‑5 (MQ=255)
                     gCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtcgtc                                                                                                                                  <  1:3871359/150‑1 (MQ=255)
                     gCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtcgtc                                                                                                                                  <  2:3055187/150‑1 (MQ=255)
                     gCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtcgtc                                                                                                                                  <  2:1931985/150‑1 (MQ=255)
                      ccATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatgg                                                                                                                                 <  1:2254002/150‑6 (MQ=255)
                        aTGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGaa                                                                                                                               >  2:920635/1‑150 (MQ=255)
                         tGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAAc                                                                                                                              <  1:4766323/150‑1 (MQ=255)
                           ttCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTg                                                                                                                            >  2:5369028/1‑150 (MQ=255)
                           ttCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTg                                                                                                                            >  2:815001/1‑150 (MQ=255)
                                  aaGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTCCATCGTCGACCTGGCCCCgg                                                                                                                     >  2:3311249/1‑150 (MQ=255)
                                    gCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGaa                                                                                                                   >  1:996968/1‑150 (MQ=255)
                                       tATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAggc                                                                                                                <  1:2932358/150‑1 (MQ=255)
                                       tATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAggc                                                                                                                >  1:2631921/1‑150 (MQ=255)
                                          gCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAggcggc                                                                                                             <  2:5024543/150‑1 (MQ=255)
                                            cGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCag                                                                                                           >  2:1354790/1‑150 (MQ=255)
                                             gTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagc                                                                                                          <  1:3584146/150‑1 (MQ=255)
                                             gTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagc                                                                                                          >  2:2407750/1‑150 (MQ=255)
                                               cATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagcga                                                                                                        <  1:5227830/150‑1 (MQ=255)
                                                aTCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagcgag                                                                                                       >  1:4830930/1‑150 (MQ=255)
                                                  ccccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagcgagcg                                                                                                     <  2:4821180/150‑1 (MQ=255)
                                                   cccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAgcgc                                                                                                    <  1:3911510/150‑1 (MQ=255)
                                                    ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                                                                                                   >  2:3021254/1‑150 (MQ=255)
                                                    ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                                                                                                   >  2:1086081/1‑150 (MQ=255)
                                                    ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                                                                                                   >  1:5360236/1‑150 (MQ=255)
                                                    ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                                                                                                   >  1:1244777/1‑150 (MQ=255)
                                                     cAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGccc                                                                                                  >  1:1421481/1‑150 (MQ=255)
                                                     cAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGccc                                                                                                  >  1:5017976/1‑150 (MQ=255)
                                                           cGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCt                                                                                            <  2:3538939/150‑1 (MQ=255)
                                                                     gACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGAt                                                                                  >  2:2614330/1‑150 (MQ=255)
                                                                       cATCGGCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGCGGGGTATCGCCCTCCCGCTGACCGGTTTTGCCGAACCGGTCTCCCAGGACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGc                                                                                <  2:214068/150‑1 (MQ=255)
                                                                        aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                                                               >  1:3431845/1‑150 (MQ=255)
                                                                        aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                                                               >  2:2642887/1‑150 (MQ=255)
                                                                        aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                                                               >  2:2574013/1‑150 (MQ=255)
                                                                        aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                                                               >  1:925621/1‑150 (MQ=255)
                                                                         tCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCgg                                                                              >  1:63891/1‑150 (MQ=255)
                                                                           gTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGaa                                                                            <  2:2050816/150‑1 (MQ=255)
                                                                                cAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgt                                                                       >  1:1295277/1‑150 (MQ=255)
                                                                                 aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg                                                                      <  2:1314254/150‑1 (MQ=255)
                                                                                 aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg                                                                      >  1:4469778/1‑150 (MQ=255)
                                                                                   ccACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtggt                                                                    <  1:4051852/150‑1 (MQ=255)
                                                                                   ccACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtggt                                                                    <  2:5229757/150‑1 (MQ=255)
                                                                                   ccACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtggt                                                                    <  2:2789145/150‑1 (MQ=255)
                                                                                       cTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTc                                                                >  2:2525696/1‑150 (MQ=255)
                                                                                        tGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCa                                                               <  2:4627111/150‑1 (MQ=255)
                                                                                        tGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCa                                                               >  2:1232418/1‑150 (MQ=255)
                                                                                          cacTTCTGGGACAATATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGACGGCGGCAGCGAGCGCCCAGAGCTGGACTCCGATGCGGAAGCCGTGGTGTTCAtc                                                             <  1:2011720/148‑1 (MQ=255)
                                                                                           ccTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCa                                                                                                 <  1:768930/113‑1 (MQ=255)
                                                                                            cTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAtcgt                                                           >  1:4440923/1‑150 (MQ=255)
                                                                                             ttCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAg                                                                                                                                           >  1:2452294/1‑69 (MQ=255)
                                                                                             ttCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAtcgtc                                                          <  2:4070382/150‑1 (MQ=255)
                                                                                                tGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGaa                                                       <  2:401085/150‑1 (MQ=255)
                                                                                                tGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGaa                                                       <  1:375142/150‑1 (MQ=255)
                                                                                                   gACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCTGTGGTGTTCATCGTCGAAGGc                                                    <  2:1470426/150‑1 (MQ=255)
                                                                                                   gACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGc                                                    >  2:3418602/1‑150 (MQ=255)
                                                                                                       aGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAtc                                                             <  2:758818/137‑1 (MQ=255)
                                                                                                        gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGgtgt                                                                  <  1:2130965/131‑1 (MQ=255)
                                                                                                        gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGaagt                                               >  2:4209206/1‑148 (MQ=255)
                                                                                                        gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACt                                               >  1:5304111/1‑150 (MQ=255)
                                                                                                        gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACt                                               >  1:2577486/1‑150 (MQ=255)
                                                                                                           gcgcATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGa                                            <  2:4526266/150‑1 (MQ=255)
                                                                                                             gcctgggggTTATCGCCCGCCCGCGGACCGGGTGGGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACa                                          <  1:3944727/144‑1 (MQ=255)
                                                                                                               aTGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATc                                        >  1:3555614/1‑150 (MQ=255)
                                                                                                                  tgGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTGCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCAcc                                     >  2:4322796/1‑150 (MQ=255)
                                                                                                                    ggTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGt                                   <  2:4440923/150‑1 (MQ=255)
                                                                                                                     gTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTc                                  <  2:3475071/150‑1 (MQ=255)
                                                                                                                      ttATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCg                                 >  2:5266140/1‑150 (MQ=255)
                                                                                                                       tATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGa                                <  1:976250/150‑1 (MQ=255)
                                                                                                                        aTCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGaa                               >  1:1303056/1‑150 (MQ=255)
                                                                                                                          cgcccgccCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGTAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAgg                             <  2:4710504/150‑1 (MQ=255)
                                                                                                                                   gcTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAacacca                    >  1:139250/1‑150 (MQ=255)
                                                                                                                                   gcTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAacacca                    >  1:1095310/1‑150 (MQ=255)
                                                                                                                                    cTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAacaccac                   >  2:37513/1‑150 (MQ=255)
                                                                                                                                              ttGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAACACCACACCCTGGTAc         >  2:4693568/1‑150 (MQ=255)
                                                                                                                                                gCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAACACCACACCCTGGTACCt       <  2:1933032/150‑1 (MQ=255)
                                                                                                                                                  cGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAACACCACACCCTGGTACCtgg     >  1:3677606/1‑150 (MQ=255)
                                                                                                                                                   gAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAACACCACACCCTGGTACCtggt    <  1:4015401/150‑1 (MQ=255)
                                                                                                                                                    aaaCGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAACACCACACCCTGGTACCtggtg   <  2:3881331/150‑1 (MQ=255)
                                                                                                                                                    aaaCGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAACACCACACCCTGGTACCtggtg   <  1:3522754/150‑1 (MQ=255)
                                                                                                                                                    aaaCGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCACCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAACACCACACCCTGGTACCtggtg   >  2:1457591/1‑150 (MQ=255)
                                                                                                                                                     aaCGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAACACCACACCCTGGTACCtggtgg  <  2:5119242/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATCACCGTCGAAGGCAAACACCACACCCTGGTACCTGGTGG  >  NC_002947/4003831‑4004129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: