Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 605,374 G→C 5.8% V229V (GTG→GTC thiL → thiamin monophosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947605,3740GC 5.8% 229.0 / 3.8 72V229V (GTG→GTCthiLthiamin monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (31/35);  new base C (2/2);  total (34/38)
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.58e-01

CGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCTGGCGCCCCTGG  >  NC_002947/605226‑605522
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                      <  3:204728/149‑1 (MQ=255)
                                tCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAgt                                                                                                                      >  2:232614/1‑149 (MQ=255)
                                tCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAgt                                                                                                                      >  8:88194/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                                                                                                                                <  7:147280/121‑1 (MQ=255)
                                                  aCCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGgcgc                                                                                                                                >  8:147280/1‑121 (MQ=255)
                                                  aCCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAgg                                                                                                     >  8:85285/1‑148 (MQ=255)
                                                  aCCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAgg                                                                                                     >  7:80192/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                                                                                                                                             <  4:126765/81‑1 (MQ=255)
                                                                             gCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAg                                                                                                                                             >  3:126765/1‑81 (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                                                                                                                                                      >  5:265376/1‑49 (MQ=255)
                                                                                                    cGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTg                                                                                                                                                      <  6:265376/49‑1 (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)
                                                                                                                                                gCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCTGGCGccc      >  3:221024/1‑149 (MQ=255)
                                                                                                                                                gCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCTGGCGccc      >  8:87604/1‑149 (MQ=255)
                                                                                                                                                 cGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCTGGCGcccc     <  8:108542/149‑1 (MQ=255)
                                                                                                                                                    ggCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCTGGCGCCCctgg  <  3:176079/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCGATGGCCTGCTCGCTGACTGCGGGCATATCGCCAAGGCCTCCGGCGTGGCGCTGGAGGTGAACCTGGCGCAGGTGCCAGTGTCTCCTGCCTTAGAGGCATTCCTCGGGCGCGAAGCGGCGTTGCAGGCAGCACTCACCGGTGGCGACGACTATGTACTGGCCTTCACCTTGCCGCCCGAAGCGCTGGCGCCCCTGG  >  NC_002947/605226‑605522

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: