Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 605,369 G→A 7.0% G228S (GGC→AGC)  thiL → thiamin monophosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947605,3690GA 7.0% 209.6 / 3.9 72G228S (GGC→AGC) thiLthiamin monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base A (3/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

CGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCTGGC  >  NC_002947/605226‑605514
                                                                                                                                               |                                                                                                                                                 
cgcAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTg                                                                                                                                              <  6:249543/149‑1 (MQ=255)
  cAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGgc                                                                                                                                            <  6:312696/149‑1 (MQ=255)
  cAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGgc                                                                                                                                            <  5:287186/149‑1 (MQ=255)
     gCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCt                                                                                                                                         <  8:98113/149‑1 (MQ=255)
     gCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCt                                                                                                                                         <  2:271561/149‑1 (MQ=255)
        gAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTg                                                                                                                                        >  7:243587/1‑147 (MQ=255)
         aGCCGCTGCTGGCGCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAg                                                                                                                                     <  4:188923/149‑1 (MQ=255)
           ccGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGt                                                                                                                                   >  5:226129/1‑149 (MQ=255)
             gctgctGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTTGCGCTGGAGGTg                                                                                                                                  >  8:55635/1‑148 (MQ=255)
              ctgctgGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGaa                                                                                                                                >  1:250615/1‑149 (MQ=255)
                    gCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCAGCGCCACGCCGGAggcct                                                                                                                                 >  4:148570/1‑139 (MQ=255)
                    gCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCAGCGCCACGCCGGAggcct                                                                                                                                 <  3:148570/142‑4 (MQ=255)
                        aTTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAg                                                                                                                      <  7:293402/149‑1 (MQ=255)
                         ttACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAgg                                                                                                                     <  1:98724/149‑1 (MQ=255)
                               gTCACCGCTGCCGCAGCGCACCCTCGGGGCTCTGCTGCGGGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAg                                                                                                               <  4:101810/149‑1 (MQ=255)
                               gTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGt                                                                                                                    <  8:289555/144‑1 (MQ=255)
                               gTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGt                                                                                                                    >  7:289555/1‑144 (MQ=255)
                               gTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCa                                                                                                                <  8:114845/148‑1 (MQ=255)
                               gTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCAGCGCCACGCCGGAggcct                                                                                                                                 >  3:282007/1‑128 (MQ=255)
                               gTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCAGCGCCACGCCGGAggcct                                                                                                                                 <  4:282007/131‑4 (MQ=255)
                                tCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAgt                                                                                                              >  8:88194/1‑149 (MQ=255)
                                tCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAgt                                                                                                              <  3:204728/149‑1 (MQ=255)
                                tCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAgt                                                                                                              >  2:232614/1‑149 (MQ=255)
                                     gATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGgctc                                                                                                         <  6:108938/149‑1 (MQ=255)
                                         ccGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGc                                                                                                     <  1:161737/149‑1 (MQ=255)
                                             aGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTa                                                                                                 <  3:119548/149‑1 (MQ=255)
                                                  aCCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGgcgc                                                                                                                        >  8:147280/1‑121 (MQ=255)
                                                  aCCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGgcgc                                                                                                                        <  7:147280/121‑1 (MQ=255)
                                                  aCCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAgg                                                                                             >  7:80192/1‑148 (MQ=255)
                                                  aCCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAgg                                                                                             >  8:85285/1‑148 (MQ=255)
                                                  aCCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGc                                                                                            <  5:110150/149‑1 (MQ=255)
                                                    ccTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCAt                                                                                          >  7:173783/1‑149 (MQ=255)
                                                      tCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTc                                                                                        >  7:78205/1‑149 (MQ=255)
                                                            ctctGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCa                                                                                                                >  5:175099/1‑119 (MQ=255)
                                                            ctctGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCa                                                                                                                <  6:175099/119‑1 (MQ=255)
                                                                  tgcGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGgc                                                                                                                          <  1:146314/103‑1 (MQ=255)
                                                                  tgcGTGGCAAGGCCACTGCGGCCCTGGACAATTCCGATGGCATGCTCTCTGACTGCGCGAATATCGCCAAGGCCTCCGGCGTGTCGCAGGAGGTGAACCTGgc                                                                                                                          >  2:146314/1‑103 (MQ=255)
                                                                     gTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAgcg                                                                         <  8:80192/149‑1 (MQ=255)
                                                                        gCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAgcggcg                                                                      >  1:46775/1‑149 (MQ=255)
                                                                          aaGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGtt                                                                    <  4:297267/149‑1 (MQ=255)
                                                                             gCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAg                                                                                                                                     >  3:126765/1‑81 (MQ=255)
                                                                             gCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAg                                                                                                                                     <  4:126765/81‑1 (MQ=255)
                                                                             gCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTgca                                                                 <  1:232614/149‑1 (MQ=255)
                                                                             gCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGAATATCTCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACATGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTAGGGCGAGAAGCGTCGTTg                                                                   >  4:199930/1‑147 (MQ=255)
                                                                                    cGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTgcaggcag                                                            <  8:78205/147‑1 (MQ=255)
                                                                                        acTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCa                                                      >  8:224336/2‑149 (MQ=255)
                                                                                         cTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCAc                                                     >  1:237676/1‑149 (MQ=255)
                                                                                              cATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTg                                                <  7:224336/149‑1 (MQ=255)
                                                                                                 ttCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGc                                              >  8:140906/1‑148 (MQ=255)
                                                                                                   ccGATTGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCgtattcgc                                            <  6:226129/148‑8 (MQ=255)
                                                                                                   ccGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTTGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGcga                                            <  7:55635/148‑1 (MQ=255)
                                                                                                    cGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTg                                                                                                                                              <  6:265376/49‑1 (MQ=255)
                                                                                                    cGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTg                                                                                                                                              >  5:265376/1‑49 (MQ=255)
                                                                                                    cGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGcgac                                           >  5:126290/1‑148 (MQ=255)
                                                                                                     gATGGCCTTCTCGCAGACTTAGGGAATATCGGCAAGGACTCCTTCGTGGCGCTGGAGGTGAACCTGGCGCAGTTGCCAGTGTCTCCTGCCTTAGATGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGcgacga                                         <  4:212699/149‑1 (MQ=255)
                                                                                                     gATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGcgacga                                         <  7:300532/149‑1 (MQ=255)
                                                                                                     gATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGcgacga                                         <  5:269720/149‑1 (MQ=255)
                                                                                                           cTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATg                                    >  3:71379/1‑148 (MQ=255)
                                                                                                                 gCTGACTGCGGGCGTATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTgcaggcag                                                            <  7:173398/118‑1 (MQ=255)
                                                                                                                 gCTGACTGCGGGCGTATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTgcaggcag                                                            >  8:173398/1‑118 (MQ=255)
                                                                                                                 gCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTg                               <  7:85285/147‑1 (MQ=255)
                                                                                                                 gCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGc                             >  5:17228/1‑149 (MQ=255)
                                                                                                                          gggCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCtt                    <  6:17228/149‑1 (MQ=255)
                                                                                                                           ggCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTg                   <  2:46775/149‑1 (MQ=255)
                                                                                                                               tatCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTgccgc               >  7:108542/1‑149 (MQ=255)
                                                                                                                                atCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTgccgcc              >  5:5278/1‑149 (MQ=255)
                                                                                                                                 tCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGccc             >  4:176079/1‑149 (MQ=255)
                                                                                                                                  cGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGccc             >  3:276717/1‑148 (MQ=255)
                                                                                                                                        ggCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAgcgc      <  2:250615/149‑1 (MQ=255)
                                                                                                                                         gCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCt     <  7:140906/149‑1 (MQ=255)
                                                                                                                                            tCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCTGgc  >  6:61842/1‑149 (MQ=255)
                                                                                                                                            tCCAGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCTGgc  >  1:134613/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCTGGC  >  NC_002947/605226‑605514

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: