Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 94,007 G→A 5.5% G102S (GGC→AGC)  trpI → HTH‑type transcriptional regulator TrpI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294794,0070GA 5.5% 274.6 / 2.6 89G102S (GGC→AGC) trpIHTH‑type transcriptional regulator TrpI
Reads supporting (aligned to +/- strand):  ref base G (39/45);  new base A (3/2);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAAGGACGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGCTGGTGTATGCCGAGC  >  NC_002947/93863‑94151
                                                                                                                                                |                                                                                                                                                  
gTCAAGGACGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAg                                                                                                                                               >  1:173875/1‑149 (MQ=255)
 tCAAGGACGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGc                                                                                                                                              >  5:27040/1‑149 (MQ=255)
  cAAGGACGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCCCGAGCAGcc                                                                                                                                             >  8:204484/1‑149 (MQ=255)
   aaGGACGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGc                                                                                                                                              >  5:103271/1‑147 (MQ=255)
   aaGGACGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGc                                                                                                                                              <  6:103271/147‑1 (MQ=255)
   aaGGACGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCct                                                                                                                                            <  6:212685/149‑1 (MQ=255)
   aaGGACGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCct                                                                                                                                            >  1:319667/1‑149 (MQ=255)
       aCGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCCCGAGCAGcccagca                                                                                                                                        <  7:204484/149‑6 (MQ=255)
       aCGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGGGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCctgctg                                                                                                                                        <  8:233316/149‑1 (MQ=255)
             gTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCgg                                                                                                                                 <  5:108265/149‑1 (MQ=255)
                  aTCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGtt                                                                                                                            >  3:26786/1‑149 (MQ=255)
                   tCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGttt                                                                                                                           <  5:194829/149‑1 (MQ=255)
                   tCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGtt                                                                                                                            >  6:194829/1‑148 (MQ=255)
                    cAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTa                                                                                                                          >  7:175313/1‑149 (MQ=255)
                     aaaCTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTAt                                                                                                                         >  4:212388/1‑149 (MQ=255)
                         tCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCgg                                                                                                                                 <  4:202027/137‑1 (MQ=255)
                         tCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCgg                                                                                                                                 >  3:202027/1‑137 (MQ=255)
                         tCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCcg                                                                                                                     >  5:166007/1‑149 (MQ=255)
                         tCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCcg                                                                                                                     >  8:167015/1‑149 (MQ=255)
                           acaGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCcgcg                                                                                                                   <  6:27040/149‑1 (MQ=255)
                               aTGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCggt                                                                                                                                >  4:288275/1‑132 (MQ=255)
                               aTGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCggt                                                                                                                                <  3:288275/132‑1 (MQ=255)
                                     gTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCgg                                                                                                         >  2:192698/1‑149 (MQ=255)
                                      tATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTggg                                                                                                             >  5:50196/1‑144 (MQ=255)
                                      tATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTggg                                                                                                             <  6:50196/144‑1 (MQ=255)
                                        tGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGtt                                                                                                       <  2:101129/148‑1 (MQ=255)
                                          cGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGc                                                                                                     <  7:147321/148‑1 (MQ=255)
                                           gTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAg                                                                                                   >  5:86020/1‑149 (MQ=255)
                                                cgcgATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGTCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCgg                                                                                               <  1:86201/149‑1 (MQ=255)
                                                cgcgATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGc                                                                                                 <  5:30253/146‑1 (MQ=255)
                                                        cAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCt                                                                                       >  7:163522/1‑148 (MQ=255)
                                                        cAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCt                                                                                       <  8:163522/148‑1 (MQ=255)
                                                        cAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCt                                                                                       >  4:335240/1‑148 (MQ=255)
                                                        cAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTg                                                                                      >  7:137770/1‑149 (MQ=255)
                                                                aGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTgcgc                                                                              <  3:335240/149‑1 (MQ=255)
                                                                    cTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGa                                                                                     <  6:328742/138‑1 (MQ=255)
                                                                    cTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGa                                                                                     >  5:328742/1‑138 (MQ=255)
                                                                    cTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTg                                                                           >  2:153673/1‑148 (MQ=255)
                                                                       cGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCAc                                                                        <  2:319667/148‑1 (MQ=255)
                                                                       cGACCGGCTGCGCAGCGTGTGTTCCGAGCGCAGCCGTGATGGGAGCGATGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGTCTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGc                                                                          <  1:94556/147‑1 (MQ=255)
                                                                          ccGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGt                                                                    <  3:212388/149‑1 (MQ=255)
                                                                          ccGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGGTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGt                                                                    >  3:53529/1‑149 (MQ=255)
                                                                            ggCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGt                                                                    >  6:49390/1‑147 (MQ=255)
                                                                            ggCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGt                                                                    <  5:49390/147‑1 (MQ=255)
                                                                             gCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCt                                                                                       <  4:290501/127‑1 (MQ=255)
                                                                             gCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCt                                                                                       >  3:290501/1‑127 (MQ=255)
                                                                              cTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTcggc                                                                <  6:261046/149‑1 (MQ=255)
                                                                                 cgcAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTcggcgg                                                              <  4:26786/148‑1 (MQ=255)
                                                                                 cgcAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTcggcgg                                                              <  2:173875/148‑1 (MQ=255)
                                                                                  gcAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGt                                                            >  7:44978/1‑149 (MQ=255)
                                                                                   cAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTCCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTg                                                           <  1:134470/149‑1 (MQ=255)
                                                                                     gCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCtt                                                                            <  2:312599/130‑1 (MQ=255)
                                                                                     gCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAg                                                         >  5:274008/1‑149 (MQ=255)
                                                                                     gCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAg                                                         >  4:28640/1‑149 (MQ=255)
                                                                                      cGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGGCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAgg                                                        <  8:137770/149‑1 (MQ=255)
                                                                                      cGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAgg                                                        <  6:166007/149‑1 (MQ=255)
                                                                                         gtgtTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTgcgc                                                                              >  6:241456/1‑124 (MQ=255)
                                                                                         gtgtTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTgcgc                                                                              <  5:241456/124‑1 (MQ=255)
                                                                                          tgtTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGc                                                                                         >  5:53434/1‑112 (MQ=255)
                                                                                          tgtTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGc                                                                                         <  6:53434/112‑1 (MQ=255)
                                                                                            ttCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGAc                                                   <  6:86020/148‑1 (MQ=255)
                                                                                            ttCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGAc                                                   <  8:44978/148‑1 (MQ=255)
                                                                                                gAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGa                                              <  7:167015/149‑1 (MQ=255)
                                                                                                gAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGa                                              <  4:53529/149‑1 (MQ=255)
                                                                                                gAGCTCAGCCGTGATGCGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGa                                              >  5:52898/1‑149 (MQ=255)
                                                                                                    tCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCcg                                          >  5:72852/1‑149 (MQ=255)
                                                                                                          gTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGc                                                                                                 <  8:137322/88‑1 (MQ=255)
                                                                                                          gTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGa                                                          <  2:244864/127‑1 (MQ=255)
                                                                                                          gTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCAGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGa                                                          >  1:244864/1‑127 (MQ=255)
                                                                                                              tgtgAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCcgcg                                        >  2:332567/1‑141 (MQ=255)
                                                                                                              tgtgAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCcgcg                                        <  1:332567/141‑1 (MQ=255)
                                                                                                                 gAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGc                             >  3:100509/1‑149 (MQ=255)
                                                                                                                  aGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCt                            >  1:264815/1‑149 (MQ=255)
                                                                                                                     gAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGa                         <  8:175313/149‑1 (MQ=255)
                                                                                                                      aGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGAt                        <  6:56294/149‑1 (MQ=255)
                                                                                                                          gcCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTgcgc                                                                              <  3:316397/91‑1 (MQ=255)
                                                                                                                          gcCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTgcgc                                                                              >  4:316397/1‑91 (MQ=255)
                                                                                                                          gcCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGt                    <  3:28640/149‑1 (MQ=255)
                                                                                                                             aTTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGc                 >  3:279315/1‑149 (MQ=255)
                                                                                                                              ttCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGCt                <  4:279315/149‑1 (MQ=255)
                                                                                                                                  cagcaggCTGCCCGAGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGCTGgtg            >  6:96333/6‑149 (MQ=255)
                                                                                                                                       ggCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGCTGGTGTATGc       >  2:66082/1‑149 (MQ=255)
                                                                                                                                       ggCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGCTGGTGTATGc       <  1:192698/149‑1 (MQ=255)
                                                                                                                                        gctgccCGAGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTg                           <  1:167424/128‑1 (MQ=255)
                                                                                                                                        gctgccCGAGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTg                           >  2:167424/1‑128 (MQ=255)
                                                                                                                                            ctCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGCTGGTGTATGCCgcgc  >  8:71136/1‑149 (MQ=255)
                                                                                                                                            ctCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGCTGGTGTATGCCGa    <  6:72852/147‑1 (MQ=255)
                                                                                                                                            ctCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGCTGGTGTATGCCGa    <  2:264815/147‑1 (MQ=255)
                                                                                                                                            ctCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGCTGGTGTATGCCGAGc  >  7:293908/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GTCAAGGACGGGCGTGGCATCAAACTCACAGATGCCGGTATGCGTCTGCGCGATGCCAGCGGTGAGGCCTTCGACCGGCTGCGCAGCGTGTGTTCCGAGCTCAGCCGTGATGTGAGCGAGGCGCCATTCGTGCTGGGCTGCTCGGGCAGCCTGCTGGCGCGGTGGTTTATTCCGCGGCTGGGGCGGTTGCAGGCGGATTTGCCTGAGCTGCGCTTGCACCTGTCGGCGGGTGAGGGCGACCTTGACCCGCGGCGGCCCGGGCTGGATGCGTTGCTGGTGTATGCCGAGC  >  NC_002947/93863‑94151

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: