Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,003,937 T→C 6.2% M52T (ATG→ACG)  allE → S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,003,9370TC 6.2% 657.3 / 19.6 209M52T (ATG→ACG) allES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (96/100);  new base C (8/5);  total (104/105)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 209 total aligned reads displayed.

AACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTG  >  NC_002947/4003789‑4004085
                                                                                                                                                    |                                                                                                                                                    
   cATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATgcgc                                                                                                                                                  <  1:4719703/150‑1 (MQ=255)
       cATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgt                                                                                                                                              <  2:3911073/150‑1 (MQ=255)
        aTACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgtg                                                                                                                                             >  1:2789145/1‑150 (MQ=255)
        aTACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtgtg                                                                                                                                             >  2:3054397/1‑150 (MQ=255)
        aTACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCAtggg                                                                                                                                             >  2:4015401/1‑148 (MQ=255)
         tACTTCGCCCCCCACGGTGGGCACCAGGCTCAAACCGAGATGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGACGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTgg                                                                                                                                            >  2:2011720/1‑150 (MQ=255)
          aCTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTggg                                                                                                                                           <  1:3054397/150‑1 (MQ=255)
          aCTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTggg                                                                                                                                           <  2:1271128/150‑1 (MQ=255)
           cTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGt                                                                                                                                          >  1:4893990/1‑150 (MQ=255)
            ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                                                         >  2:5064745/1‑150 (MQ=255)
            ttCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGtt                                                                                                                                         >  1:4627158/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                                                                                                                                         >  2:3836936/1‑150 (MQ=255)
             tCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTa                                                                                                                                        <  2:2533162/150‑1 (MQ=255)
               gCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATc                                                                                                                                      <  2:3297606/150‑1 (MQ=255)
                 cccccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgc                                                                                                                                    <  1:1606906/150‑1 (MQ=255)
                 cccccACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgc                                                                                                                                    <  1:1582951/150‑1 (MQ=255)
                     cACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATcgcccgc                                                                                                                                >  2:3926111/1‑150 (MQ=255)
                               aCCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCg                                                                                                                      >  2:3681602/1‑150 (MQ=255)
                                 ccGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGt                                                                                                                    >  1:1359778/1‑150 (MQ=255)
                                      tCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGc                                                                                                               >  1:4526266/1‑150 (MQ=255)
                                       cAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGcc                                                                                                              <  1:828689/150‑1 (MQ=255)
                                        aaaCCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCg                                                                                                             >  1:1283513/1‑150 (MQ=255)
                                            cGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAAc                                                                                                         <  1:291817/150‑1 (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)
                                                  gctgACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTtctc                                                                                                   >  1:2726046/1‑150 (MQ=255)
                                                   ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                                                                  <  1:3341939/150‑1 (MQ=255)
                                                   ctgactgaCCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTcc                                                                                                  <  1:4612988/150‑1 (MQ=255)
                                                    tgaccgatCGCGCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTtctcgc                                                                                                 >  2:174779/1‑148 (MQ=255)
                                                     gactgacCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCa                                                                                                <  2:3153213/150‑1 (MQ=255)
                                                           ccGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAt                                                                                          <  2:3115479/150‑1 (MQ=255)
                                                            cgcgCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatatt                                                                                         <  1:1887276/147‑2 (MQ=255)
                                                            cGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtc                                                                                         <  2:3612626/150‑1 (MQ=255)
                                                            cGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtc                                                                                         >  2:3831530/1‑150 (MQ=255)
                                                               gCCATGTTCACCGAAGCCTATGCGGTTATTCCGAAAGGCACATTCAGCGATATCGTCACGAGCTTTCTTCCCTTCTGGGACAAGACGCGGCTGTGGATCATCGCCCGCCCGCTATCAGGTTTTGCCGAGACTTTCTCGCAatatattatg                                                                                      <  2:3466720/150‑5 (MQ=255)
                                                               gCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtcgtc                                                                                      <  2:1931985/150‑1 (MQ=255)
                                                               gCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACAtcgtc                                                                                      <  1:3871359/150‑1 (MQ=255)
                                                                   tGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTGCATCGTCGAAc                                                                                  >  1:2627360/1‑150 (MQ=255)
                                                                     ttCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTg                                                                                >  2:5369028/1‑150 (MQ=255)
                                                                        aCCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGcc                                                                             >  1:4821180/1‑150 (MQ=255)
                                                                          cGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGAGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGcccc                                                                           <  1:4461369/150‑1 (MQ=255)
                                                                            aaGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTggg                                                                                                                                           >  1:4283542/1‑84 (MQ=255)
                                                                            aaGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTggg                                                                                                                                           <  2:4283542/84‑1 (MQ=255)
                                                                              gCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGaa                                                                       >  1:996968/1‑150 (MQ=255)
                                                                              gCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGaa                                                                       <  2:5240346/150‑1 (MQ=255)
                                                                                 tATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAggc                                                                    >  2:2835337/1‑150 (MQ=255)
                                                                                    gCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAggcggc                                                                 <  2:5024543/150‑1 (MQ=255)
                                                                                      cGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCag                                                               >  2:1354790/1‑150 (MQ=255)
                                                                                       gTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagc                                                              >  2:3035391/1‑150 (MQ=255)
                                                                                         cATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagcga                                                            <  1:5227830/150‑1 (MQ=255)
                                                                                          aTCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCa                                                                                                >  2:635896/1‑113 (MQ=255)
                                                                                            ccccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCagcgagcg                                                         <  2:4821180/150‑1 (MQ=255)
                                                                                             cccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAgcgc                                                        <  1:3911510/150‑1 (MQ=255)
                                                                                              ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                                                       >  1:1244777/1‑150 (MQ=255)
                                                                                              ccAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGcc                                                       >  2:1086081/1‑150 (MQ=255)
                                                                                               cAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGccc                                                      >  1:5017976/1‑150 (MQ=255)
                                                                                               cAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGccc                                                      >  1:1421481/1‑150 (MQ=255)
                                                                                                     cGTGATGCGTGACATCGTCACCAGCCACCTGACCTTCTGGGACAAGATGCGCATGTGGGTTATCTCCCGCCCGCTGACCGGTTTTTCCGAAACGTTCTCCCAGTACATATTCTAAATGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCt                                                <  1:4397717/150‑1 (MQ=255)
                                                                                                       tGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTgg                                              >  1:5322214/1‑150 (MQ=255)
                                                                                                           gCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGc                                          <  2:3993133/150‑1 (MQ=255)
                                                                                                            cGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGcc                                         >  1:2044599/1‑150 (MQ=255)
                                                                                                               gACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGAt                                      >  2:2614330/1‑150 (MQ=255)
                                                                                                                 cATCGGCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGCGGGGTATCGCCCTCCCGCTGACCGGTTTTGCCGAACCGGTCTCCCAGGACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGc                                    <  2:214068/150‑1 (MQ=255)
                                                                                                                  aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                   >  2:2574013/1‑150 (MQ=255)
                                                                                                                  aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                   >  1:3431845/1‑150 (MQ=255)
                                                                                                                  aTCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCg                                   >  1:925621/1‑150 (MQ=255)
                                                                                                                   tCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCgg                                  >  1:63891/1‑150 (MQ=255)
                                                                                                                     gTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGaa                                >  2:1974657/1‑150 (MQ=255)
                                                                                                                     gTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGaa                                <  2:2050816/150‑1 (MQ=255)
                                                                                                                     gTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGaa                                <  1:5102929/150‑1 (MQ=255)
                                                                                                                          cAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgt                           >  1:1295277/1‑150 (MQ=255)
                                                                                                                           aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGACCTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg                          <  1:3609734/150‑1 (MQ=255)
                                                                                                                           aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg                          <  2:3322695/150‑1 (MQ=255)
                                                                                                                           aGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtg                          <  2:1314254/150‑1 (MQ=255)
                                                                                                                            gCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtgg                         <  2:3059872/150‑1 (MQ=255)
                                                                                                                            gCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtgg                         <  2:4451819/150‑1 (MQ=255)
                                                                                                                             ccACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtggt                        <  2:5229757/150‑1 (MQ=255)
                                                                                                                             ccACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtggt                        <  1:4051852/150‑1 (MQ=255)
                                                                                                                             ccACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCgtggt                        <  2:2789145/150‑1 (MQ=255)
                                                                                                                               aCCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGgtgt                      >  2:84324/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)
                                                                                                                                    cccTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAtc                 <  1:2997624/150‑1 (MQ=255)
                                                                                                                                     ccTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCa                                                     >  2:768930/1‑113 (MQ=255)
                                                                                                                                       ttCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAg                                                                                               >  1:2452294/1‑69 (MQ=255)
                                                                                                                                       ttCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAg                                                                                               <  2:2452294/69‑1 (MQ=255)
                                                                                                                                       ttCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAtcgtc              <  2:4070382/150‑1 (MQ=255)
                                                                                                                                          tGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGaa           >  1:401085/1‑150 (MQ=255)
                                                                                                                                           gggACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAg          <  2:3110517/150‑1 (MQ=255)
                                                                                                                                             gACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCTGTGGTGTTCATCGTCGAAGGc        <  2:1470426/150‑1 (MQ=255)
                                                                                                                                                 aGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCAtc                 >  1:758818/1‑137 (MQ=255)
                                                                                                                                                  gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGgtgt                      <  1:2130965/131‑1 (MQ=255)
                                                                                                                                                  gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGaagt   >  2:4209206/1‑148 (MQ=255)
                                                                                                                                                  gATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACt   >  1:2577486/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGGGTTATCGCCCGCCCGCTGACCGGTTTTGCCGAAACGTTCTCCCAGTACATCGTCGAACTGGCCCCGGAAGGCGGCAGCGAGCGCCCAGAGCTGGACGCCGATGCGGAAGCCGTGGTGTTCATCGTCGAAGGCGAACTG  >  NC_002947/4003789‑4004085

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: