Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,816,594 T→A 7.0% I21N (ATT→AAT)  PP_2472 → MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,816,5940TA 7.0% 211.5 / 3.9 71I21N (ATT→AAT) PP_2472MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (28/38);  new base A (2/3);  total (30/41)
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.31e-01

CAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGTG  >  NC_002947/2816447‑2816740
                                                                                                                                                     |                                                                                                                                                  
cAAGACAGCGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACGGGGCAAACGCTACTTCACCAtt                                                                                                                                                   <  7:60018/149‑1 (MQ=255)
 aaGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTg                                                                                                                                                  <  3:30105/149‑1 (MQ=255)
    aCAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTg                                                                                                                                               <  4:119590/149‑1 (MQ=255)
        gggAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGt                                                                                                                                           <  5:51554/149‑1 (MQ=255)
        gggAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGt                                                                                                                                           <  3:15633/149‑1 (MQ=255)
         ggAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCAATGGTGAAGta                                                                                                                                          <  6:149248/149‑2 (MQ=255)
         ggAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGGTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCAATGGTGAAGta                                                                                                                                          <  2:257691/149‑2 (MQ=255)
            agagATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGc                                                                                                                                       <  5:97514/149‑1 (MQ=255)
               gaTCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAg                                                                                                                                    <  5:5252/149‑1 (MQ=255)
                   ccGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTtgt                                                                                                                                <  1:187574/149‑1 (MQ=255)
                    cGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTtgtg                                                                                                                               >  1:110413/1‑149 (MQ=255)
                     gATCACCGCACGGCGCGTCGTCACCTTTCGTCCATGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATTCCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTtgtgt                                                                                                                              <  3:102760/149‑1 (MQ=255)
                      aTCACCGCACGGCGCGTCGACACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTtgtgt                                                                                                                              <  5:341730/149‑1 (MQ=255)
                           cGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGt                                                                                                                         >  6:289290/1‑148 (MQ=255)
                             cACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTa                                                                                                                        >  2:296817/1‑147 (MQ=255)
                                ggCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTT‑‑CCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCcac                                                                                                                 <  5:44487/149‑1 (MQ=255)
                                 gcgcGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCca                                                                                                                  <  7:138034/149‑1 (MQ=255)
                                         tcACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGtgct                                                                                                          >  2:75002/1‑149 (MQ=255)
                                           aCCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAACCCATAACGACGAACTGCCCCCCAGACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGtgctg                                                                                                         <  6:282717/148‑1 (MQ=255)
                                             cTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCgg                                                                                                      >  5:260749/1‑149 (MQ=255)
                                                 cGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTAt                                                                                                   <  5:256387/148‑1 (MQ=255)
                                                  gTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGc                                                                                                                                       <  7:66652/111‑1 (MQ=255)
                                                  gTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGc                                                                                                                                       >  8:66652/1‑111 (MQ=255)
                                                  gTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAACGCCACACGTGCTGCGGAATTg                                                                                                 >  8:285145/1‑149 (MQ=255)
                                                   tCCAGGGCATAAGTTGAAGTCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTaa                                                                                                                       <  2:50264/126‑1 (MQ=255)
                                                   tCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTgg                                                                                                >  5:135552/1‑149 (MQ=255)
                                                   tCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAAACAAGCCATAACGCCGAAATGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTaa                                                                                                                       >  1:50264/1‑126 (MQ=255)
                                                     cAGGGAAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCAAAACGACGAAAAGCCCCCCATACCCGGCAAACGCTACTTAACCATTGGAGAAGTGAGCGAGTTCTGTGCCGTAAAGCCACACGTGAAGCGGTATTggg                                                                                               >  8:141447/1‑148 (MQ=255)
                                                      aGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAg                                                                                             <  4:111519/149‑1 (MQ=255)
                                                       gggCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGc                                                                                            <  6:260749/149‑1 (MQ=255)
                                                        ggCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAg                                                                                             >  4:366534/1‑147 (MQ=255)
                                                        ggCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCa                                                                                           <  2:168517/149‑1 (MQ=255)
                                                         gCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGAATTGGGAGCAg                                                                                          >  1:333996/1‑149 (MQ=255)
                                                           aGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAg                                                                                          >  2:152658/1‑147 (MQ=255)
                                                                  gAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGAc                                                                                 >  6:335423/1‑149 (MQ=255)
                                                                   aaGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGAc                                                                                 <  4:365888/148‑1 (MQ=255)
                                                                     ggCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAg                                                                              >  8:87341/1‑149 (MQ=255)
                                                                      gCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTg                                                                                                 <  1:210344/129‑1 (MQ=255)
                                                                      gCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTg                                                                                                 >  2:210344/1‑129 (MQ=255)
                                                                      gCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGc                                                                             <  1:296817/149‑1 (MQ=255)
                                                                               gatgCCTATGCTGGAAGCAATCCATAACGACGAACAGCCCCCCATACCGGGCAAACGCTACTTCACCTTTGGTGAAGTTAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGTGTTCGACCAGCTGAACCCt                                                                     <  7:141447/145‑1 (MQ=255)
                                                                                aGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCGCCCCATTCCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCTCACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCtgt                                                                   <  7:157152/149‑1 (MQ=255)
                                                                                aGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCtgt                                                                   >  4:364785/1‑149 (MQ=255)
                                                                                 ggCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCtgt                                                                   <  2:110413/148‑1 (MQ=255)
                                                                                  gCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCtgtg                                                                  >  4:337572/1‑148 (MQ=255)
                                                                                       tGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTTTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAgc                                                              <  5:289290/147‑1 (MQ=255)
                                                                                       tGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAgcgg                                                            <  6:345682/149‑1 (MQ=255)
                                                                                        gCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAgcggc                                                           >  4:282050/1‑149 (MQ=255)
                                                                                            gAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGcggcgg                                                         <  1:189954/147‑1 (MQ=255)
                                                                                            gAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCgggg                                                       >  3:307652/1‑149 (MQ=255)
                                                                                              aCCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACAGATTGTTAAATTGAGCAATTTTTTTCCCTAAAACCCACACTTGATGCGGTATTGGGACCGGGATTGCGACAAGCTCAACCATGTGAACGGGCgggg                                                       >  1:300150/1‑148 (MQ=255)
                                                                                              aCCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCt                                                                     >  6:187460/1‑133 (MQ=255)
                                                                                              aCCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCt                                                                     <  5:187460/133‑1 (MQ=255)
                                                                                                cAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAAc                                                   >  1:311043/1‑149 (MQ=255)
                                                                                                    ccATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACcggc                                               >  3:84511/1‑149 (MQ=255)
                                                                                                      aTAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGTGGCGGGGGAACcggcgg                                             <  1:75002/149‑1 (MQ=255)
                                                                                                      aTAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACcggcgg                                             <  7:87341/149‑1 (MQ=255)
                                                                                                              gaACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCa                                     <  3:337572/149‑1 (MQ=255)
                                                                                                                aCTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAgc                                   <  2:300150/149‑1 (MQ=255)
                                                                                                                           aTACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACcggcg                                              >  6:313895/1‑128 (MQ=255)
                                                                                                                           aTACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACcggcg                                              <  5:313895/128‑1 (MQ=255)
                                                                                                                                gggCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtga                    <  4:84511/149‑1 (MQ=255)
                                                                                                                                gggCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtga                    <  4:307652/149‑1 (MQ=255)
                                                                                                                                  gCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTATGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtgat                   >  1:267585/1‑148 (MQ=255)
                                                                                                                                      aCGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATcc              <  2:311043/149‑1 (MQ=255)
                                                                                                                                         tcaCTTCACCAATGGTGAAGTAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGc            >  2:221546/3‑147 (MQ=255)
                                                                                                                                            cTTCACCAATGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTGTTGGGAGCAGGa                                                                                        >  6:149098/1‑69 (MQ=255)
                                                                                                                                            cTTCACCAATGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTGTTGGGAGCAGGa                                                                                        <  5:149098/69‑1 (MQ=255)
                                                                                                                                             ttCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGAt       >  4:278615/1‑149 (MQ=255)
                                                                                                                                                aCCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCg    >  8:194393/1‑149 (MQ=255)
                                                                                                                                                 ccATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGt   <  3:282050/149‑1 (MQ=255)
                                                                                                                                                   cATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGTg  <  3:123652/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGTG  >  NC_002947/2816447‑2816740

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: