Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,951,589 T→C 7.6% Q456Q (CAA→CAG)  ppk ← polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,951,5890TC 7.6% 259.8 / 7.3 91Q456Q (CAA→CAG) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (39/45);  new base C (5/2);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGC  >  NC_002947/5951443‑5951737
                                                                                                                                                  |                                                                                                                                                     
tACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCtt                                                                                                                                                      <  1:252372/148‑1 (MQ=255)
  cAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGc                                                                                                                                                    <  1:241797/148‑1 (MQ=255)
   aGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAg                                                                                                                                                  <  4:352566/149‑1 (MQ=255)
     gcgggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggc                                                                                                                                                >  4:39781/1‑149 (MQ=255)
       gggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                              <  6:312064/149‑1 (MQ=255)
              gCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGcc                                                                                                                                        <  8:121922/148‑1 (MQ=255)
              gCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTggtgc                                                                                                                                       >  4:347003/1‑146 (MQ=255)
               ccGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAt                                                                                                                                      <  8:248382/149‑1 (MQ=255)
               ccGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAt                                                                                                                                      >  8:235290/1‑149 (MQ=255)
               ccGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTggtgcg                                                                                                                                      <  3:347003/149‑5 (MQ=255)
                 ggCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATca                                                                                                                                    <  1:373194/149‑2 (MQ=255)
                 ggCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCt                                                                                                                                    >  6:244913/1‑149 (MQ=255)
                  gCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                              <  7:116195/138‑1 (MQ=255)
                  gCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                              >  8:116195/1‑138 (MQ=255)
                   cGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCATAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGc                                                                                                                                  <  2:151642/149‑1 (MQ=255)
                    gTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTg                                                                                                                                   >  4:115660/1‑147 (MQ=255)
                       gTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCACCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGt                                                                                                                              <  3:371012/149‑1 (MQ=255)
                        tAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGtt                                                                                                                             <  2:101212/149‑1 (MQ=255)
                           ttGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGa                                                                                                                          <  3:149884/149‑1 (MQ=255)
                                 gTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTc                                                                                                                    <  6:117153/149‑1 (MQ=255)
                                 gTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTc                                                                                                                    <  7:118797/149‑1 (MQ=255)
                                    cccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCATCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGt                                                                                                                  <  1:269272/148‑1 (MQ=255)
                                      cAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCggaa                                                                                                                >  1:141846/1‑145 (MQ=255)
                                        gATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCaaaa                                                                                                             <  8:123511/149‑1 (MQ=255)
                                         aTGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATAATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAAc                                                                                                            >  5:178635/1‑149 (MQ=255)
                                            cgcgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAAcgc                                                                                                          <  3:212291/148‑1 (MQ=255)
                                                 aaCGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGcgcacg                                                                                                     >  2:366736/1‑148 (MQ=255)
                                                      aCGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggc                                                                                                                                                <  5:229918/100‑1 (MQ=255)
                                                      aCGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggc                                                                                                                                                >  6:229918/1‑100 (MQ=255)
                                                      aCGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTc                                                                                               <  4:323303/149‑1 (MQ=255)
                                                              gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGAt                                                                                            <  7:77678/144‑1 (MQ=255)
                                                              gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                       <  1:366736/149‑1 (MQ=255)
                                                              gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                       >  8:231006/1‑149 (MQ=255)
                                                              gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                       >  2:215195/1‑149 (MQ=255)
                                                              gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                       <  3:39781/149‑1 (MQ=255)
                                                               cccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                       >  8:66277/1‑148 (MQ=255)
                                                                ccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCa                                                                                     <  5:244913/149‑1 (MQ=255)
                                                                   gCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGt                                                                                   >  5:326520/1‑148 (MQ=255)
                                                                    cTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGCAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGt                                                                                   <  1:204305/147‑1 (MQ=255)
                                                                    cTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGt                                                                                   <  6:112931/147‑1 (MQ=255)
                                                                    cTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGt                                                                                   <  1:132296/147‑1 (MQ=255)
                                                                     tCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAc                                                                                >  1:318346/1‑149 (MQ=255)
                                                                       gcggcgCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcc                                                                               <  3:212717/148‑1 (MQ=255)
                                                                          gcgcAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcctcct                                                                           >  3:232128/1‑149 (MQ=255)
                                                                             cAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGc                                                                        >  8:336414/1‑149 (MQ=255)
                                                                              aGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCCCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCCCCGCAGTGACCTCCTTGcc                                                                       >  4:29679/1‑149 (MQ=255)
                                                                                 aTCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGtt                                                                    <  4:14128/149‑1 (MQ=255)
                                                                                 aTCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGtt                                                                    <  5:123771/149‑1 (MQ=255)
                                                                                 aTCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGt                                                                     >  2:338157/1‑148 (MQ=255)
                                                                                  tCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTa                                                                   >  2:326305/1‑149 (MQ=255)
                                                                                       atcatcTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCATTGACCTCCTTGCCGTTACgcgc                                                              <  7:235290/149‑1 (MQ=255)
                                                                                                gCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑ccac                                                    >  8:127722/1‑149 (MQ=255)
                                                                                                gCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCGTGCCGTTACGCGCCTCAT‑ccac                                                    >  7:343772/1‑149 (MQ=255)
                                                                                                  gTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCA‑CAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAg                                                 >  1:251640/1‑149 (MQ=255)
                                                                                                   tGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAg                                                 >  5:240390/1‑149 (MQ=255)
                                                                                                   tGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAg                                                 >  8:264768/1‑149 (MQ=255)
                                                                                                           gAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcctc                                                                             <  1:260244/114‑1 (MQ=255)
                                                                                                           gAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcctc                                                                             >  2:260244/1‑114 (MQ=255)
                                                                                                           gAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTc                                           <  7:66277/147‑1 (MQ=255)
                                                                                                               ccAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGAt                                     <  3:111289/149‑1 (MQ=255)
                                                                                                                cAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGt                                            >  2:36959/1‑141 (MQ=255)
                                                                                                                cAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGt                                            <  1:36959/141‑1 (MQ=255)
                                                                                                                   ccACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTcc                                 <  6:326520/149‑1 (MQ=255)
                                                                                                                     acacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGa                               <  1:326305/149‑1 (MQ=255)
                                                                                                                      cacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAg                              <  2:318346/149‑1 (MQ=255)
                                                                                                                        ccGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCCACCAGGGCGTCGACGATTTCCGAGt                             <  1:215195/149‑1 (MQ=255)
                                                                                                                            aGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGcgcac                                                                                                      <  7:247732/72‑1 (MQ=255)
                                                                                                                            aGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGcgcac                                                                                                      >  8:247732/1‑72 (MQ=255)
                                                                                                                               tCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGccccc                     >  6:367669/1‑149 (MQ=255)
                                                                                                                               tCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGccccc                     >  5:61138/1‑149 (MQ=255)
                                                                                                                                 aCTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGa                   <  7:127722/149‑1 (MQ=255)
                                                                                                                                 aCTACCGCACCAGCCGCTTGCAGGCGGCTGGCAATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGccccc                     >  6:129228/1‑147 (MQ=255)
                                                                                                                                   tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGtt                            >  8:18181/1‑138 (MQ=255)
                                                                                                                                   tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGtt                            <  7:18181/138‑1 (MQ=255)
                                                                                                                                   tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGc                 >  7:60767/1‑149 (MQ=255)
                                                                                                                                   tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGAGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGc                 >  1:66197/1‑149 (MQ=255)
                                                                                                                                   tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTACGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGACGACGAATTCCGAGTTAGCCACCGAGc                 >  3:72710/1‑149 (MQ=255)
                                                                                                                                    ggatggccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCg                <  6:251626/143‑1 (MQ=255)
                                                                                                                                     ccGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCAGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCagcggagcgg               >  8:140428/1‑149 (MQ=255)
                                                                                                                                       tggccAGCCGCCTGCAAGCGGCTGGTCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTa             >  7:101465/4‑149 (MQ=255)
                                                                                                                                        cACCAGCCGCTTGCAGGCGGCTGGCCAGCTGCAGGTTGGAGTCTGCGTCCAAACGCGCCCGCAACTCGATGACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAG‑GCGCCAGGGCGTCGCCGAGTTCCGAGGTAGCCGCCGAGCGGTag            >  4:222394/1‑148 (MQ=255)
                                                                                                                                         aCCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTa             >  4:245433/1‑147 (MQ=255)
                                                                                                                                         aCCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACa           <  3:115660/149‑1 (MQ=255)
                                                                                                                                          ccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAg          >  7:309005/1‑149 (MQ=255)
                                                                                                                                          ccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAg          >  3:168474/1‑149 (MQ=255)
                                                                                                                                           cAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg         <  5:367669/149‑1 (MQ=255)
                                                                                                                                           cAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg         >  1:225462/1‑149 (MQ=255)
                                                                                                                                            aGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAggg        >  3:127221/1‑149 (MQ=255)
                                                                                                                                             gCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg         >  7:29176/1‑147 (MQ=255)
                                                                                                                                             gCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg         <  8:29176/147‑1 (MQ=255)
                                                                                                                                                gCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTg    <  1:338157/149‑1 (MQ=255)
                                                                                                                                                 cTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTAC‑GGGTCTGgc  <  2:251640/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGC  >  NC_002947/5951443‑5951737

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