Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,978,102 T→A 5.5% D9V (GAC→GTC)  PP_5240 ← two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,978,1020TA 5.5% 272.3 / 4.6 91D9V (GAC→GTC) PP_5240two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base T (37/48);  new base A (3/2);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

CGGCATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGCTG  >  NC_002947/5977957‑5978250
                                                                                                                                                  |                                                                                                                                                    
cGGCATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTcga                                                                                                                                                   <  3:159994/149‑1 (MQ=255)
cGGCATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTcga                                                                                                                                                   >  6:115816/1‑149 (MQ=255)
   cATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTcgacga                                                                                                                                                <  8:162176/149‑1 (MQ=255)
   cATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTcgacga                                                                                                                                                <  4:81276/149‑1 (MQ=255)
          aGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCCGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCacac                                                                                                                                         <  3:114136/149‑1 (MQ=255)
           ggTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCCGAGGCATCGTCGTCGACGATGCacac                                                                                                                                         <  1:4600/149‑1 (MQ=255)
            gTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACAccc                                                                                                                                       >  2:41342/1‑149 (MQ=255)
             tCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCct                                                                                                                                     <  6:5623/149‑1 (MQ=255)
              ccAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCct                                                                                                                                     >  1:34099/1‑149 (MQ=255)
               cAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCctc                                                                                                                                    >  2:167786/1‑149 (MQ=255)
                  cAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGc                                                                                                                                 <  2:34099/149‑1 (MQ=255)
                  cAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGc                                                                                                                                 >  2:207892/1‑149 (MQ=255)
                    acacacAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCt                                                                                                                                >  3:47546/1‑148 (MQ=255)
                        acaGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATgc                                                                                                                           <  8:339885/149‑1 (MQ=255)
                        acaGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATgc                                                                                                                           >  6:69283/1‑149 (MQ=255)
                          aGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTGGAAGATGCACACCCCTCTGGACATgcgc                                                                                                                         >  4:80740/1‑149 (MQ=255)
                            cccGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGcgccg                                                                                                                       >  5:112596/1‑149 (MQ=255)
                             ccGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGcgccgc                                                                                                                      >  6:18797/1‑149 (MQ=255)
                             ccGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGcgccgc                                                                                                                      >  5:79450/1‑149 (MQ=255)
                             ccGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGcgccgc                                                                                                                      <  8:318503/149‑1 (MQ=255)
                               gCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCgccgcc                                                                                                                     >  2:136802/1‑148 (MQ=255)
                               gCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCgccgcc                                                                                                                     >  5:276893/1‑148 (MQ=255)
                                 cTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAg                                                                                                                  >  6:26138/1‑149 (MQ=255)
                                 cTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAg                                                                                                                  <  5:18797/149‑1 (MQ=255)
                                 cTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCCCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAg                                                                                                                  >  6:129108/1‑149 (MQ=255)
                                    gcccgccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAgg                                                                                                                 <  6:120761/147‑1 (MQ=255)
                                     cccgccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTc                                                                                                              <  6:29637/149‑1 (MQ=255)
                                      ccgccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCg                                                                                                             <  3:370209/149‑1 (MQ=255)
                                      ccgccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCg                                                                                                             <  6:79450/149‑1 (MQ=255)
                                      ccgccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCg                                                                                                             <  1:167786/149‑1 (MQ=255)
                                       cgccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGc                                                                                                            <  1:257483/149‑1 (MQ=255)
                                        gccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCa                                                                                                           <  6:276893/149‑1 (MQ=255)
                                          cAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGa                                                                                                         >  6:151653/1‑149 (MQ=255)
                                          cAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGa                                                                                                         <  6:161534/149‑1 (MQ=255)
                                             tGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACcg                                                                                                      >  4:118624/1‑149 (MQ=255)
                                                     ccAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCCCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCga                                                                                              <  5:129108/149‑1 (MQ=255)
                                                        ggAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCgagga                                                                                           <  6:366243/149‑1 (MQ=255)
                                                         gAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCgaggag                                                                                          <  2:5589/149‑1 (MQ=255)
                                                         gAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCgaggag                                                                                          >  7:111498/1‑149 (MQ=255)
                                                           aTGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGACGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGCCCGCGCTCCGAGGAGGt                                                                                        >  6:15910/1‑149 (MQ=255)
                                                            tGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGtt                                                                                       <  3:118624/149‑1 (MQ=255)
                                                              cGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTtgt                                                                                     <  6:112596/149‑1 (MQ=255)
                                                              cGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTtgt                                                                                     <  2:293697/149‑1 (MQ=255)
                                                               gTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTtgtg                                                                                    >  4:94399/1‑149 (MQ=255)
                                                               gTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTtgtg                                                                                    >  8:282983/1‑149 (MQ=255)
                                                                    cggcggAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGc                                                                                                            >  8:44324/1‑120 (MQ=255)
                                                                     ggcggAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACCCCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGc                                                                              >  1:147581/1‑149 (MQ=255)
                                                                           aaaCGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCCGGGTTTTGCGCACCGAGGCATAGTCGTCGACGATGCCCCCCCCGGGGGGCAAGCGCCGCCCAGGCGGGCCGCCCCCGCGCCGGGGGGGGTGTGCATGGCCGccac                                                                        >  7:241154/1‑146 (MQ=255)
                                                                              cGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTg                                                                                                                                  <  4:221826/88‑1 (MQ=255)
                                                                              cGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTg                                                                                                                                  >  3:221826/1‑88 (MQ=255)
                                                                              cGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATgg                                                                               >  6:242784/1‑139 (MQ=255)
                                                                              cGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATgg                                                                               <  5:242784/139‑1 (MQ=255)
                                                                              cGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTc                                                                     <  1:201518/149‑1 (MQ=255)
                                                                               gACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGAt                                                                                                                                               <  7:300709/74‑1 (MQ=255)
                                                                               gACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGAt                                                                                                                                               >  8:300709/1‑74 (MQ=255)
                                                                               gACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGTGGTTGTGCATGGCCGCAAGCTcc                                                                    >  8:309207/1‑149 (MQ=255)
                                                                               gACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTcc                                                                    <  2:126811/149‑1 (MQ=255)
                                                                                     gTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGccc                                                                                                                    >  6:320788/1‑95 (MQ=255)
                                                                                     gTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGccc                                                                                                                    <  5:320788/95‑1 (MQ=255)
                                                                                     gTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCCTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCagcagc                                                              >  3:357853/1‑149 (MQ=255)
                                                                                       ctcGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGtt                                                                                       <  5:260864/122‑1 (MQ=255)
                                                                                       ctcGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGtt                                                                                       >  6:260864/1‑122 (MQ=255)
                                                                                         cGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTc                                                                     <  3:249033/138‑1 (MQ=255)
                                                                                         cGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTc                                                                     >  4:249033/1‑138 (MQ=255)
                                                                                          gAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCagcgagcg                                                         >  5:1780/1‑149 (MQ=255)
                                                                                          gAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCagcgag                                                           <  8:111498/147‑1 (MQ=255)
                                                                                           aaCCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCgg                                                        <  1:41342/149‑1 (MQ=255)
                                                                                           aaCCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCgg                                                        <  7:282983/149‑1 (MQ=255)
                                                                                           aaCCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGGCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCgg                                                        <  5:115816/149‑1 (MQ=255)
                                                                                                 gccgAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGc                                                  <  3:94399/149‑1 (MQ=255)
                                                                                                     aGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCt                                              <  1:354388/149‑1 (MQ=255)
                                                                                                          agcagGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTc                                         <  5:69283/149‑1 (MQ=255)
                                                                                                            cagGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCat                                       <  6:1780/149‑1 (MQ=255)
                                                                                                               gTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCatgat                                    >  1:226986/1‑149 (MQ=255)
                                                                                                               gTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCatgat                                    <  5:151653/149‑1 (MQ=255)
                                                                                                                   gCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCaa                                                                         <  5:178843/108‑1 (MQ=255)
                                                                                                                   gCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCaa                                                                         >  6:178843/1‑108 (MQ=255)
                                                                                                                   gCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttc                                <  6:178090/149‑1 (MQ=255)
                                                                                                                   gCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttc                                <  2:226986/149‑1 (MQ=255)
                                                                                                                    ccAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttct                               >  4:178585/1‑149 (MQ=255)
                                                                                                                     cAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttctt                              >  3:122378/1‑149 (MQ=255)
                                                                                                                      aGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttcttc                             <  5:26138/149‑1 (MQ=255)
                                                                                                                      aGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttcttc                             <  2:147581/149‑1 (MQ=255)
                                                                                                                      aGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttcttc                             <  4:122378/149‑1 (MQ=255)
                                                                                                                                    ccGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTgg                <  5:121850/149‑1 (MQ=255)
                                                                                                                                        tgCATCGTCGACGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACGGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGgtcg             >  6:222648/2‑148 (MQ=255)
                                                                                                                                         gCATCGTCGTCGACGATGCACACCCGTCTGCTCTTGCGCGGCCCAGGCTCGCAGACCGCGCGCCGAGGAGGTTGTGCATGGCCGCAAGCTCAAGCAGCGTGGGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCTCCTTGgtcgtc           <  5:179/149‑1 (MQ=255)
                                                                                                                                          cATCGTCGACGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGcc                                                                             <  7:67487/82‑1 (MQ=255)
                                                                                                                                          cATCGTCGACGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGcc                                                                             >  8:67487/1‑82 (MQ=255)
                                                                                                                                          cATCGTCGACGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCca                                                                   <  2:176393/92‑1 (MQ=255)
                                                                                                                                          cATCGTCGACGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCca                                                                   >  1:176393/1‑92 (MQ=255)
                                                                                                                                            tcgtcgtcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCAcc        <  1:136802/149‑1 (MQ=255)
                                                                                                                                                cgtcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGc    <  4:357853/149‑1 (MQ=255)
                                                                                                                                                  tcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTcttctt                           >  7:52002/1‑124 (MQ=255)
                                                                                                                                                  tcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTcttctt                           <  8:52002/124‑1 (MQ=255)
                                                                                                                                                  tcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGCTg  >  1:345980/1‑149 (MQ=255)
                                                                                                                                                  tcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGCTg  >  6:192439/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGGCATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGCTG  >  NC_002947/5977957‑5978250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: