Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 605,131 T→C 6.6% G148G (GGT→GGC thiL → thiamin monophosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947605,1310TC 6.6% 261.1 / 6.5 90G148G (GGT→GGCthiLthiamin monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (43/41);  new base C (4/2);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

GCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACC  >  NC_002947/604984‑605278
                                                                                                                                                   |                                                                                                                                                   
gctggcGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCTTCACCGTGAGCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGt                                                                                                                                                     <  5:169689/148‑1 (MQ=255)
gctggcGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTg                                                                                                                                                    >  3:298620/1‑149 (MQ=255)
      ggcCTATGCGGACGGCCTGGGTGGCAGGGCCCACGGTTGCCGCATTAGCGTTTTCGGTGGTTACACAAGGCGCGGCCCGTTTAGCATCACCGTGACCGTGTTCGGCCGGGTGCCGGCTGGCCAGGCCTTGCGCCGAAGCGGTGCAAGGc                                                                                                                                              <  8:64585/149‑1 (MQ=255)
      ggcCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGc                                                                                                                                              <  6:261898/149‑1 (MQ=255)
         cTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCACCGCTGCGGCGCAAg                                                                                                                                                <  7:1275/144‑1 (MQ=255)
         cTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCACCGCTGCGGCGCAAg                                                                                                                                                >  8:5821/1‑144 (MQ=255)
         cTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCACCGCTGCGGCGCAAg                                                                                                                                                <  7:5821/144‑1 (MQ=255)
         cTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCACCGCTGCGGCGCAAg                                                                                                                                                >  8:1275/1‑144 (MQ=255)
          tATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCggg                                                                                                                                          >  5:287641/1‑149 (MQ=255)
            tGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGt                                                                                                                                         >  5:232449/1‑148 (MQ=255)
               cGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGAt                                                                                                                                      >  2:185910/1‑148 (MQ=255)
               cGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGAt                                                                                                                                      <  1:185910/148‑1 (MQ=255)
                   ggCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgtt                                                                                                                                 <  7:292190/149‑1 (MQ=255)
                   ggCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgtt                                                                                                                                 <  8:198394/149‑1 (MQ=255)
                     ccTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgttgt                                                                                                                               <  6:232449/149‑1 (MQ=255)
                     ccTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgttgt                                                                                                                               <  1:183784/149‑1 (MQ=255)
                     cattGGGCGCTGGGCCCACCGTGGCCGCCGGCGCGTTTTCGGTGGTGACACCCGGCGCGGCCGTTTGCGCATCACCGTGACCGTTTGCGGCCGCGTGCGGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCATGGCCGGGTTATGtgttgt                                                                                                                               <  4:298620/147‑1 (MQ=255)
                      cTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtg                                                                                                                              >  6:288644/1‑149 (MQ=255)
                      cTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtg                                                                                                                              >  1:37761/1‑149 (MQ=255)
                      cTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtg                                                                                                                              >  1:286335/1‑149 (MQ=255)
                         ggTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTgtgtgt                                                                                                                           >  1:66210/1‑149 (MQ=255)
                             gCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCggt                                                                                                                       <  6:287641/149‑1 (MQ=255)
                              cATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCggtg                                                                                                                      <  4:115085/149‑1 (MQ=255)
                               aTGGGCCACCGGTGGCGCATGTGCCTGATCGGTGGTGACACCACGCGCGGCCCGGTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCGAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCggtg                                                                                                                      <  8:131448/148‑1 (MQ=255)
                               aTGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCACCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCggtgg                                                                                                                     >  1:149614/1‑149 (MQ=255)
                                    ccACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGc                                                                                                                >  5:221854/1‑149 (MQ=255)
                                       ccGTTGCGGCATGTGCCTGATCGGTGGTGACACCACGCGCGGGCCGTTGAGCATCACCGTGACGGTGTTGGGCCGGGTGCCGGCTGGCAAGGCCTTGCGCCGCTGCGGTGCTATGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTgg                                                                                                             <  2:149614/149‑1 (MQ=255)
                                         gTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTggg                                                                                                            >  2:130434/1‑148 (MQ=255)
                                             ccGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTg                                                                                                              >  2:141414/1‑142 (MQ=255)
                                             ccGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTg                                                                                                              <  1:141414/142‑1 (MQ=255)
                                             ccGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAgg                                                                                                       >  4:131946/1‑149 (MQ=255)
                                             ccGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGGTGAGCATCACCGTGACGGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTTCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAg                                                                                                        <  7:162453/148‑1 (MQ=255)
                                              cGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGgc                                                                                                      >  2:180735/1‑149 (MQ=255)
                                               gCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGgcc                                                                                                     >  5:189449/1‑149 (MQ=255)
                                                 aTGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGgccgc                                                                                                   <  6:92818/149‑1 (MQ=255)
                                                   gAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGccgccg                                                                                                 >  1:148133/1‑149 (MQ=255)
                                                      ccTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgttgt                                                                                                                               >  5:250331/1‑116 (MQ=255)
                                                      ccTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgttgt                                                                                                                               <  6:250331/116‑1 (MQ=255)
                                                      ccTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTgg                                                                                                             <  4:196548/134‑1 (MQ=255)
                                                      ccTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTgg                                                                                                             >  3:196548/1‑134 (MQ=255)
                                                      ccTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTg                                                                                              >  8:123750/1‑149 (MQ=255)
                                                      ccTGATCGGTGGTGACACCACGCGCAGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTg                                                                                              >  8:61211/1‑149 (MQ=255)
                                                       cTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGc                                                                                             <  7:123750/149‑1 (MQ=255)
                                                         gATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCt                                                                                           >  5:96636/1‑149 (MQ=255)
                                                               tggtgACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTg                                                                                                              >  5:161385/1‑124 (MQ=255)
                                                               tggtgACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTg                                                                                                              <  6:161385/124‑1 (MQ=255)
                                                                 gtgACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTg                                                                                   <  6:172793/149‑1 (MQ=255)
                                                                 gtgACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTg                                                                                   <  8:192332/149‑1 (MQ=255)
                                                                      accacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCt                                                                              <  2:286335/149‑1 (MQ=255)
                                                                      accacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCt                                                                              >  1:209825/1‑149 (MQ=255)
                                                                        cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCttg                                                                            <  1:130434/149‑3 (MQ=255)
                                                                        cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTgg                                                                            <  2:66210/149‑1 (MQ=255)
                                                                        cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTgg                                                                            <  2:44507/149‑1 (MQ=255)
                                                                        cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTgg                                                                            <  6:189449/149‑1 (MQ=255)
                                                                        cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTgg                                                                            <  2:209825/149‑1 (MQ=255)
                                                                        cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTgg                                                                            <  1:44626/149‑1 (MQ=255)
                                                                             gcgGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAg                                                                       <  5:288644/149‑1 (MQ=255)
                                                                                  ccGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAgcgc                                                                    >  1:131119/1‑147 (MQ=255)
                                                                                   cGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAg                                                                 <  1:180735/149‑1 (MQ=255)
                                                                                   cGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAg                                                                 <  6:96636/149‑1 (MQ=255)
                                                                                        aGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGgcgc                                                             <  4:92275/148‑1 (MQ=255)
                                                                                          cATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAgg                                                                >  7:124767/1‑143 (MQ=255)
                                                                                          cATCACCGTGACCGTGTTCGGCCGCGGGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAgg                                                                <  8:124767/143‑1 (MQ=255)
                                                                                           aTCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGc                                                         >  2:307338/1‑149 (MQ=255)
                                                                                            tCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGcc                                                        <  2:131119/149‑1 (MQ=255)
                                                                                             cACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCcg                                                       >  4:66713/1‑149 (MQ=255)
                                                                                               ccgtgaccgtgTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCcgcg                                                     <  5:181549/149‑1 (MQ=255)
                                                                                               ccgtgaccgtgTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCcgcg                                                     >  6:181549/1‑149 (MQ=255)
                                                                                                     ccgtgTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGgcgc                                                             <  7:141264/135‑1 (MQ=255)
                                                                                                     ccgtgTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGgcgc                                                             >  8:141264/1‑135 (MQ=255)
                                                                                                      cgtgTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGcc                                               >  3:118981/1‑148 (MQ=255)
                                                                                                               ccGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCgctgc                                     >  2:292806/1‑149 (MQ=255)
                                                                                                                 gcgTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGctgctg                                   >  4:177148/1‑149 (MQ=255)
                                                                                                                   gTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGc                                 >  2:221141/1‑149 (MQ=255)
                                                                                                                      ccGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCa                              >  5:304491/1‑149 (MQ=255)
                                                                                                                        ggCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCAtt                            <  2:37761/149‑1 (MQ=255)
                                                                                                                         gCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCAGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTa                           >  6:185680/1‑149 (MQ=255)
                                                                                                                            cacccggCCTTGCACCGCTGCGGCGCAAGGCCAGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTg                        >  8:62738/6‑149 (MQ=255)
                                                                                                                             gCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTgg                       >  3:230423/1‑149 (MQ=255)
                                                                                                                             gCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTg                        <  7:61211/148‑1 (MQ=255)
                                                                                                                             acccggCCTTGCACCGCTGCGGCGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTgg                       >  5:291477/5‑149 (MQ=255)
                                                                                                                              ccAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGt                      >  4:137961/1‑149 (MQ=255)
                                                                                                                                  gCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCAcc                  >  2:19032/1‑149 (MQ=255)
                                                                                                                                     ttGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGAt               >  4:254403/1‑149 (MQ=255)
                                                                                                                                       gcgcCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATg              <  5:141558/148‑1 (MQ=255)
                                                                                                                                        cgccgcAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGc             <  6:304491/148‑1 (MQ=255)
                                                                                                                                         gccgcAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCg           >  3:223617/1‑149 (MQ=255)
                                                                                                                                         gccgcAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCg           >  3:257532/1‑149 (MQ=255)
                                                                                                                                           cgcAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCa         <  1:19032/149‑1 (MQ=255)
                                                                                                                                            gcAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGcagcag        >  4:155917/1‑149 (MQ=255)
                                                                                                                                                  gTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCAcc  >  4:249149/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACC  >  NC_002947/604984‑605278

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: