Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,003,942 A→C 6.1% M54L (ATG→CTG)  allE → S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,003,9420AC 6.1% 688.0 / 23.7 207M54L (ATG→CTG) allES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (92/102);  new base C (8/5);  total (100/107)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 211 total aligned reads displayed.

ATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATC  >  NC_002947/4003793‑4004091
                                                                                                                                                     |                                                                                                                                                     
aTTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCa                                                                                                                                                       >  2:375142/1‑150 (MQ=255)
 ttCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAt                                                                                                                                                      <  2:2018834/150‑1 (MQ=255)
 ttCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAt                                                                                                                                                      >  2:1606906/1‑150 (MQ=255)
   cATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgt                                                                                                                                                    <  1:210409/150‑1 (MQ=255)
    aTACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgtg                                                                                                                                                   >  2:2827515/1‑150 (MQ=255)
    aTACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgtg                                                                                                                                                   >  2:3054301/1‑150 (MQ=255)
    aTACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtggg                                                                                                                                                   >  2:4015401/1‑148 (MQ=255)
      aCTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTggg                                                                                                                                                 <  2:1271128/150‑1 (MQ=255)
      aCTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTggg                                                                                                                                                 <  1:3054397/150‑1 (MQ=255)
       cTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGt                                                                                                                                                >  1:4893990/1‑150 (MQ=255)
        ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                                                               >  2:1987956/1‑150 (MQ=255)
        ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                                                               <  1:1555934/150‑1 (MQ=255)
        ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                                                               >  1:841712/1‑150 (MQ=255)
        ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                                                               >  1:3153213/1‑150 (MQ=255)
        ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                                                               >  1:4627158/1‑150 (MQ=255)
           gCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATc                                                                                                                                            <  2:3297606/150‑1 (MQ=255)
             cccccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgc                                                                                                                                          <  1:1582951/150‑1 (MQ=255)
             cccccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgc                                                                                                                                          <  1:1606906/150‑1 (MQ=255)
               cccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgccc                                                                                                                                        <  1:427785/150‑1 (MQ=255)
                 cACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgcccgc                                                                                                                                      >  2:3926111/1‑150 (MQ=255)
                        ggCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGa                                                                                                                               <  2:4503876/150‑1 (MQ=255)
                           aCCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCg                                                                                                                            >  2:3681589/1‑150 (MQ=255)
                             ccGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGt                                                                                                                          >  1:1359778/1‑150 (MQ=255)
                              cGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGtt                                                                                                                         <  2:4449070/150‑1 (MQ=255)
                               ggCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGttt                                                                                                                        >  1:3475071/1‑150 (MQ=255)
                                  tCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGc                                                                                                                     <  1:4371376/150‑1 (MQ=255)
                                   cAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGcc                                                                                                                    <  1:828689/150‑1 (MQ=255)
                                    aaaCCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCg                                                                                                                   >  1:1283513/1‑150 (MQ=255)
                                        cGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAAc                                                                                                               >  1:1060396/1‑150 (MQ=255)
                                           gctgctGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGtt                                                                                                            <  1:590523/150‑1 (MQ=255)
                                            ctcctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtc                                                                                                           <  1:3749418/147‑1 (MQ=255)
                                             tcctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtct                                                                                                          >  1:861330/3‑150 (MQ=255)
                                               ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                                                                        <  1:2827515/150‑1 (MQ=255)
                                               ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                                                                        <  1:4612988/150‑1 (MQ=255)
                                               ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                                                                        <  1:3341939/150‑1 (MQ=255)
                                               ctgaccgaTCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtctcg                                                                                                        >  2:3070560/1‑149 (MQ=255)
                                                tgaccgatCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtctcgc                                                                                                       >  2:175181/1‑148 (MQ=255)
                                                 gactgacCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCa                                                                                                      <  2:3153213/150‑1 (MQ=255)
                                                       ccGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAt                                                                                                <  2:3115458/150‑1 (MQ=255)
                                                        cgcgCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatatt                                                                                               <  1:1887276/147‑2 (MQ=255)
                                                        cGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCACAGTACAtc                                                                                               >  2:3831942/1‑150 (MQ=255)
                                                          tGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtcgt                                                                                             <  1:4559789/150‑1 (MQ=255)
                                                           gCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatg                                                                                            <  2:3466720/150‑5 (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:1932025/150‑1 (MQ=255)
                                                            ccATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatgg                                                                                           <  1:2254002/150‑6 (MQ=255)
                                                               tGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTGCATCGTCGAAc                                                                                        >  1:2627360/1‑150 (MQ=255)
                                                                 ttCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTg                                                                                      >  2:5369028/1‑150 (MQ=255)
                                                                    aCCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGcc                                                                                   >  1:4821180/1‑150 (MQ=255)
                                                                     ccGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGccc                                                                                  <  2:4902127/150‑1 (MQ=255)
                                                                        aaGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTggg                                                                                                                                                 <  2:4283542/84‑1 (MQ=255)
                                                                          gCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGaa                                                                             >  1:996968/1‑150 (MQ=255)
                                                                          gCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGaa                                                                             <  2:5240346/150‑1 (MQ=255)
                                                                             tATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAggc                                                                          <  1:2932358/150‑1 (MQ=255)
                                                                                   gTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagc                                                                    >  2:2407750/1‑150 (MQ=255)
                                                                                     cATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagcga                                                                  <  1:5227830/150‑1 (MQ=255)
                                                                                      aTCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCa                                                                                                      >  2:635896/1‑113 (MQ=255)
                                                                                       tCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagcgagc                                                                >  2:5038068/1‑150 (MQ=255)
                                                                                          ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                                                             >  2:3021254/1‑150 (MQ=255)
                                                                                          ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                                                             >  1:1244777/1‑150 (MQ=255)
                                                                                          ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                                                             >  2:1086081/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)
                                                                                                   tGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTgg                                                    >  1:5322214/1‑150 (MQ=255)
                                                                                                       gCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGc                                                <  2:3993133/150‑1 (MQ=255)
                                                                                                       gCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGc                                                <  1:2735671/150‑1 (MQ=255)
                                                                                                           gACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGAt                                            >  2:2614330/1‑150 (MQ=255)
                                                                                                            aCATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATg                                           <  2:3845615/150‑1 (MQ=255)
                                                                                                             cATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGc                                          <  1:3942904/150‑1 (MQ=255)
                                                                                                              aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                         >  1:3431845/1‑150 (MQ=255)
                                                                                                              aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                         >  2:2574013/1‑150 (MQ=255)
                                                                                                              aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                         >  2:2642887/1‑150 (MQ=255)
                                                                                                              aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                         >  1:925621/1‑150 (MQ=255)
                                                                                                               tCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCgg                                        >  1:63891/1‑150 (MQ=255)
                                                                                                                 gTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGaa                                      <  1:5102929/150‑1 (MQ=255)
                                                                                                                      cAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgt                                 >  1:1295277/1‑150 (MQ=255)
                                                                                                                       aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGACCTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg                                <  1:3609734/150‑1 (MQ=255)
                                                                                                                       aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg                                >  1:4469778/1‑150 (MQ=255)
                                                                                                                       aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg                                <  2:3322695/150‑1 (MQ=255)
                                                                                                                       aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg                                <  2:1314254/150‑1 (MQ=255)
                                                                                                                         ccACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtggt                              <  2:5229757/150‑1 (MQ=255)
                                                                                                                         ccACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACAGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtggt                              <  1:4051780/150‑1 (MQ=255)
                                                                                                                           aCCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGgtgt                            >  2:84324/1‑150 (MQ=255)
                                                                                                                             cTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTc                          >  2:2525696/1‑150 (MQ=255)
                                                                                                                              tGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCa                         >  1:3571439/1‑150 (MQ=255)
                                                                                                                              tGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCa                         <  2:4627111/150‑1 (MQ=255)
                                                                                                                              tGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCa                         >  2:1232418/1‑150 (MQ=255)
                                                                                                                                 ccTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCa                                                           >  2:768930/1‑113 (MQ=255)
                                                                                                                                  cTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAtcgt                     >  1:4440923/1‑150 (MQ=255)
                                                                                                                                   ttCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAg                                                                                                     >  1:2452294/1‑69 (MQ=255)
                                                                                                                                   ttCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAtcgtc                    <  2:4070382/150‑1 (MQ=255)
                                                                                                                                              gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGgtgt                            >  2:2130965/1‑131 (MQ=255)
                                                                                                                                              gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGaagt         >  2:4209206/1‑148 (MQ=255)
                                                                                                                                              gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACt         >  1:2577486/1‑150 (MQ=255)
                                                                                                                                              gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACt         >  2:1462324/1‑150 (MQ=255)
                                                                                                                                                 gcgcATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGa      <  2:5193642/150‑1 (MQ=255)
                                                                                                                                                 gcgcAGGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGa      <  1:2642887/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTGGACATC  >  NC_002947/4003793‑4004091

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: