Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,621,878 G→A 7.2% L264L (CTG→CTA) ‡ PP_1421 → sensor histidine kinase TctE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,621,8780GA 7.2% 199.6 / 4.3 69L264L (CTG→CTA) ‡PP_1421sensor histidine kinase TctE
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base A (2/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

GAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAGAATG  >  NC_002947/1621736‑1622023
                                                                                                                                              |                                                                                                                                                   
gAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGgcggcg                                                                                                                                               >  3:302536/1‑149 (MQ=255)
   cGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCCAGCGGCGttc                                                                                                                                            <  1:35519/149‑2 (MQ=255)
     ttGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                          <  5:346752/149‑1 (MQ=255)
        ccGGTGGTGCTGGTGCAGCGCGAGCTTGGCCCGCTGGAGCGCTCCTTGAACCACTATACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCTGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGc                                                                                                                                       <  4:204840/149‑1 (MQ=255)
         cGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCCAGCGGCGttcgcagtt                                                                                                                                      <  7:212943/149‑8 (MQ=255)
                gcaggtgcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCg                                                                                                                               >  7:261468/1‑149 (MQ=255)
                       cagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTggg                                                                                                                        >  6:85878/1‑149 (MQ=255)
                             gAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTgcg                                                                                                                  <  2:271686/149‑1 (MQ=255)
                                     cccGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGcc                                                                                                           >  1:375997/1‑148 (MQ=255)
                                     cccGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCg                                                                                                          >  7:41675/1‑149 (MQ=255)
                                     cccGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGAGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGcc                                                                                                           >  6:143988/1‑148 (MQ=255)
                                      ccGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGcc                                                                                                           >  4:265480/1‑147 (MQ=255)
                                       cGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAg                                                                                                        <  4:159011/149‑1 (MQ=255)
                                         cTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGcc                                                                                                      >  8:189347/1‑149 (MQ=255)
                                          tGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCg                                                                                                     >  3:108784/1‑149 (MQ=255)
                                            gTGCGCGCCTTGAA‑CACTTTCCCGAGCCCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAg                                                                                                  <  8:125943/149‑1 (MQ=255)
                                                gcgcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAg                                                                                              <  6:192921/149‑1 (MQ=255)
                                                   ccTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTgg                                                                                           >  6:126222/1‑149 (MQ=255)
                                                    cTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGACCTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTggc                                                                                          <  6:262588/149‑1 (MQ=255)
                                                    cTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTggc                                                                                          <  5:124419/149‑1 (MQ=255)
                                                     ttGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTggcg                                                                                         >  7:250861/1‑149 (MQ=255)
                                                     ttGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAAGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTggcg                                                                                         >  1:235911/1‑149 (MQ=255)
                                                       gAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTggcgg                                                                                        <  5:357777/148‑1 (MQ=255)
                                                        aaCCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTggcggc                                                                                       <  7:189347/148‑1 (MQ=255)
                                                        aaCCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTggcgg                                                                                        >  1:130322/1‑147 (MQ=255)
                                                        aaCCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCa                                                                                      >  4:275488/1‑149 (MQ=255)
                                                          ccACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAg                                                                                                  <  3:245031/136‑1 (MQ=255)
                                                          ccACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAg                                                                                                  >  4:245031/1‑136 (MQ=255)
                                                           cACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGAc                                                                                   >  4:59012/1‑149 (MQ=255)
                                                           cACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGAc                                                                                   >  2:330324/1‑149 (MQ=255)
                                                                tACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTg                                                                               <  3:349826/148‑1 (MQ=255)
                                                                 aCCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTg                                                                               >  5:77716/1‑147 (MQ=255)
                                                                  ccGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGa                                                                             <  2:130322/148‑1 (MQ=255)
                                                                   cGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGt                                                                           <  5:85878/149‑1 (MQ=255)
                                                                    gAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTc                                                                          <  8:250861/149‑1 (MQ=255)
                                                                    gAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAAGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTc                                                                          <  2:235911/149‑1 (MQ=255)
                                                                     aGCGCTTGCGCGGCCATTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCg                                                                         <  2:302208/149‑1 (MQ=255)
                                                                      gcgcTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATCGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCgg                                                                        >  3:56908/1‑149 (MQ=255)
                                                                         cTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAg                                                                     >  3:30742/1‑149 (MQ=255)
                                                                            gcgcGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGcc                                                                   <  4:302536/148‑1 (MQ=255)
                                                                              gcgGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAg                                                                                                        <  1:196052/110‑1 (MQ=255)
                                                                              gcgGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAg                                                                                                        >  2:196052/1‑110 (MQ=255)
                                                                              gcgGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAg                                                                >  6:228176/1‑149 (MQ=255)
                                                                                          aGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAgg                                                    <  5:143988/149‑1 (MQ=255)
                                                                                           gcgcCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGAAAGGc                                                   >  4:247877/1‑149 (MQ=255)
                                                                                                aGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGAcc                                              <  4:30742/149‑1 (MQ=255)
                                                                                                 ggCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGAccc                                             <  2:375997/149‑1 (MQ=255)
                                                                                                   cgcAGGTCATTGCCGCCGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCCCCGCCAGGCTGACCCAc                                           >  5:359371/1‑149 (MQ=255)
                                                                                                     cAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCt                                         <  5:228176/149‑1 (MQ=255)
                                                                                                      aGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTg                                        >  7:70830/1‑149 (MQ=255)
                                                                                                        ttCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑Gc                                     <  3:253557/149‑1 (MQ=255)
                                                                                                               ccGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGACGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAg                              >  2:257167/1‑149 (MQ=255)
                                                                                                                  aCGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCtg                           <  1:330324/149‑1 (MQ=255)
                                                                                                                    gcgGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCtgtt                         <  4:108784/149‑1 (MQ=255)
                                                                                                                      ggCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGATAGGCTGACCCACCTG‑GCCAACCAGCtgttgt                       <  3:247877/149‑1 (MQ=255)
                                                                                                                         ccACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGCCAGGCTGACACACCTG‑GCCAACCAGCTGTTGTccc                    >  4:220576/1‑149 (MQ=255)
                                                                                                                           aCGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTg                  >  5:73578/1‑149 (MQ=255)
                                                                                                                                ccttcaaCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTGGcgcg             >  3:126981/6‑149 (MQ=255)
                                                                                                                                cTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTGGcgcg             <  3:59012/149‑1 (MQ=255)
                                                                                                                                   tcaaCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTAGCGCGGGt          >  4:55929/3‑149 (MQ=255)
                                                                                                                                   cGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGTGAGTGGCGGCAGACATTGGAGTCGGCAGCCCAGGGCACCCACAGCATGACCAACATG‑GCAAACAAGCTGTTGTCCATGGCGCggg           >  1:113989/1‑149 (MQ=255)
                                                                                                                                   cGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGAGAGTGGCGGCAGACATTGGAGTCGGCAGCCCAGGGCCCCGACAGGATGACCAACATG‑GCCAACAAGATGTTGTCCATGGCGCggg           >  1:109185/1‑149 (MQ=255)
                                                                                                                                     aaCGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGCGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTGGCGCGGGTCg        >  8:295242/1‑149 (MQ=255)
                                                                                                                                     aaCGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTGGCGCGGGTCg        <  2:81512/149‑1 (MQ=255)
                                                                                                                                      aCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTGGCGCGGGTCga       <  4:126981/149‑1 (MQ=255)
                                                                                                                                          cgcTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAGAAt   >  8:316517/1‑149 (MQ=255)
                                                                                                                                           gcTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAGAATg  >  7:174765/1‑149 (MQ=255)
                                                                                                                                           gcTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAGAATg  <  6:77716/149‑1 (MQ=255)
                                                                                                                                            cTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGCGCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAGAATg  <  1:257167/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAG‑GAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTG‑GCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAGAATG  >  NC_002947/1621736‑1622023

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