Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,443,315 G→C 6.4% P489A (CCG→GCG)  PP_4781 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,443,3150GC 6.4% 199.1 / 3.5 62P489A (CCG→GCG) PP_4781sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base C (2/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 = 9.46e-01

GTGCGCCAGCTTCTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGA  >  NC_002947/5443169‑5443462
                                                                                                                                                  |                                                                                                                                                   
gTGCGCCAGCTTCTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCt                                                                                                                                                   >  4:166974/1‑149 (MQ=255)
            cTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAgg                                                                                                                                        <  2:17622/148‑1 (MQ=255)
               cccGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGCCTGCGGGCGCAgccgg                                                                                                                                     >  1:301640/1‑148 (MQ=255)
                 cGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGcggccg                                                                                                                                  >  4:30184/1‑149 (MQ=255)
                 cGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGCCTGCGGGCGCAgccggcac                                                                                                                                  <  2:301640/149‑3 (MQ=255)
                       cGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGc                                                                                                                            <  4:308197/149‑1 (MQ=255)
                       cGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGc                                                                                                                            <  6:115500/149‑1 (MQ=255)
                        gAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTg                                                                                                                             >  1:359994/1‑147 (MQ=255)
                              cccGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCg                                                                                                                      >  1:142546/1‑148 (MQ=255)
                              cccGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGCCTGCGGGCGCAgccggca                                                                                                                                   <  4:272178/135‑2 (MQ=255)
                              cccGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGCCTGCGGGCGCAgccggca                                                                                                                                   >  3:272178/1‑134 (MQ=255)
                                cGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGAc                                                                                                                   <  6:143002/149‑1 (MQ=255)
                                 gTCCGGCAGGTCGGCCGGAATGACCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACa                                                                                                                  <  2:186684/149‑1 (MQ=255)
                                    cGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAggg                                                                                                               <  5:103609/149‑1 (MQ=255)
                                    cGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAggg                                                                                                               >  6:1989/1‑149 (MQ=255)
                                    cGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGG‑CTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGt                                                                                                              <  5:201820/149‑1 (MQ=255)
                                     ggCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGcggccg                                                                                                                                  >  6:373206/1‑129 (MQ=255)
                                     ggCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGcggccg                                                                                                                                  <  5:373206/129‑1 (MQ=255)
                                       cAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTaa                                                                                                            <  3:166974/149‑1 (MQ=255)
                                                       tCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCa                                                                                             >  8:49706/1‑148 (MQ=255)
                                                        ccAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCata                                                                                           >  5:194155/1‑149 (MQ=255)
                                                              tCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTaa                                                                                                            >  6:62195/1‑126 (MQ=255)
                                                              tCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTaa                                                                                                            <  5:62195/126‑1 (MQ=255)
                                                                aGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCcatcat                                                                                   >  7:165210/1‑149 (MQ=255)
                                                                  gTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGc                                                                                 <  3:185162/149‑1 (MQ=255)
                                                                    aCAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCt                                                                               >  6:79194/1‑149 (MQ=255)
                                                                      aGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTg                                                                             <  2:142546/149‑1 (MQ=255)
                                                                          cTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTggg                                                                           >  2:294965/1‑147 (MQ=255)
                                                                               cccAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGAccc                                                                    >  4:212977/1‑149 (MQ=255)
                                                                                 cAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCca                                                                  >  4:213843/1‑149 (MQ=255)
                                                                                 cAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCca                                                                  >  1:286565/1‑149 (MQ=255)
                                                                                  aGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCcac                                                                 <  5:1989/149‑1 (MQ=255)
                                                                                     gcCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGc                                                              >  2:37522/1‑149 (MQ=255)
                                                                                      ccTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGcc                                                             <  2:286565/149‑1 (MQ=255)
                                                                                       cTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGCCCCCACAGCCg                                                            >  8:62008/1‑149 (MQ=255)
                                                                                        tGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgt                                                           >  3:246339/1‑149 (MQ=255)
                                                                                           ccGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgtgt                                                         >  1:21155/1‑148 (MQ=255)
                                                                                            cGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGa                                                       >  3:276251/1‑149 (MQ=255)
                                                                                             ggCATAACCGGCGCGCAAACCGTCCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACATCCGTGTGAt                                                      <  3:286360/149‑1 (MQ=255)
                                                                                             ggCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGAt                                                      <  6:194155/149‑1 (MQ=255)
                                                                                              gcttAACGGGCGCGAAACCCTTGCACCTGGGCGCGATGGCTTAAGGTTCCGGGCTGCGCCGGATGGGGGCGGGCTGGATTTTCGGCAAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATg                                                     <  1:294965/146‑1 (MQ=255)
                                                                                              gCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATg                                                     <  7:49706/149‑1 (MQ=255)
                                                                                                    ccGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACt                                                >  3:43749/1‑148 (MQ=255)
                                                                                                    ccGGCGCGCAAACCCTGCAACAGCGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTg                                               >  8:205101/1‑149 (MQ=255)
                                                                                                           gcAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCata                                                                                           >  8:186396/1‑98 (MQ=255)
                                                                                                           gcAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCata                                                                                           <  7:186396/98‑1 (MQ=255)
                                                                                                           gcAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCAtt                                        >  5:332951/1‑149 (MQ=255)
                                                                                                            cAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCAttt                                       <  1:37522/149‑1 (MQ=255)
                                                                                                             aaaCCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGCGCAGACCCCACAGCCGTGTGATGGTACTGAGCCAttt                                       >  4:341530/1‑148 (MQ=255)
                                                                                                                cccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGc                                   <  3:30184/149‑1 (MQ=255)
                                                                                                                cccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGc                                   <  7:205101/149‑1 (MQ=255)
                                                                                                                  cTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                 >  1:196/1‑149 (MQ=255)
                                                                                                                  cTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                 <  2:196/149‑1 (MQ=255)
                                                                                                                   tGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGcccc                                >  4:361206/1‑149 (MQ=255)
                                                                                                                   tGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGcccc                                >  1:45168/1‑149 (MQ=255)
                                                                                                                        cAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCg                           <  3:212977/149‑1 (MQ=255)
                                                                                                                                  aGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGCCCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGtgtg                 >  8:180800/1‑149 (MQ=255)
                                                                                                                                   ggCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGc                <  2:45168/149‑1 (MQ=255)
                                                                                                                                    gCTGAACGTTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGc                <  4:43749/148‑1 (MQ=255)
                                                                                                                                    gCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGCCAGGGggtgaa               >  6:256572/1‑147 (MQ=255)
                                                                                                                                     cTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCAt              <  8:165210/149‑1 (MQ=255)
                                                                                                                                      tGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCAt              >  3:275234/1‑148 (MQ=255)
                                                                                                                                                 ggCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcga  <  5:256572/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTGCGCCAGCTTCTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGA  >  NC_002947/5443169‑5443462

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