Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,914,433 G→C 8.1% R118R (CGC→CGG)  PP_3454 ← transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,914,4330GC 8.1% 192.8 / 6.8 62R118R (CGC→CGG) PP_3454transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base C (3/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01

GCCAGCACCACCAGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCTGCGACTCCT  >  NC_002947/3914287‑3914581
                                                                                                                                                  |                                                                                                                                                    
gcCAGCACCACCAGCAGGTTGTACTCACCGCTGGACAGCTCCCCCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTgcggcg                                                                                                                                                    >  3:224352/1‑149 (MQ=255)
     caccacCAGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCATTCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCgcagca                                                                                                                                               <  1:205196/149‑1 (MQ=255)
       ccaccaGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCagcagc                                                                                                                                             <  1:114910/149‑1 (MQ=255)
         accaGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACGCTGCTGCGCCGCAGCAgtgt                                                                                                                                           <  5:50459/149‑1 (MQ=255)
          ccaGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGCCGCAGCAgtgtc                                                                                                                                          <  7:153821/149‑2 (MQ=255)
           cagcagGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTgc                                                                                                                                         >  4:182787/1‑149 (MQ=255)
           cagcagGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTgc                                                                                                                                         >  4:149131/1‑149 (MQ=255)
              cagGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCa                                                                                                                                      >  2:108825/1‑149 (MQ=255)
                 gTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCAccc                                                                                                                                   >  8:233830/1‑149 (MQ=255)
                        tCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCcag                                                                                                                            >  3:109310/1‑149 (MQ=255)
                         cACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCcagc                                                                                                                           >  8:80993/1‑149 (MQ=255)
                          aCCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCcag                                                                                                                            >  8:187758/1‑147 (MQ=255)
                                 gACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCAcgcg                                                                                                                   >  3:33689/1‑149 (MQ=255)
                                  aCAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTgc                                                                                                                                         <  4:65769/126‑1 (MQ=255)
                                  aCAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTgc                                                                                                                                         >  3:65769/1‑126 (MQ=255)
                                   cAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCATCCGGGCCAGCAGCACGCGTg                                                                                                                 <  4:141046/149‑1 (MQ=255)
                                     gCTCCACCTCCTCGCCGCGCCAACTCCCCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGc                                                                                                               >  5:207622/1‑149 (MQ=255)
                                      cTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCATGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCt                                                                                                              >  8:191637/1‑149 (MQ=255)
                                       tCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGc                                                                                                               <  7:252839/147‑1 (MQ=255)
                                       tCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGc                                                                                                               >  8:252839/1‑147 (MQ=255)
                                          aCCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTca                                                                                                            >  6:132040/1‑147 (MQ=255)
                                             tcctcGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGc                                                                                                               <  7:113573/141‑1 (MQ=255)
                                             tcctcGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGc                                                                                                               >  8:113573/1‑141 (MQ=255)
                                                  gccgcGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTgg                                                                                                  >  5:116032/1‑149 (MQ=255)
                                                     gcgcCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCa                                                                                               <  6:116032/149‑1 (MQ=255)
                                                     gcgcCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCa                                                                                               <  5:132040/149‑1 (MQ=255)
                                                      cgcCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCAc                                                                                              <  3:149131/149‑1 (MQ=255)
                                                      cgcCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCAc                                                                                              <  4:109310/149‑1 (MQ=255)
                                                       gcCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACg                                                                                             <  4:213243/149‑1 (MQ=255)
                                                          aaCTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAg                                                                                          >  3:171358/1‑149 (MQ=255)
                                                           aCTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCAcg                                                                                                                     >  1:116131/1‑121 (MQ=255)
                                                           aCTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCAcg                                                                                                                     <  2:116131/121‑1 (MQ=255)
                                                               aCCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAgtcgtc                                                                                     >  2:216022/1‑149 (MQ=255)
                                                                    gcgcTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTg                                                                                                                                          >  6:153406/1‑91 (MQ=255)
                                                                    gcgcTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTg                                                                                                                                          <  5:153406/91‑1 (MQ=255)
                                                                          ggCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCAttt                                                                          >  4:41160/1‑149 (MQ=255)
                                                                           gCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTc                                                                         <  7:233830/149‑1 (MQ=255)
                                                                                  cGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGcc                                                                   >  3:155545/1‑148 (MQ=255)
                                                                                     tACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAg                                                                >  7:19876/1‑148 (MQ=255)
                                                                                        gCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCAcc                                                            >  1:82635/1‑149 (MQ=255)
                                                                                          aGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCTGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTg                                                          <  3:182787/149‑1 (MQ=255)
                                                                                           ggCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTgg                                                         <  4:33689/149‑1 (MQ=255)
                                                                                              caccaATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCa                                                                             <  1:105604/126‑1 (MQ=255)
                                                                                              caccaATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCa                                                      >  4:128430/1‑149 (MQ=255)
                                                                                              caccaATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCa                                                      >  3:89081/1‑149 (MQ=255)
                                                                                                 caATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATgg                                                   <  7:191637/149‑1 (MQ=255)
                                                                                                   atcaggatcaggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCt                                                 <  4:171358/149‑1 (MQ=255)
                                                                                                   atcaggatcaggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCt                                                 <  7:187758/149‑1 (MQ=255)
                                                                                                      aggatcaggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTAtc                                              <  4:89081/149‑1 (MQ=255)
                                                                                                          tcaggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCt                                                           <  7:241557/132‑1 (MQ=255)
                                                                                                          tcaggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCt                                                           >  8:241557/1‑132 (MQ=255)
                                                                                                             ggTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCgg                                       >  2:254640/1‑149 (MQ=255)
                                                                                                               tCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGc                                     >  4:162816/1‑149 (MQ=255)
                                                                                                                 gTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGcc                                    <  2:82635/148‑1 (MQ=255)
                                                                                                                  ttGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGt                                  <  1:108825/149‑1 (MQ=255)
                                                                                                                    gTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATgg                                                   <  2:10782/130‑1 (MQ=255)
                                                                                                                    gTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATgg                                                   >  1:10782/1‑130 (MQ=255)
                                                                                                                                   gTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCa                 <  4:155545/149‑1 (MQ=255)
                                                                                                                                    ttGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGc                                     <  4:115401/128‑1 (MQ=255)
                                                                                                                                       acgctgctgCGCCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGt             >  5:92575/4‑149 (MQ=255)
                                                                                                                                         gctgctgCGCCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCt           >  4:84620/2‑149 (MQ=255)
                                                                                                                                         acTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCt           >  6:226201/1‑149 (MQ=255)
                                                                                                                                          ctgctgCGCCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCTg          >  2:99724/1‑149 (MQ=255)
                                                                                                                                                  gcgcAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCTGCGACTCCt  <  8:19876/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCAGCACCACCAGCAGGTTGTACTCACCGCTGGACAGCTCCACCTCCTCGCCGCGCCAACTCACCGTGCGCTCGGCCAGGTCGATACGCAGGCCACCAATCAGGATCAGGTCGTTGTCCAGGCGCGGTTCGTTGACACTGCTGCGGCGCAGCAGCGTGCGCACCCGGGCCAGCAGCACGCGTGGCTCACAGGGTTTGGTCACGTAGTCGTCGGCCCCCATTTCCAGGCCTAGCACCTGGTCATGGCTATCGTCGCGGGCCGTCAGCATGAGGATTGGCAGGGTCTGCGACTCCT  >  NC_002947/3914287‑3914581

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