Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 498,339 A→G 6.3% G737G (GGA→GGG)  PP_0409 → two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947498,3390AG 6.3% 180.3 / 2.5 63G737G (GGA→GGG) PP_0409two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base G (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

GGTGCAAGTACCCGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCTGGGGGGGAGCATGACGCTGGAAAGCGAGCCGGGCGAAGGCACCA  >  NC_002947/498198‑498487
                                                                                                                                              |                                                                                                                                                    
ggTGCAAGTACCCGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTgc                                                                                                                                               >  4:356927/1‑149 (MQ=255)
 gTGCAAGTACCCGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAgct                                                                                                                                                 <  6:304993/146‑3 (MQ=255)
 gTGCAAGTACCCGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTgcg                                                                                                                                              <  3:84671/149‑1 (MQ=255)
 gTGCAAGTACCCGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTgc                                                                                                                                               <  3:33313/148‑1 (MQ=255)
     aaGTACCCGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCAt                                                                                                                                           >  8:330027/1‑148 (MQ=255)
     aaGTACCCGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCAt                                                                                                                                           >  6:177694/1‑148 (MQ=255)
       gTACCCGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGAc                                                                                                                                        >  7:155843/1‑149 (MQ=255)
          cccGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACgg                                                                                                                                      <  5:177694/148‑1 (MQ=255)
          cccGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACgg                                                                                                                                      <  7:330027/148‑1 (MQ=255)
          cccGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGAAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGt                                                                                                                                     <  4:47740/149‑1 (MQ=255)
           ccGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGt                                                                                                                                     >  8:382124/1‑148 (MQ=255)
             gATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTg                                                                                                                                                >  5:261751/1‑135 (MQ=255)
             gATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTg                                                                                                                                                <  6:261751/135‑1 (MQ=255)
              aTAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGt                                                                                                                                 <  8:155843/149‑1 (MQ=255)
                          cATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGcc                                                                                                                       >  5:353086/1‑147 (MQ=255)
                          cATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCaa                                                                                                                     <  5:121661/149‑1 (MQ=255)
                           aTTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAg                                                                                                                    >  6:128909/1‑149 (MQ=255)
                            ttGAGTGATGCCAAGGCCACCGGGCTGTTTCATATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAgg                                                                                                                   <  5:368691/149‑1 (MQ=255)
                                gTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCATGGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCtt                                                                                                                >  1:73401/1‑149 (MQ=255)
                                gTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCttt                                                                                                               >  4:7268/1‑149 (MQ=255)
                                     gCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCcg                                                                                                          >  1:410483/1‑149 (MQ=255)
                                      ccAAGGCCACCGGGCTGTTTCGTATCATCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCcgc                                                                                                         >  5:181254/1‑149 (MQ=255)
                                       cAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCcgc                                                                                                         >  1:176850/1‑148 (MQ=255)
                                                       tttCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGccc                                                                                          <  2:73401/147‑1 (MQ=255)
                                                         tCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCAcc                                                                                       <  2:410483/148‑1 (MQ=255)
                                                           gTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCa                                                                                    <  3:356927/149‑1 (MQ=255)
                                                                ctccA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTgg                                                                               >  5:5620/1‑149 (MQ=255)
                                                                 tccA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGc                                                                              >  2:188489/1‑149 (MQ=255)
                                                                 tccA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGc                                                                              >  1:245610/1‑149 (MQ=255)
                                                                   cA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGGCAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGa                                                                                                       <  6:183126/122‑1 (MQ=255)
                                                                   cA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGGCAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGa                                                                                                       <  6:189084/122‑1 (MQ=255)
                                                                   cA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGGCAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGa                                                                                                       >  5:189084/1‑122 (MQ=255)
                                                                   cA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGGCAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGa                                                                                                       >  5:183126/1‑122 (MQ=255)
                                                                   cA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGcc                                                                             <  1:188489/148‑1 (MQ=255)
                                                                   cA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGc                                                                              <  2:245610/147‑1 (MQ=255)
                                                                    a‑ggaggCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGtgat                                                                                                                              <  8:63923/98‑1 (MQ=255)
                                                                         ggCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCtggtgg                                                                       <  7:382124/149‑1 (MQ=255)
                                                                         ggCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCtggtgg                                                                       <  3:360026/149‑1 (MQ=255)
                                                                           cgcTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTggtgggt                                                                     <  3:7268/149‑1 (MQ=255)
                                                                                 ccAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTgcgc                                                               <  8:26328/149‑1 (MQ=255)
                                                                                 ccAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTgcgc                                                               <  8:150383/149‑1 (MQ=255)
                                                                                     tgtgATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGc                                                           >  4:365358/1‑149 (MQ=255)
                                                                                         aTGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCggg                                                        <  2:176850/148‑1 (MQ=255)
                                                                                         aTGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGt                                                       >  8:137240/1‑149 (MQ=255)
                                                                                         aTGCGTCACGCCCAGGCGAACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCGGGTGGGTGTGCGCGAGCGGGt                                                       >  8:178595/1‑149 (MQ=255)
                                                                                           gCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGc                                                     >  1:375719/1‑149 (MQ=255)
                                                                                            cGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGc                                                     <  5:128909/148‑1 (MQ=255)
                                                                                            cGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGc                                                     <  2:375719/148‑1 (MQ=255)
                                                                                             gTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGgtgt                                                                   <  5:20641/133‑1 (MQ=255)
                                                                                             gTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGgtgt                                                                   >  6:20641/1‑133 (MQ=255)
                                                                                              tCACTCCCAGGCGCACAGCTTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGa                                                  <  6:181254/149‑1 (MQ=255)
                                                                                               cACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGAt                                                 >  7:234162/1‑149 (MQ=255)
                                                                                                  gCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCt                                              >  4:153541/1‑149 (MQ=255)
                                                                                                    ccAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCTgg                                            <  6:353086/149‑1 (MQ=255)
                                                                                                          gcACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCTgggggg                                        >  6:135910/1‑147 (MQ=255)
                                                                                                            acaGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCTGGGGGGGAg                                     >  1:69403/1‑148 (MQ=255)
                                                                                                            acaGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCTGGGGGGGAg                                     >  1:298035/1‑148 (MQ=255)
                                                                                                                        ttGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGc                                                                                                     >  4:101829/1‑72 (MQ=255)
                                                                                                                        ttGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGc                                                                                                     <  3:101829/72‑1 (MQ=255)
                                                                                                                             cTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCTGGGGGGGAGCATGACGCTGGAAAGc                     >  5:93783/1‑147 (MQ=255)
                                                                                                                                gTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCTGGGGGGGAGCATGACGCTGGAAAGCGAGcc                >  4:345986/1‑149 (MQ=255)
                                                                                                                                      gAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCTGGGGGGGAGCATGACGCTGGAAAGCGAGCCGGGCg           >  3:75100/1‑148 (MQ=255)
                                                                                                                                         ggCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGACCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGGTGCTGGGGGGGAGAATGACGCTGGAAAGCGAGCCGGGCGAAgg       >  3:86258/1‑149 (MQ=255)
                                                                                                                                              aCAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCTGGGGGGGAGCATGACGCTGGAAAGCGAGCCGGGCGAAGGCACCa  <  3:345986/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGTGCAAGTACCCGATAATCTGCCAGCATTGAGTGATGCCAAGGCCACCGGGCTGTTTCGTATCCTCCA‑GGAGGCGCTGACCAATGTGATGCGTCACGCCCAGGCGCACAGCGTAGAGATTGAGCTGGTACGTGAGGGCGGACAGTTGCGCATGACGGTCAGTGATGATGGCCAAGGCTTTTGCCGCGAGCAGACCCGGCCCACCTCATTTGGCCTGGTGGGTGTGCGCGAGCGGGTGCTGATGCTGGGGGGGAGCATGACGCTGGAAAGCGAGCCGGGCGAAGGCACCA  >  NC_002947/498198‑498487

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