Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,108,470 C→G 5.4% R69R (CGC→CGG)  PP_4494 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,108,4700CG 5.4% 237.1 / 3.5 74R69R (CGC→CGG) PP_4494LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (35/35);  new base G (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

CCGGCATGAAAGCTACACCCGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGT  >  NC_002947/5108323‑5108618
                                                                                                                                                   |                                                                                                                                                    
ccGGCATGAAAGCTACACCCGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGgccgcc                                                                                                                                                     >  2:158861/1‑149 (MQ=255)
        acaGCCACACCCGCTCCGCCGCCGAACTGTCTCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGGTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCa                                                                                                                                              <  4:302855/146‑1 (MQ=255)
        aaaGCTACACCCGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGGCGGCCGACATCGGTCAGCTGTAcctg                                                                                                                                             <  7:217014/149‑3 (MQ=255)
        aaaGCTACACCCGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCa                                                                                                                                              <  3:113177/148‑1 (MQ=255)
           gCTACACCCGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCACGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAgcgc                                                                                                                                          <  8:307013/149‑1 (MQ=255)
             tACACCCGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGa                                                                                                                                        >  6:212855/1‑149 (MQ=255)
               cacCCGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACTGATGTCGGCCGCCTGTACCAGCGCGa                                                                                                                                        <  1:5799/147‑1 (MQ=255)
                acCCGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCt                                                                                                                                     >  7:146086/1‑149 (MQ=255)
                 cccGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCt                                                                                                                                     >  1:318182/1‑148 (MQ=255)
                  ccGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGa                                                                                                                                   <  4:91439/149‑1 (MQ=255)
                    gcgcCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGc                                                                                                                                 <  2:222151/149‑1 (MQ=255)
                    gcgcCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGAAGTCGGCCGCCTGTACCAGCGCGAGCTGAGc                                                                                                                                 >  4:34902/1‑149 (MQ=255)
                     cgccgccgcCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCg                                                                                                                                <  2:188406/149‑1 (MQ=255)
                        cgacgccgAACTTTCGCTGACCCAAAGCGCGGTCATCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTTAACCTGTTCCTCAGCTAGGGGAGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGa                                                                                                                               <  6:168109/144‑1 (MQ=255)
                          ccgccgAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGGCGGCCGACATCGGTCAGCTGTAcct                                                                                                                                              <  6:64643/130‑2 (MQ=38)
                          ccgccgAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGGCGGCCGACATCGGGCAGCTGTAcct                                                                                                                                              >  5:64643/1‑129 (MQ=255)
                          ccgccgAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCa                                                                                                                           <  7:179381/149‑1 (MQ=255)
                            gccgAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGa                                                                                                                                   >  4:148091/1‑139 (MQ=255)
                            gccgAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGa                                                                                                                                   <  3:148091/139‑1 (MQ=255)
                              cgAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTgg                                                                                                                       >  7:227364/1‑149 (MQ=255)
                              cgAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTgg                                                                                                                       <  3:49558/149‑1 (MQ=255)
                                aaCTGTCGCTGACCCAAAGCGCGGTCAGCTGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTgggg                                                                                                                     <  7:224263/149‑1 (MQ=255)
                                   tGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGgcg                                                                                                                   >  5:286730/1‑148 (MQ=255)
                                            cccAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAgc                                                                                                         >  5:212770/1‑149 (MQ=255)
                                            cccAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAgc                                                                                                         <  6:212770/149‑1 (MQ=255)
                                             ccAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAgcg                                                                                                        <  3:34902/149‑1 (MQ=255)
                                              cAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGc                                                                                                                                 <  4:285016/123‑1 (MQ=255)
                                              cAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGc                                                                                                                                 >  3:285016/1‑123 (MQ=255)
                                                  gcgcGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCAcc                                                                                                   >  6:188177/1‑149 (MQ=255)
                                                      ggTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGc                                                                                               <  6:104915/149‑1 (MQ=255)
                                                        tCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCAGCCTGTACCAGCGCGAGCTGAGCGAGGCACTCGGGCGCATCCGCAGCGCCACCTTGCAg                                                                                             >  8:296704/1‑149 (MQ=255)
                                                            ccgccAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCa                                                                                          >  1:138553/1‑148 (MQ=255)
                                                            ccgccAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCAc                                                                                         >  8:194374/1‑149 (MQ=255)
                                                                cAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGc                                                                                     <  2:55985/149‑1 (MQ=255)
                                                                   gTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTa                                                                                  <  5:212855/149‑1 (MQ=255)
                                                                         gCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTa                                                                                  >  5:249127/1‑143 (MQ=255)
                                                                         gCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTa                                                                                  <  6:249127/143‑1 (MQ=255)
                                                                         gCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATc                                                                            >  8:239090/1‑149 (MQ=255)
                                                                         gCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATc                                                                            <  6:286730/149‑1 (MQ=255)
                                                                           ccTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCgg                                                                          <  5:188177/149‑1 (MQ=255)
                                                                           ccTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCgg                                                                          <  8:329837/149‑1 (MQ=255)
                                                                           ccTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCgg                                                                          <  7:332340/149‑1 (MQ=255)
                                                                               gAGCAGCAGTTGGTCCGGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGa                                                                        <  4:230793/147‑1 (MQ=255)
                                                                               gAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAg                                                                       >  5:172436/1‑148 (MQ=255)
                                                                                                 aCCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGc                                                    >  2:131144/1‑149 (MQ=255)
                                                                                                 aCCCTGTTCAGCCGCGAGGGGCGCAAGGTACAGCAGACAGATGTCGGCCGCCTATACAAGCGCGAGCAGAGCGAGGAACTGGGGCGCCTCCGCAGCGCCACCTTGCAGGCACTCGCCAATCAATCGGGAGTAGGCACCATGCGCCTgg                                                     >  5:98580/1‑148 (MQ=255)
                                                                                                   ccTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCa                                                  >  2:26748/1‑149 (MQ=255)
                                                                                                      gTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCAccc                                               >  1:204247/1‑149 (MQ=255)
                                                                                                       ttCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCCCCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCt                                              >  5:209573/1‑149 (MQ=255)
                                                                                                                         ctggTACAGGCGGCCGACATCGGTCAGCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACCTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGa                             >  7:266964/3‑148 (MQ=255)
                                                                                                                         cAGGTACAGCTGCCCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGa                             >  5:265218/1‑148 (MQ=255)
                                                                                                                         cAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGa                             >  2:243865/1‑148 (MQ=255)
                                                                                                                            gTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATg                         >  3:233983/1‑149 (MQ=255)
                                                                                                                             tACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATgg                        <  8:227364/149‑1 (MQ=255)
                                                                                                                              aCAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGc                       >  2:324805/1‑149 (MQ=255)
                                                                                                                               cAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCt                      <  7:194374/149‑1 (MQ=255)
                                                                                                                               cAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCt                      <  7:239090/149‑1 (MQ=255)
                                                                                                                                aGCTGAACGATGTCGGCCGCCTGTACCAGCGCCAGCAGAGAGAGGCACTGGGGCGCATCCGCCGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGCCACTTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATggatc                     >  4:360360/1‑146 (MQ=255)
                                                                                                                                 gCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGt                    >  6:246876/1‑149 (MQ=255)
                                                                                                                                  cTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGtt                   >  5:23594/1‑149 (MQ=255)
                                                                                                                                    gACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGc                 >  3:294274/1‑149 (MQ=255)
                                                                                                                                     aCCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCAcc                                                >  1:146521/1‑117 (MQ=255)
                                                                                                                                     aCCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCAcc                                                <  2:146521/117‑1 (MQ=255)
                                                                                                                                     aCCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGt                    >  1:118880/1‑145 (MQ=255)
                                                                                                                                     aCCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGt                    <  2:118880/145‑1 (MQ=255)
                                                                                                                                     aCCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGcc                >  1:226685/1‑149 (MQ=255)
                                                                                                                                      ccGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCg               <  5:246876/149‑1 (MQ=255)
                                                                                                                                       cGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGcc                >  3:222293/1‑147 (MQ=255)
                                                                                                                                       cGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGa              <  1:158861/149‑1 (MQ=255)
                                                                                                                                       cGATGGCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGa              <  2:138553/149‑1 (MQ=255)
                                                                                                                                         aTGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTa                                                                     <  8:92223/92‑1 (MQ=255)
                                                                                                                                         aTGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTa                                                                     >  7:92223/1‑92 (MQ=255)
                                                                                                                                         aTGGCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAgg            <  2:318182/149‑1 (MQ=255)
                                                                                                                                              ggCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCa       >  5:182127/1‑149 (MQ=255)
                                                                                                                                              ggCCGCCGGTACCAGCGCGATCTGAGAGAGGCACTGGGGCGCATCCGCAGCGACACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACGCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGc        <  4:233983/148‑1 (MQ=255)
                                                                                                                                                   ccTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGt  >  2:173646/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGGCATGAAAGCTACACCCGCGCCGCCGCCGAACTGTCGCTGACCCAAAGCGCGGTCAGCCGCCAGGTACAGGCCCTGGAGCAGCAGTTGGGCCTGACCCTGTTCCGCCGCGAGGGGCGCCAGGTACAGCTGACCGATGTCGGCCGCCTGTACCAGCGCGAGCTGAGCGAGGCACTGGGGCGCATCCGCAGCGCCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGT  >  NC_002947/5108323‑5108618

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