Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,628,386 A→C 6.7% A16A (GCA→GCC)  clpP → ATP‑dependent Clp protease proteolytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,628,3860AC 6.7% 217.9 / 6.7 75A16A (GCA→GCC) clpPATP‑dependent Clp protease proteolytic subunit
Reads supporting (aligned to +/- strand):  ref base A (33/37);  new base C (3/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

TGAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAACT  >  NC_002947/2628241‑2628529
                                                                                                                                                 |                                                                                                                                                
tgAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAggc                                                                                                                                               >  4:158611/1‑149 (MQ=255)
  aGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAggcgg                                                                                                                                             <  6:179368/149‑1 (MQ=255)
     aaCGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCt                                                                                                                                          >  3:208915/1‑149 (MQ=255)
      aCGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTg                                                                                                                                         <  4:237232/149‑1 (MQ=255)
        ggTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGt                                                                                                                                       <  1:138808/149‑1 (MQ=255)
        ggTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGt                                                                                                                                       >  4:144136/1‑149 (MQ=255)
         gTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTc                                                                                                                                      <  3:305518/149‑1 (MQ=255)
            ttctCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTccc                                                                                                                                    >  7:214061/1‑148 (MQ=255)
                 tcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGt                                                                                                                              >  4:1663/1‑149 (MQ=255)
                   gTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATgg                                                                                                                               <  8:316726/146‑1 (MQ=255)
                   gTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATgg                                                                                                                               <  6:340276/146‑1 (MQ=255)
                       gAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTggatgtca                                                                                                                                  >  4:316188/1‑133 (MQ=255)
                       gAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTggatgtca                                                                                                                                  <  3:316188/139‑7 (MQ=255)
                          aCCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCa                                                                                                                     <  8:214061/149‑1 (MQ=255)
                            cTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGt                                                                                                                   >  3:58636/1‑149 (MQ=255)
                             tCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTc                                                                                                                  >  4:105405/1‑149 (MQ=255)
                               cacaAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCg                                                                                                                >  7:259842/1‑149 (MQ=255)
                               cacaAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGAAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCg                                                                                                                >  4:301975/1‑149 (MQ=255)
                                acaAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGc                                                                                                               >  7:310728/1‑149 (MQ=255)
                                acaAGCCAGCCTTCGAGCTGGCATATGCGTATTAAAGCCATGACAATTTGGGAATTAACGCAGGACATGTCCCGCAATTAATATTTTAAGCAGAGCTATTAAATACAGGCCGCAGGCGGCCTTGTACCGATTGTTATAGAGCAGTCCGc                                                                                                               >  7:53692/1‑149 (MQ=255)
                                 caAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGcc                                                                                                              <  2:221050/149‑1 (MQ=255)
                                    gccagccTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGt                                                                                                           >  8:100071/1‑149 (MQ=255)
                                     cacgccTTCGAGGTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCTGGACATGTCCCGCAATTGTTATATTGAGGAGAGCTCTCTCATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCGGTg                                                                                                          <  6:260993/146‑1 (MQ=255)
                                         ccTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGa                                                                                                      <  3:144136/149‑1 (MQ=255)
                                          tcgaaggCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGaa                                                                                                     >  4:246898/7‑149 (MQ=255)
                                             aaggCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAAcgc                                                                                                  <  3:246898/146‑1 (MQ=255)
                                              gAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAAcgcg                                                                                                 >  3:208896/1‑149 (MQ=255)
                                                gCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATc                                                                                                                          >  1:271520/1‑122 (MQ=255)
                                                gCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATc                                                                                                                          <  2:271520/122‑1 (MQ=255)
                                                    gCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAg                                                                                                                       <  7:148632/121‑1 (MQ=255)
                                                    gCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAg                                                                                                                       >  8:148632/1‑121 (MQ=255)
                                                    gCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTAc                                                                                            <  3:158611/148‑1 (MQ=255)
                                                    gCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTAc                                                                                            <  8:53692/148‑1 (MQ=255)
                                                    gCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTAc                                                                                            <  4:208896/148‑1 (MQ=255)
                                                              aTTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCg                                                                                                                >  2:183616/1‑118 (MQ=255)
                                                              aTTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCg                                                                                                                <  1:183616/118‑1 (MQ=255)
                                                                tCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTcgc                                                                              >  5:215709/1‑149 (MQ=255)
                                                                  aaGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCgcgc                                                                            <  3:105405/149‑1 (MQ=255)
                                                                  aaGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCgcgc                                                                            >  5:190705/1‑149 (MQ=255)
                                                                  aaGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCgcgc                                                                            >  5:277889/1‑149 (MQ=255)
                                                                    gCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCt                                                                          <  4:208915/149‑1 (MQ=255)
                                                                     ccATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTg                                                                         <  6:277889/149‑1 (MQ=255)
                                                                      cATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGt                                                                        <  3:1663/149‑1 (MQ=255)
                                                                        tGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGtt                                                                       >  6:166679/1‑148 (MQ=255)
                                                                            tATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAgg                                                                  <  3:301975/149‑1 (MQ=255)
                                                                              tttGGGAGTGAGCGCAGGACATGTACCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAg                                                                >  8:256392/1‑149 (MQ=255)
                                                                                    aGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCgtgt                                                           <  7:81241/148‑1 (MQ=255)
                                                                                    aGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCgtgt                                                           >  8:81241/1‑148 (MQ=255)
                                                                                            cAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCt                                                  <  8:259842/149‑1 (MQ=255)
                                                                                            cAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCt                                                  <  7:100071/149‑1 (MQ=255)
                                                                                             aGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTg                                                 >  4:289646/1‑149 (MQ=255)
                                                                                                 cttGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACTTCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCTAGCAGTCCGCCCGTTGCGAACGCTCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGtt                                              <  1:268114/146‑1 (MQ=255)
                                                                                                 cATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTg                                             <  1:222323/149‑1 (MQ=255)
                                                                                                 cATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTg                                             <  6:215709/149‑1 (MQ=255)
                                                                                                    gTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGAAGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCTCGCATGTTCAAGGAGCGTGTGATCTTCCTGGtt                                              >  2:268114/1‑145 (MQ=255)
                                                                                                       ccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCgtg                                                            <  6:248720/128‑1 (MQ=255)
                                                                                                       ccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCgtg                                                            >  5:248720/1‑128 (MQ=255)
                                                                                                           aaTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTaga                                   <  6:355689/149‑1 (MQ=255)
                                                                                                            aTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                  >  4:174902/1‑149 (MQ=255)
                                                                                                            aTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                  <  8:310728/149‑1 (MQ=255)
                                                                                                            aTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                  >  2:204839/1‑149 (MQ=255)
                                                                                                              tCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGa                                >  1:130403/1‑149 (MQ=255)
                                                                                                                  atatTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGAc                               <  4:58636/146‑1 (MQ=255)
                                                                                                                  atatTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTa                             >  3:229084/1‑148 (MQ=255)
                                                                                                                  atatTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTa                             <  5:166679/148‑1 (MQ=255)
                                                                                                                   tatTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTAc                            <  6:190705/148‑1 (MQ=255)
                                                                                                                   tatTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACa                           >  5:58599/1‑149 (MQ=255)
                                                                                                                           agagCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAAc                   <  4:229084/149‑1 (MQ=255)
                                                                                                                               ctctGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCtcg               >  8:283011/1‑149 (MQ=255)
                                                                                                                                  tGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCtcgtcg            >  4:320861/1‑149 (MQ=255)
                                                                                                                                      gaccAGGCCGCCTGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgc        >  3:165339/3‑149 (MQ=255)
                                                                                                                                        ccAGGCCGCCTGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTTACAGCCGCGTGTTGAAGGACCGTGTGATGTTCCTGGTTGGCCCGGTTGAGGAGTACATGGCCACCGTGGTCGTTgcg       >  1:111923/1‑149 (MQ=255)
                                                                                                                                         cAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgcgc      <  7:283011/148‑1 (MQ=255)
                                                                                                                                            gCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAAct  <  3:174902/149‑1 (MQ=255)
                                                                                                                                            gCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAAct  <  1:204839/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATC‑TACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAACT  >  NC_002947/2628241‑2628529

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: ‑