Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 385,475 A→C 9.2% intergenic (+15/+20) PP_0320 → / ← ltaE methyl‑accepting chemotaxis transducer/low specificity l‑threonine aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947385,4750AC 9.2% 157.2 / 3.2 59intergenic (+15/+20)PP_0320/ltaEmethyl‑accepting chemotaxis transducer/low specificity l‑threonine aldolase
Reads supporting (aligned to +/- strand):  ref base A (30/23);  new base C (2/4);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCCC  >  NC_002947/385333‑385622
                                                                                                                                              |                                                                                                                                                        
cGGTGGAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCAcggcgg                                                                                                                                                    <  3:3984/149‑1 (MQ=255)
     cAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCg                                                                                                                                               >  6:67935/1‑149 (MQ=255)
     cAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCg                                                                                                                                               >  3:130646/1‑149 (MQ=255)
      aGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCCGTGCGGcccctg                                                                                                                                              <  1:216534/149‑3 (MQ=255)
      aGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCgg                                                                                                                                              <  3:62522/149‑1 (MQ=255)
      aGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCgg                                                                                                                                              >  4:176256/1‑149 (MQ=255)
      aGGGGATGTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCCGGGGCCGCCGTGCGGcccctg                                                                                                                                              <  8:329899/149‑3 (MQ=255)
       ggggATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCCGGGGCCGCCGTGCGGCCCCTGc                                                                                                                                             >  6:226920/1‑149 (MQ=255)
          gATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCCGGGGCCGCCGTGCGGCCCCTGcgat                                                                                                                                          <  5:226920/149‑4 (MQ=255)
           aTTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCtctc                                                                                                                                         <  5:87335/149‑1 (MQ=255)
             ttCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGc                                                                                                                                       >  5:134860/1‑149 (MQ=255)
             ttCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGc                                                                                                                                       >  3:9782/1‑149 (MQ=255)
              tCCTATTTGGCGCGGGAGACGGCGGCGGAGGTGCAAGGGGGCCGCGAGGAGTGCCAGGCGTTGCGGGGGTTGGCTGATGACCTGGCGCGGCAGAGGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGc                                                                                                                                       <  1:2455/148‑1 (MQ=255)
                cGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGCCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGCCGGCCCCGGCTCTCGCt                                                                                                                                      >  7:137077/1‑147 (MQ=255)
                cGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCt                                                                                                                                      >  4:215391/1‑147 (MQ=255)
                                aCGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTg                                                                                                                    >  1:214259/1‑149 (MQ=255)
                                  ggcggcggAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCg                                                                                                                  >  6:121191/1‑149 (MQ=255)
                                           ggTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATccgc                                                                                                         <  1:256398/149‑1 (MQ=255)
                                              gCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGCGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGATTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAAAGATCGTGCGGATATCagccgc                                                                                                      >  4:271824/1‑149 (MQ=255)
                                                aaGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCcgccgc                                                                                                      >  7:176815/1‑147 (MQ=255)
                                                 aGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGCCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGGGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAg                                                                                                   <  3:271824/149‑1 (MQ=255)
                                                         gcgAGGAGTGCCAGGCGTTGCGTGAGATGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTAAGGCTCTAGATCGCAGGGGCCGCCCGGCGGCCCCGGCTCTCGCTCCGCCAGCAATCAACGTGCGGAAAACCGCCGCCAGCT‑‑CA‑CGCa                                                                                        >  6:223799/1‑149 (MQ=255)
                                                             ggagTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCcg                                                                                    >  1:73938/1‑149 (MQ=255)
                                                                gTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCgcgc                                                                                 <  5:67935/149‑1 (MQ=255)
                                                                   ccAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGct‑tc                                                                             >  2:314562/1‑149 (MQ=255)
                                                                     aGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGct‑tctt                                                                           >  1:9507/1‑149 (MQ=255)
                                                                       gctttGCGGGGGTTGGCTGATCCCGGGGCGGGGCAGATGGGTGGGCTCTGCCTCTAGATCCCCGGGCCCCCCGGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGct‑tctt                                                                           <  3:176256/144‑1 (MQ=255)
                                                                       gcgtggcgtgGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCg                                                                         <  4:338797/149‑1 (MQ=255)
                                                                       gCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCg                                                                         <  4:130646/149‑1 (MQ=255)
                                                                       gCGTTGCGTGGGGTGGCTGATGACCTGGCGCGGCAGAGGGGGGGGTGCCGGCTCTAGATCGCAGGGCCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCg                                                                         <  7:166874/149‑1 (MQ=255)
                                                                            gCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGcc                                                                                                                                 >  2:272539/1‑92 (MQ=255)
                                                                            gCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGcc                                                                                                                                 <  1:272539/92‑1 (MQ=255)
                                                                             cgggggTGGGCTTCTGACCTGGCGCGGACTATGGGTGGGTTCAGGCTCTCGATCGCCGGGGCCGCAGGGGGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATcc                                                                   <  4:9782/146‑1 (MQ=255)
                                                                              gTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATccc                                                                  <  3:215391/149‑1 (MQ=255)
                                                                                      gCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑ca                                                          <  5:182308/148‑1 (MQ=255)
                                                                                      gCTGATGACCTGGCGCGGCAGATGGGGGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCTCT‑TCTTCGGTATCCCACGAG‑cac                                                         <  3:345312/149‑1 (MQ=255)
                                                                                        tgatgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTGTCCCACGAG‑CACAt                                                       >  2:62237/1‑149 (MQ=255)
                                                                                         gatgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACAt                                                       >  3:278806/1‑148 (MQ=255)
                                                                                           tgcccgggCGCGGCAGATGGGTGGGTGCCGGCTCTAGATCGCAGGGCCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGaa                                                    <  3:249552/143‑1 (MQ=255)
                                                                                            gaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGaa                                                    <  3:102752/148‑1 (MQ=255)
                                                                                               cTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCCCCCGCGCT‑TCGTCGGTATCCCACGAG‑CACATGACgcgc                                                >  4:295997/1‑149 (MQ=255)
                                                                                                  gcgcGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAgcgcgcg                                             >  3:40037/1‑149 (MQ=255)
                                                                                                       gCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCAcc                                        <  4:40037/149‑1 (MQ=255)
                                                                                                          gATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCt                                     >  4:316975/1‑149 (MQ=255)
                                                                                                               gtgggtTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGa                                >  7:302552/1‑149 (MQ=255)
                                                                                                                  ggtTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATa                              >  8:283912/1‑148 (MQ=255)
                                                                                                                   gtTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATa                              >  4:195071/1‑147 (MQ=255)
                                                                                                                    ttCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATAAAg                           >  8:314708/1‑149 (MQ=255)
                                                                                                                          ctctAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCAGCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATAAAGGta                        >  1:81362/1‑146 (MQ=255)
                                                                                                                                aTCGCAGGGGCCGCACAGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCt‑tcttc                                                                          >  2:20047/1‑91 (MQ=255)
                                                                                                                                aTCGCAGGGGCCGCACAGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCt‑tcttc                                                                          <  1:20047/91‑1 (MQ=255)
                                                                                                                                  cGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTGGGAGCCCACCCGCGCTGAGATCGGTAACCCACGAGCCACA‑GAAGCGCGCCCCCCCGCGGCCGATAAAGGTATAGAAgc                 >  1:116402/1‑149 (MQ=255)
                                                                                                                                   gCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCa             <  7:283912/148‑1 (MQ=255)
                                                                                                                                    cAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGc           >  1:132659/1‑149 (MQ=255)
                                                                                                                                     aGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGGTGCGGAGAAAGGTATAGAAGCGCCAGcc          >  4:244010/1‑149 (MQ=255)
                                                                                                                                      ggggCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGccc         >  4:12526/1‑149 (MQ=255)
                                                                                                                                        ggCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCtt       <  2:214259/149‑1 (MQ=255)
                                                                                                                                          ccGCAGGGCGGTCCCGGCTCTCGCTCAGCCTCCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGc     <  5:223799/149‑1 (MQ=255)
                                                                                                                                             cACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGcccc  <  3:195071/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                        
CGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCT‑‑CA‑CGCACCCGCGCT‑TCTTCGGTATCCCACGAG‑CACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCCC  >  NC_002947/385333‑385622

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