Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,360,289 G→C 9.6% A278G (GCC→GGC)  PP_3835 ← two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,360,2890GC 9.6% 210.8 / 11.5 73A278G (GCC→GGC) PP_3835two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base C (2/5);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

ATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCGTTAA  >  NC_002947/4360148‑4360436
                                                                                                                                                  |                                                                                                                                                    
aTGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACg                                                                                                                                              <  3:93127/149‑1 (MQ=255)
aTGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACg                                                                                                                                              <  3:246321/149‑1 (MQ=255)
  gcgcGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACggg                                                                                                                                            >  2:294676/1‑149 (MQ=255)
   cgcgGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACGGGc                                                                                                                                           <  1:17893/149‑1 (MQ=255)
   cgcgGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GACAGGCTGATGCGGCCATGCCCGTGATGCCATCGACggg                                                                                                                                            <  8:95728/148‑1 (MQ=255)
     cgGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCCGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACGGGCAt                                                                                                                                         <  8:254052/149‑1 (MQ=255)
            aTCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACGGGCAtg                                                                                                                                        >  2:320934/1‑142 (MQ=255)
            aTCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACGGGCAtg                                                                                                                                        <  1:320934/143‑2 (MQ=255)
               tctGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAg                                                                                                                               >  8:122928/1‑149 (MQ=255)
               tctGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAg                                                                                                                               >  2:112403/1‑149 (MQ=255)
                 tgtgGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGt                                                                                                                             >  4:171019/1‑149 (MQ=255)
                   tgGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTg                                                                                                                           <  6:159238/149‑1 (MQ=255)
                      gTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAg                                                                                                                         >  3:205182/1‑148 (MQ=255)
                      gTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGc                                                                                                                        <  5:46244/149‑1 (MQ=255)
                            cTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACAGGAGGACCACATCGACA‑‑GTAAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCaa                                                                                                                      >  1:107605/1‑148 (MQ=255)
                            cTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAacac                                                                                                                   >  2:94118/1‑148 (MQ=255)
                               cgccgcTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCGCATCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACACCTTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACAcc                                                                                                                  <  1:275659/148‑1 (MQ=255)
                                       tCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTg                                                                                                       >  4:325523/1‑149 (MQ=255)
                                         ttCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGcc                                                                                                     >  6:30133/1‑149 (MQ=255)
                                           ccACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCg                                                                                                   >  7:47091/1‑149 (MQ=255)
                                             aCCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGCCCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACGGGCAtggc                                                                                                                                      >  6:186782/1‑109 (MQ=38)
                                             aCCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACGGGCAtggc                                                                                                                                      <  5:186782/112‑4 (MQ=38)
                                                tCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAAATTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCAGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATc                                                                                              >  8:269671/1‑149 (MQ=255)
                                                 cGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAac                                                                                                                     >  4:161665/1‑125 (MQ=255)
                                                 cGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAac                                                                                                                     <  3:161665/125‑1 (MQ=255)
                                                  gAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATc                                                                                              <  5:187878/147‑1 (MQ=255)
                                                   aaGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCAc                                                                                           >  2:244368/1‑149 (MQ=255)
                                                         aCGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGGGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGAt                                                                                     <  6:180898/149‑1 (MQ=255)
                                                                 gcccttcagcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCAc                                                                                           >  6:123021/1‑135 (MQ=255)
                                                                 gcccttcagcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCAc                                                                                           <  5:123021/135‑1 (MQ=255)
                                                                 gcccttcagcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACgg                                                                               >  2:281810/1‑147 (MQ=255)
                                                                  cccttcagcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGc                                                                              <  1:294676/147‑1 (MQ=255)
                                                                  cccttcagcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAg                                                                            >  1:173089/1‑149 (MQ=255)
                                                                  cccttcagcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑CTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAg                                                                            >  6:283707/1‑149 (MQ=255)
                                                                   ccttcagcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAgg                                                                           >  2:99116/1‑149 (MQ=255)
                                                                     ttcagcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTc                                                                         >  2:317370/1‑149 (MQ=255)
                                                                      tcagcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGATcc                                                                        >  5:76700/1‑149 (MQ=255)
                                                                        agcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACCTCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCg                                                                       <  4:205182/148‑1 (MQ=255)
                                                                        agcccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCgg                                                                      <  5:283707/149‑1 (MQ=255)
                                                                          cccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGGTGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGtt                                                                    <  3:268955/149‑1 (MQ=255)
                                                                          cccttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGtt                                                                    <  2:115534/149‑1 (MQ=255)
                                                                             ttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATc                                                                                                               <  5:191207/103‑1 (MQ=255)
                                                                             ttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATc                                                                                                               >  6:191207/1‑103 (MQ=255)
                                                                             ttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCa                                                                             >  1:114646/1‑137 (MQ=255)
                                                                             ttcaACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCa                                                                             <  2:114646/137‑1 (MQ=255)
                                                                               caACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgca                                                                 <  3:325523/147‑1 (MQ=255)
                                                                               caACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgca                                                                 <  3:171019/147‑1 (MQ=255)
                                                                                     cGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAacac                                                          <  1:317370/148‑1 (MQ=255)
                                                                                      gCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATTTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCt                                                        <  1:94118/149‑1 (MQ=255)
                                                                                          aacaac‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCaa                                                    <  5:30133/149‑1 (MQ=255)
                                                                                          aacaac‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCaa                                                    <  7:269671/149‑1 (MQ=255)
                                                                                              ac‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTc                                                 >  8:23077/1‑148 (MQ=255)
                                                                                               c‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCaa                                               >  5:207971/1‑149 (MQ=255)
                                                                                               c‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCa                                                >  1:187591/1‑148 (MQ=255)
                                                                                                     gACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAacaaca                                                                                                                                <  1:162092/65‑1 (MQ=255)
                                                                                                     gACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAacaaca                                                                                                                                >  2:162092/1‑65 (MQ=255)
                                                                                                        cacaTCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGc                                         >  2:238042/1‑149 (MQ=255)
                                                                                                         acaTCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTc                                                                         >  2:216261/1‑116 (MQ=255)
                                                                                                         acaTCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTc                                                                         <  1:216261/116‑1 (MQ=255)
                                                                                                         acaTCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCt                                        <  2:173089/149‑1 (MQ=255)
                                                                                                         acaTCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCt                                        >  3:338075/1‑149 (MQ=255)
                                                                                                           aTCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTg                                      <  8:47091/149‑1 (MQ=255)
                                                                                                            tCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAAGAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTgg                                     <  1:281810/149‑1 (MQ=255)
                                                                                                                    gTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGc                               <  2:187591/149‑1 (MQ=255)
                                                                                                                        ggCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCaa                                                                                                                      <  5:354060/58‑1 (MQ=255)
                                                                                                                        ggCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCaa                                                                                                                      >  6:354060/1‑58 (MQ=255)
                                                                                                                           tGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTc                         >  3:88878/1‑148 (MQ=255)
                                                                                                                             aTGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCaa                       <  2:107605/148‑1 (MQ=255)
                                                                                                                                 ggCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGATCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAgg                  <  6:76700/149‑1 (MQ=255)
                                                                                                                                    cATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCa               >  1:325066/1‑149 (MQ=255)
                                                                                                                                     aTGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCAc              >  3:9365/1‑149 (MQ=255)
                                                                                                                                       gCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCAcgc            >  8:17789/1‑149 (MQ=255)
                                                                                                                                                tGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTACCGtt    >  3:111797/1‑148 (MQ=255)
                                                                                                                                                 ggCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCGTTaa  >  5:74782/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAAC‑‑‑TCGACCACATCGACA‑‑GTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCGTTAA  >  NC_002947/4360148‑4360436

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: