Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,200,987 T→C 5.2% A175A (GCA→GCG ybgJ ← allophanate hydrolase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,200,9870TC 5.2% 1485.0 / 39.8 446A175A (GCA→GCGybgJallophanate hydrolase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (189/233);  new base C (11/12);  total (200/246)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 457 total aligned reads displayed.

ACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCGG  >  NC_002947/5200838‑5201136
                                                                                                                                                     |                                                                                                                                                     
acaTCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACt                                                                                                                                                       >  1:3023378/1‑150 (MQ=255)
     accaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCgg                                                                                                                                                  >  1:6215956/1‑150 (MQ=255)
      ccaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCggg                                                                                                                                                 >  2:6058249/1‑150 (MQ=255)
       caccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGt                                                                                                                                                <  1:3153376/150‑1 (MQ=255)
        accCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTa                                                                                                                                               <  2:4911150/150‑1 (MQ=255)
        accCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTa                                                                                                                                               <  1:6997060/150‑1 (MQ=255)
        accCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTa                                                                                                                                               >  1:3293090/1‑150 (MQ=255)
         cccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAg                                                                                                                                              <  1:2822487/150‑1 (MQ=255)
          ccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAgg                                                                                                                                             <  1:2358151/150‑1 (MQ=255)
           cctGGCCAGGAATTCGCTGCGCGAAACCGGCACGCAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGtggct                                                                                                                                            <  2:321636/148‑5 (MQ=255)
            ctGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgcc                                                                                                                                           <  2:5487810/149‑1 (MQ=255)
            aTGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAg                                                                                                                                              >  1:3750707/1‑147 (MQ=255)
              ggCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgccgc                                                                                                                                         <  1:7428913/150‑1 (MQ=255)
              ggCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgccgc                                                                                                                                         >  1:3618765/1‑150 (MQ=255)
               gCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgccgcc                                                                                                                                        <  1:901846/150‑1 (MQ=255)
                ccAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGccgccg                                                                                                                                       >  2:333931/1‑150 (MQ=255)
                ccAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGccgccg                                                                                                                                       >  1:2171197/1‑150 (MQ=255)
                  aGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTc                                                                                                                                     <  1:1325035/150‑1 (MQ=255)
                  aGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTc                                                                                                                                     <  1:5246265/150‑1 (MQ=255)
                    gAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTg                                                                                                                                   <  2:4565582/150‑1 (MQ=255)
                      aTTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGgc                                                                                                                                 <  2:5091811/150‑1 (MQ=255)
                       ttCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGgcg                                                                                                                                <  2:7105888/150‑1 (MQ=255)
                            tGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                           >  2:3500088/1‑150 (MQ=255)
                            tGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                           <  2:2098828/150‑1 (MQ=255)
                            ggcgcGAAACCGGCACGAACCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCCATCAGGTTCCAGCCACCCGGCGATACAGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                           <  2:2645484/149‑1 (MQ=255)
                                   aaCCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATAcc                                                                                                                    <  1:4166573/150‑1 (MQ=255)
                                    aCCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATAccc                                                                                                                   >  2:5802974/1‑150 (MQ=255)
                                        gCACGAAGCGCACCCGGTCACCGTGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCTGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGACGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGc                                                                                                               <  1:3103014/150‑1 (MQ=255)
                                            gAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGcc                                                                                                           >  2:3854161/1‑150 (MQ=255)
                                                   aCCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCgccgcc                                                                                                    >  1:2737763/1‑150 (MQ=255)
                                                       ggTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCAcgc                                                                                                <  1:2015125/150‑1 (MQ=255)
                                                           aCCGGGTTGCAAGCGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCtt                                                                                            <  2:5414488/150‑1 (MQ=255)
                                                                gTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTg                                                                                       <  1:4264502/150‑1 (MQ=255)
                                                                gTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTg                                                                                       <  2:4781867/150‑1 (MQ=255)
                                                                gTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTg                                                                                       >  2:6539495/1‑150 (MQ=255)
                                                                  tGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTggg                                                                                     <  2:843427/150‑1 (MQ=255)
                                                                  tGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTggg                                                                                     <  1:3522670/150‑1 (MQ=255)
                                                                  tGCAAGAGGCTGGAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTggg                                                                                     <  2:5238801/150‑1 (MQ=255)
                                                                   gCAAGAGGCTGTAACCCTCCCGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgggg                                                                                    <  2:4258207/150‑1 (MQ=255)
                                                                    cAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGt                                                                                   >  1:968018/1‑150 (MQ=255)
                                                                      agagGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTAc                                                                                 >  2:5677146/1‑150 (MQ=255)
                                                                      agagGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTAc                                                                                 <  2:1847288/150‑1 (MQ=255)
                                                                            tGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCACTACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTCCTCAGAc                                                                           >  2:5153482/1‑150 (MQ=255)
                                                                            tGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGAc                                                                           >  2:2309836/1‑150 (MQ=255)
                                                                              tAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGa                                                                         <  1:2152713/150‑1 (MQ=255)
                                                                                  ccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGc                                                                     >  1:5012448/1‑150 (MQ=255)
                                                                                  ccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGc                                                                     >  1:545762/1‑150 (MQ=255)
                                                                                    tcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTg                                                                   <  2:1269012/150‑1 (MQ=255)
                                                                                    tcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTg                                                                   >  2:6258139/1‑150 (MQ=255)
                                                                                    tcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTg                                                                   >  1:724620/1‑150 (MQ=255)
                                                                                     cacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTgg                                                                  <  2:1965811/150‑1 (MQ=255)
                                                                                        gctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCa                                                               >  1:6072180/1‑150 (MQ=255)
                                                                                         ctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAg                                                              <  2:5930932/150‑1 (MQ=255)
                                                                                          tcacgGTCGAACAAACGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAgg                                                             >  2:4645593/1‑150 (MQ=255)
                                                                                               gTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtc                                                        >  1:3414537/1‑150 (MQ=255)
                                                                                               gTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtc                                                        <  2:2782827/150‑1 (MQ=255)
                                                                                                tCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcg                                                       >  1:6978807/1‑150 (MQ=255)
                                                                                                  gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                     <  1:837640/150‑1 (MQ=255)
                                                                                                  gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                     >  2:2953479/1‑150 (MQ=255)
                                                                                                    acaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTAGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGa                                                   >  1:1517775/1‑150 (MQ=255)
                                                                                                     cgggcGCACCCCGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtc                                                        <  2:7358807/140‑1 (MQ=255)
                                                                                                     caacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGAc                                                  <  1:5731462/150‑1 (MQ=255)
                                                                                                       acCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATTCCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCa                                                <  1:6747073/150‑1 (MQ=255)
                                                                                                           cACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgg                                                                                      >  2:1427341/1‑108 (MQ=255)
                                                                                                           cACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgg                                                                                      <  1:1427341/108‑1 (MQ=255)
                                                                                                           cACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGcc                                            >  1:1007889/1‑150 (MQ=255)
                                                                                                            aCCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGccc                                           >  1:4808239/1‑150 (MQ=255)
                                                                                                             ccGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCa                                          <  1:3989676/150‑1 (MQ=255)
                                                                                                                 ggTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGaa                                      <  1:1798227/150‑1 (MQ=255)
                                                                                                                 ggTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGaa                                      <  1:2576672/150‑1 (MQ=255)
                                                                                                                          gATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAAcc                             >  2:6942963/1‑150 (MQ=255)
                                                                                                                            tCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCg                           <  2:379453/150‑1 (MQ=255)
                                                                                                                             cAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCgg                          <  2:35575/150‑1 (MQ=255)
                                                                                                                              aGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                         >  2:4188258/1‑150 (MQ=255)
                                                                                                                               ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTg                        <  1:4888875/150‑1 (MQ=255)
                                                                                                                                   ccAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGaa                    <  2:6784478/150‑1 (MQ=255)
                                                                                                                                   ccAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGaa                    <  1:3924111/150‑1 (MQ=255)
                                                                                                                                   ccAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGaa                    <  1:447350/150‑1 (MQ=255)
                                                                                                                                     ggCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGc                  >  1:5834680/2‑150 (MQ=255)
                                                                                                                                      gCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGcc                 >  1:6473188/1‑150 (MQ=255)
                                                                                                                                      gCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGcc                 >  1:4716476/1‑150 (MQ=255)
                                                                                                                                      cctaCCCGGCAGTATCGCCGGGTGGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTCCTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGcc                 >  2:2577336/4‑150 (MQ=255)
                                                                                                                                       ctaCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGccc                >  2:1099199/3‑150 (MQ=255)
                                                                                                                                       ccACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGccc                <  1:7059057/150‑1 (MQ=255)
                                                                                                                                        taCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCa               >  1:4666555/2‑150 (MQ=255)
                                                                                                                                        taCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCa               >  1:5641924/2‑150 (MQ=255)
                                                                                                                                          cccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAgg             <  1:2309642/150‑1 (MQ=255)
                                                                                                                                           ccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGa                                                                         >  1:6599270/1‑89 (MQ=255)
                                                                                                                                           ccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAggg            >  1:4231375/1‑150 (MQ=255)
                                                                                                                                           ccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAggg            <  2:563237/150‑1 (MQ=255)
                                                                                                                                             ggCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAg                   >  2:1616215/1‑141 (MQ=255)
                                                                                                                                                gATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAac       >  1:4539743/1‑150 (MQ=255)
                                                                                                                                                   aCTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCaa         <  1:2470220/145‑1 (MQ=255)
                                                                                                                                                   aCTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCaa         >  2:2470220/1‑145 (MQ=255)
                                                                                                                                                    cTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCg   <  1:3791920/150‑1 (MQ=255)
                                                                                                                                                    cTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCg   <  2:3662562/150‑1 (MQ=255)
                                                                                                                                                     tGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCgg  >  2:4020567/1‑150 (MQ=255)
                                                                                                                                                     tGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCgg  <  1:1633613/150‑1 (MQ=255)
                                                                                                                                                     tGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCgg  <  1:1175815/150‑1 (MQ=255)
                                                                                                                                                     tGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTACCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCgg  <  2:2171170/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCGG  >  NC_002947/5200838‑5201136

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: