Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,195,990 T→C 10.5% 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,195,9900TC10.5% 202.6 / 12.3 76A175A (GCA→GCGybgJallophanate hydrolase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (30/38);  new base C (5/3);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

ACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCG  >  NC_002947/5195846‑5196138
                                                                                                                                                 |                                                                                                                                                    
accaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTCTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGcc                                                                                                                                                   <  4:51720/148‑1 (MQ=255)
  caccCATGGCCAGGAATTCGCTCCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCggg                                                                                                                                                <  4:288959/149‑1 (MQ=255)
   accCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGt                                                                                                                                               <  2:72684/149‑1 (MQ=255)
     ccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAg                                                                                                                                             <  3:304037/149‑1 (MQ=255)
     ccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAg                                                                                                                                             <  7:76643/149‑1 (MQ=255)
      cATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCAGCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAgg                                                                                                                                            <  6:194555/149‑1 (MQ=255)
      cATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAgg                                                                                                                                            <  6:21671/149‑1 (MQ=255)
      cATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCGCACGCTCACGGTCGCACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAgg                                                                                                                                            >  3:267312/1‑149 (MQ=255)
       aTGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCTACCCGGCAGTATCGCCGGGtggc                                                                                                                                            <  5:369577/149‑4 (MQ=255)
       aTGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgc                                                                                                                                           >  6:50643/1‑149 (MQ=255)
       aTGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGtggct                                                                                                                                           <  3:202741/149‑5 (MQ=255)
           ccAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCggg                                                                                                                                                <  5:37101/140‑1 (MQ=255)
           ccAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCggg                                                                                                                                                >  6:37101/1‑140 (MQ=255)
           ccAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgccgcc                                                                                                                                       <  6:245179/149‑1 (MQ=255)
            cAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGAAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGccgccg                                                                                                                                      >  1:362042/1‑149 (MQ=255)
               gAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGGCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGGGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCt                                                                                                                                   <  6:196277/149‑1 (MQ=255)
                  ttCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGgc                                                                                                                                <  2:288999/149‑1 (MQ=255)
                      cTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTc                                                                                                                            <  7:336219/149‑1 (MQ=255)
                         cgcgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCg                                                                                                                           <  6:230723/147‑1 (MQ=255)
                          gcgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGc                                                                                                                         >  8:101184/1‑148 (MQ=255)
                          gcgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGATCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGc                                                                                                                         >  6:249807/1‑148 (MQ=255)
                           cgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCCGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCaa                                                                                                                       >  5:310082/1‑149 (MQ=255)
                            gAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCt                                                                                                                                   <  8:179011/136‑1 (MQ=255)
                            gAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCt                                                                                                                                   >  7:179011/1‑136 (MQ=255)
                              aaCCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATAc                                                                                                                    <  2:83697/149‑1 (MQ=255)
                                  ggCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCAc                                                                                                                <  7:101184/149‑1 (MQ=255)
                                   gCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACg                                                                                                               <  5:124349/149‑1 (MQ=255)
                                      cGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGtg                                                                                                                                              >  1:288997/1‑114 (MQ=255)
                                      cGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGtg                                                                                                                                              <  2:288997/115‑2 (MQ=255)
                                            gcACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCccgcc                                                                                                      <  1:81472/149‑1 (MQ=255)
                                                ccGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCa                                                                                                  >  7:211122/1‑149 (MQ=255)
                                                   gTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCAcgc                                                                                               <  2:344734/149‑1 (MQ=255)
                                                      aCCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                          <  3:228335/119‑1 (MQ=255)
                                                      aCCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                          >  4:228335/1‑119 (MQ=255)
                                                          ggTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCg                                                                                        <  1:157081/149‑1 (MQ=255)
                                                           gTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGt                                                                                       <  6:315127/149‑1 (MQ=255)
                                                           gTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGt                                                                                       >  5:315127/1‑149 (MQ=255)
                                                                aaGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgggg                                                                                   >  1:27934/1‑148 (MQ=255)
                                                                         tAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGa                                                                           <  5:50643/147‑1 (MQ=255)
                                                                           aCCCTCACGCTTACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGa                                                                        <  3:1583/148‑1 (MQ=255)
                                                                            cccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAgg                                                                      <  7:257212/149‑1 (MQ=255)
                                                                             ccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAggg                                                                     >  5:125628/1‑149 (MQ=255)
                                                                             ccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGCCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAggg                                                                     <  5:249807/149‑1 (MQ=255)
                                                                              ctcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGc                                                                    >  6:193783/1‑149 (MQ=255)
                                                                                 acgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTgg                                                                 <  6:47569/149‑1 (MQ=255)
                                                                                   gctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGc                                                                >  5:13493/1‑148 (MQ=255)
                                                                                   gctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGc                                                                >  5:8106/1‑148 (MQ=255)
                                                                                             gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTggg                                                                                    >  7:74487/1‑118 (MQ=255)
                                                                                             gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTggg                                                                                    <  8:74487/118‑1 (MQ=255)
                                                                                             gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAggg                                                                     >  4:272142/1‑133 (MQ=255)
                                                                                             gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgt                                                     <  2:27934/149‑1 (MQ=255)
                                                                                             gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgt                                                     >  7:129987/1‑149 (MQ=255)
                                                                                               acaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCACCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtcg                                                   >  3:74352/1‑149 (MQ=255)
                                                                                                 aacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGAc                                                 <  2:362042/149‑1 (MQ=255)
                                                                                                     gCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAgg                                             >  2:175783/1‑149 (MQ=255)
                                                                                                       aCCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGcc                                           >  7:179307/1‑149 (MQ=255)
                                                                                                       aCCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGcc                                           >  4:37097/1‑149 (MQ=255)
                                                                                                       aCCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTAGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGca                                           >  4:209049/1‑148 (MQ=255)
                                                                                                              tGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCACCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAg                                    <  4:74352/149‑1 (MQ=255)
                                                                                                                    cGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGa                                                  <  3:274559/129‑1 (MQ=255)
                                                                                                                    cGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGa                                                  >  4:274559/1‑129 (MQ=255)
                                                                                                                    cGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGaaa                              <  6:8106/149‑1 (MQ=255)
                                                                                                                         aGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCgg                         <  6:13493/149‑1 (MQ=255)
                                                                                                                          ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                        <  8:179307/149‑1 (MQ=255)
                                                                                                                           gTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTg                       >  3:291013/1‑149 (MQ=255)
                                                                                                                           gTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTg                       >  6:11036/1‑149 (MQ=255)
                                                                                                                               cAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCg                     >  4:41602/1‑147 (MQ=255)
                                                                                                                                ggcctaCCCGGCAGTATCGCCGGGTGGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCg                     >  7:248527/6‑147 (MQ=255)
                                                                                                                                 gCCCCCCTGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATCCCCACGCTTCCCGCCTCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGc                 <  4:267312/149‑1 (MQ=255)
                                                                                                                                 gCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCCCGCTGCCCGCCGCCACGCGCTTGCGGGGGGTACTCAGACGAGCGCTGGCCAGCCGCTCGTCGACCAGGCCCCTGAATCCGAAACCCGCTGCGAAGc                 >  3:19427/1‑149 (MQ=255)
                                                                                                                                  cctaCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGc                 >  3:214923/4‑149 (MQ=255)
                                                                                                                                   ctaCCCGGCAGTATCGCCGGGTGGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGcc                >  2:354798/3‑149 (MQ=255)
                                                                                                                                    taCCCGGCAGTATCGCCGGGTGGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGcc                >  2:339148/2‑148 (MQ=255)
                                                                                                                                       ccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAgg            <  6:125628/149‑1 (MQ=255)
                                                                                                                                            gATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCTAAGCCCAGGGCaaa       >  1:229151/1‑149 (MQ=255)
                                                                                                                                             aTACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAac      <  5:11036/149‑1 (MQ=255)
                                                                                                                                                 tGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCg  >  6:170541/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAACACCG  >  NC_002947/5195846‑5196138

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