Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,456,994 A→C 8.0% E22A (GAG→GCG)  PP_2150 → thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,456,9940AC 8.0% 372.4 / 18.3 122E22A (GAG→GCG) PP_2150thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (64/48);  new base C (4/6);  total (68/54)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

CCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTCC  >  NC_002947/2456846‑2457142
                                                                                                                                                      |                                                                                                                                                    
cccGCTCCTACAGGATTATCGCTGCCCCTTTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGa                                                                                                                                                      <  1:149038/149‑1 (MQ=255)
cccGCTCCTACAGGATTATCGCTGCCCCTGTATGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCATCCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGa                                                                                                                                                      <  7:233811/149‑1 (MQ=255)
 ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCACCCAGGGCCCCACCGTGGAg                                                                                                                                                     >  8:281356/1‑149 (MQ=255)
       cTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCAGGCGCTccag                                                                                                                                               >  4:26668/1‑147 (MQ=255)
        tACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                              >  4:216283/1‑149 (MQ=255)
          cAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGggt                                                                                                                                            >  1:59541/1‑149 (MQ=255)
           aGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtg                                                                                                                                           <  3:346780/149‑1 (MQ=255)
           aGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtg                                                                                                                                           >  6:14169/1‑149 (MQ=255)
            gtaTTATCGCGCCCCATGTCGTAGGGGTTTTACCCGCTAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCACACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtgg                                                                                                                                          <  6:10017/147‑1 (MQ=255)
            ggATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                              >  7:153656/1‑145 (MQ=255)
            ggATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                              <  8:153656/145‑1 (MQ=255)
             gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtggt                                                                                                                                         <  4:346091/149‑1 (MQ=255)
             gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtggt                                                                                                                                         >  3:119946/1‑149 (MQ=255)
              aTTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtgggg                                                                                                                                        <  3:26668/149‑3 (MQ=255)
               ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTcc                                                                                                                                       >  3:77866/1‑149 (MQ=255)
                  tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtgggatagg                                                                                                                                    >  1:128796/1‑143 (MQ=255)
                    gCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                              >  2:164346/1‑137 (MQ=255)
                    gCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                              <  1:164346/137‑1 (MQ=255)
                     cTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCAt                                                                                                                                 >  1:214549/1‑149 (MQ=255)
                       gCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATgg                                                                                                                               >  4:174764/1‑149 (MQ=255)
                        ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                              <  5:14169/149‑1 (MQ=255)
                          ccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                              >  3:175501/1‑147 (MQ=255)
                          ccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                              >  7:133554/1‑147 (MQ=255)
                           cTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcag                                                                                                                           <  3:216283/149‑1 (MQ=255)
                                ggAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCagcagc                                                                                                                       <  8:283623/148‑1 (MQ=255)
                                 gAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                             <  4:175521/125‑1 (MQ=255)
                                 gAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                             >  3:175521/1‑125 (MQ=255)
                                  aGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTAc                                                                                                                    <  7:280276/149‑1 (MQ=255)
                                  aGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTAc                                                                                                                    <  6:256921/149‑1 (MQ=255)
                                     gggTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                             >  1:45569/1‑121 (MQ=255)
                                     gggTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                             <  2:45569/121‑1 (MQ=255)
                                      ggTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGTAGCTACAGCa                                                                                                                >  7:106376/1‑149 (MQ=255)
                                      ggTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCa                                                                                                                <  1:292529/149‑1 (MQ=255)
                                            cccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAg                                                                                                           <  7:69209/148‑1 (MQ=255)
                                            cccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAg                                                                                                           >  6:277331/1‑148 (MQ=255)
                                            cccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGt                                                                                                          <  4:119946/149‑1 (MQ=255)
                                             ccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagc                                                                                                                          <  1:282166/132‑1 (MQ=255)
                                             ccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagc                                                                                                                          >  2:282166/1‑132 (MQ=255)
                                             ccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTa                                                                                                         >  1:26231/1‑149 (MQ=255)
                                                     agGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCcgc                                                                                                  >  2:6146/1‑148 (MQ=255)
                                                       gCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAg                                                                                               >  2:342175/1‑149 (MQ=255)
                                                         cAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGcc                                                                                             <  6:335751/149‑1 (MQ=255)
                                                                 gCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGcc                                                                                     >  8:112348/1‑149 (MQ=255)
                                                                 gCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGCGGcc                                                                                     >  7:337434/1‑149 (MQ=255)
                                                                 gCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCACCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGACTCCAGTACGTCCGCGAGCCTGGTGGcc                                                                                     >  3:108887/1‑149 (MQ=255)
                                                                     ccGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCt                                                                                 <  8:18700/149‑1 (MQ=255)
                                                                     ccGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCt                                                                                 >  6:128191/1‑149 (MQ=255)
                                                                      cGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCGGCTGAAACCAAGCCCCCCCCGTGGAGCGCATGGGTGGACCGCCAATGGCCACCAGAAACAGAAACAATAACGACCGCGCCCATGGTGCCCCGc                                                                                  >  1:98787/1‑149 (MQ=255)
                                                                            aCCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGcc                                                                          <  4:77866/149‑1 (MQ=255)
                                                                             cccTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGc                                                                           <  8:215408/147‑1 (MQ=255)
                                                                               cTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGTAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAg                                                                       >  4:293743/1‑149 (MQ=255)
                                                                               cTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAg                                                                       <  3:174764/149‑1 (MQ=255)
                                                                               cTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTAGCCCGCTGCCGGCCCAg                                                                       >  1:123108/1‑149 (MQ=255)
                                                                                tACTGTGCGCACAGCCCTACTGTACATTCTGCTGCTACTCACCGC‑‑CTCCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCACCAGCTACAGCAACCAGAACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAgg                                                                      >  7:65444/1‑149 (MQ=255)
                                                                                    gtgCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCt                                                                                 <  8:189590/134‑1 (MQ=255)
                                                                                    gtgCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCt                                                                                 >  7:189590/1‑134 (MQ=255)
                                                                                     tgCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCAACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGc                                                                   <  7:281356/147‑1 (MQ=255)
                                                                                        gcACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATg                                                              <  5:277331/149‑1 (MQ=255)
                                                                                         cacaGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATgg                                                             <  1:6146/149‑1 (MQ=255)
                                                                                            aGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGcc                                                                          >  8:108373/1‑133 (MQ=255)
                                                                                             gCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGt                                                         <  5:189432/149‑1 (MQ=255)
                                                                                              cccTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTg                                                        >  2:250282/1‑149 (MQ=255)
                                                                                                cTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTgg                                                       >  5:253585/1‑148 (MQ=255)
                                                                                                 tACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCtggt                                                                                         <  6:71533/113‑1 (MQ=255)
                                                                                                     gTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCg                                                                             <  2:134399/121‑1 (MQ=255)
                                                                                                     gTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCg                                                                             >  1:134399/1‑121 (MQ=255)
                                                                                                     gTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCa                                                 >  1:313968/1‑149 (MQ=255)
                                                                                                      ttCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAAGTGCAGATGGCGGTGGAAACCAt                                                >  3:238371/1‑149 (MQ=255)
                                                                                                        cATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCATCATCTACATCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATc                                               <  2:214549/148‑1 (MQ=255)
                                                                                                        cATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATcc                                              >  8:164635/1‑149 (MQ=255)
                                                                                                        cATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATc                                               <  3:249637/148‑1 (MQ=255)
                                                                                                             tgctgcTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCAc                                         >  2:232249/1‑149 (MQ=255)
                                                                                                                  ctACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACAt                                    >  8:47342/1‑149 (MQ=255)
                                                                                                                   tACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGCGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATc                                   <  8:337434/149‑1 (MQ=255)
                                                                                                                     cTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGa                                 <  2:123108/149‑1 (MQ=255)
                                                                                                                      tCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGac                                >  4:287510/1‑149 (MQ=255)
                                                                                                                      tCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCAAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGac                                >  8:46596/1‑149 (MQ=255)
                                                                                                                        aCCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGacca                              >  7:278168/1‑149 (MQ=255)
                                                                                                                        aCCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGacca                              >  4:289189/1‑149 (MQ=255)
                                                                                                                           gC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGAccacca                           >  8:307939/1‑149 (MQ=255)
                                                                                                                                tGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCacta                        >  1:46391/1‑149 (MQ=255)
                                                                                                                                 gCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTTCAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCt                                                                                 <  3:117054/91‑1 (MQ=255)
                                                                                                                                 gCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCt                                                                                 >  4:117054/1‑91 (MQ=255)
                                                                                                                                   aaCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGa                                 <  6:229797/137‑1 (MQ=255)
                                                                                                                                   aaCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGa                                 >  5:229797/1‑137 (MQ=255)
                                                                                                                                   aaCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTc                     >  4:66068/1‑149 (MQ=255)
                                                                                                                                    aCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTc                     >  1:168924/1‑148 (MQ=255)
                                                                                                                                    aCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCg                    <  2:293946/149‑1 (MQ=255)
                                                                                                                                      cAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCa                                          <  8:103690/125‑1 (MQ=255)
                                                                                                                                      cAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCa                                          >  7:103690/1‑125 (MQ=255)
                                                                                                                                       gtcggaccACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGAc                  >  4:58953/7‑148 (MQ=255)
                                                                                                                                            accACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGac                                <  5:85686/128‑1 (MQ=255)
                                                                                                                                            accACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGac                                >  6:85686/2‑129 (MQ=255)
                                                                                                                                             ccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACcgc           >  8:274882/1‑149 (MQ=255)
                                                                                                                                               aCCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCgcgg         >  1:214923/1‑149 (MQ=255)
                                                                                                                                                cccAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCgcggc        <  3:58953/149‑1 (MQ=255)
                                                                                                                                                 ccTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCgcggc        >  7:146042/1‑148 (MQ=255)
                                                                                                                                                     gAGCGCCTGGGTGGTCCGACCATGGGCAGTAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGAtt    <  8:106376/148‑1 (MQ=255)
                                                                                                                                                     gAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGAtt    <  2:128796/148‑1 (MQ=255)
                                                                                                                                                     gAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTc   <  7:47342/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGC‑‑CTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTCC  >  NC_002947/2456846‑2457142

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: