Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 605,126 A→T 7.2% S147C (AGC→TGC)  thiL → thiamin monophosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947605,1260AT 7.2% 265.7 / 7.8 92S147C (AGC→TGC) thiLthiamin monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base A (44/41);  new base T (4/3);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAG  >  NC_002947/604979‑605272
                                                                                                                                                   |                                                                                                                                                  
gACTGGCTGGCGGCCTATGCCGACGGACTTGGTCGCATGGACCACCGTTGCCGCATGAGCCTGATCGGTGGTGAAACCACGCGCGGCCCGTTGAGCATCACCGAGACCGTGTTCGGCCGCGTGCAGGCTGGCCAGGCCTTGCGCCGCAg                                                                                                                                                   >  8:162453/1‑149 (MQ=255)
 aCTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGc                                                                                                                                                  >  7:64585/1‑149 (MQ=255)
     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                                                                                                                                          >  8:5821/1‑144 (MQ=255)
              cTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCACCGCTGCGGCGCAAg                                                                                                                                          <  7:1275/144‑1 (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                                                                                                                                <  1:185910/148‑1 (MQ=255)
                    cGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGAt                                                                                                                                >  2:185910/1‑148 (MQ=255)
                        ggCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgtt                                                                                                                           <  8:198394/149‑1 (MQ=255)
                        ggCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgtt                                                                                                                           <  7:292190/149‑1 (MQ=255)
                          ccTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgttgt                                                                                                                         <  1:183784/149‑1 (MQ=255)
                          ccTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgttgt                                                                                                                         <  6:232449/149‑1 (MQ=255)
                          cattGGGCGCTGGGCCCACCGTGGCCGCCGGCGCGTTTTCGGTGGTGACACCCGGCGCGGCCGTTTGCGCATCACCGTGACCGTTTGCGGCCGCGTGCGGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCATGGCCGGGTTATGtgttgt                                                                                                                         <  4:298620/147‑1 (MQ=255)
                           cTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtg                                                                                                                        >  6:288644/1‑149 (MQ=255)
                           cTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtg                                                                                                                        >  1:286335/1‑149 (MQ=255)
                           cTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtg                                                                                                                        >  1:37761/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                                                                                                        <  1:141414/142‑1 (MQ=255)
                                                  ccGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTg                                                                                                        >  2:141414/1‑142 (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                                                                                                                         <  6:250331/116‑1 (MQ=255)
                                                           ccTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCtgttgt                                                                                                                         >  5:250331/1‑116 (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                                                                                                        <  6:161385/124‑1 (MQ=255)
                                                                    tggtgACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTg                                                                                                        >  5:161385/1‑124 (MQ=255)
                                                                      gtgACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTg                                                                             <  8:192332/149‑1 (MQ=255)
                                                                      gtgACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTg                                                                             <  6:172793/149‑1 (MQ=255)
                                                                           accacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCt                                                                        >  1:209825/1‑149 (MQ=255)
                                                                           accacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCt                                                                        <  2:286335/149‑1 (MQ=255)
                                                                             cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCttg                                                                      <  1:130434/149‑3 (MQ=255)
                                                                             cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTgg                                                                      <  1:44626/149‑1 (MQ=255)
                                                                             cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTgg                                                                      <  2:209825/149‑1 (MQ=255)
                                                                             cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTgg                                                                      <  2:44507/149‑1 (MQ=255)
                                                                             cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTgg                                                                      <  2:66210/149‑1 (MQ=255)
                                                                             cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTgg                                                                      <  6:189449/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                                                       >  8:141264/1‑135 (MQ=255)
                                                                                                          ccgtgTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGgcgc                                                       <  7:141264/135‑1 (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)
                                                                                                                                                   |                                                                                                                                                  
GACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAG  >  NC_002947/604979‑605272

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: