Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT39.2% 56.3 / 34.8 51S62F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (5/26);  new base T (19/1);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 2.64e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCA  >  NC_000913/4604205‑4604483
                                                                                                                                           |                                                                                                                                            
aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTtttttc                                                                                                                                              >  1:75451/1‑140 (MQ=255)
    gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                                                                          >  1:481211/1‑140 (MQ=255)
       tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                       <  1:142760/140‑1 (MQ=255)
         aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                     <  1:684803/140‑1 (MQ=255)
         aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                     <  2:730127/140‑1 (MQ=255)
         aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                     <  1:611932/140‑1 (MQ=255)
                 taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGaa                                                                                                                             >  1:8287/1‑140 (MQ=255)
                   tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGt                                                                                                                           <  2:267899/140‑1 (MQ=255)
                   tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                           <  2:195350/140‑1 (MQ=255)
                        aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                                                                      <  2:98874/140‑1 (MQ=255)
                        aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                                                                      <  2:540941/140‑1 (MQ=255)
                        aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                                                                      <  1:31553/140‑1 (MQ=255)
                          tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa                                                                                                                    <  2:75451/140‑1 (MQ=255)
                              agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATCAGAAATGGGCGCagagg                                                                                                                >  2:31553/1‑139 (MQ=255)
                               gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagagag                                                                                                              <  2:596846/140‑1 (MQ=255)
                                  tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                                                                           <  2:481211/140‑1 (MQ=255)
                                                         ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGc                                                                                    <  2:8287/140‑1 (MQ=255)
                                                              taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTa                                                                               >  1:711246/1‑140 (MQ=255)
                                                                   tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGATCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt                                                                          >  2:600555/1‑140 (MQ=255)
                                                                       tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGTGAGGGCCGAGAGGGGTTATTTTGCCTGACTGAACTGGCGATTAAGTTTCcgc                                                                      >  2:405981/1‑140 (MQ=255)
                                                                           tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAc                                                                  <  1:187864/140‑1 (MQ=255)
                                                                            ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                                 <  1:141519/140‑1 (MQ=255)
                                                                             taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCTTTTTTTTTTTTTTTGGCATGTGAAGTGAGAGGAGTGAGGGGTTGTTTTGGCTGGGTTAACTAGCGACTAAGTTTCtgcctgtcc                                                                >  2:443452/1‑134 (MQ=255)
                                                                              cagaGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg                                                               <  1:75519/137‑1 (MQ=255)
                                                                               ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                                                              <  1:112016/140‑1 (MQ=255)
                                                                                  gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTTAGAAAAGAGCGCCGAGAGGGGTTATTTTTTCCGACTGAACCGGGGGTTTAAATTCCGCGGCCCCcgcgg                                                           >  2:562953/1‑139 (MQ=255)
                                                                                  gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAATTTCCGCGTACCCcgcgg                                                           >  2:305821/1‑139 (MQ=255)
                                                                                  gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt                                                           >  2:280075/1‑140 (MQ=255)
                                                                                   tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCACTTAGAAGGGAACGCAGAGCGGGATTTTCTTTGCTGACTGAAATTTGCATTAAATTTCCGTCTGCCCGGCgct                                                          >  2:112500/1‑138 (MQ=255)
                                                                                    cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCCGCGgtt                                                         >  2:154437/1‑140 (MQ=255)
                                                                                     aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTTCATGAGGAGGGAGCGCAGAGGGGGGTTTTTTTGCCTGACTGAAATGGCGGTTTAATTTTCGCCGACCCCGCGgttt                                                        >  2:187864/1‑139 (MQ=255)
                                                                                     aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGGGGGGGTATCTTTCCTGACTGAAATGGCGGTTAAGTTTCCCCGGACCCGGCGgttt                                                        >  2:397696/1‑139 (MQ=255)
                                                                                     aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCCGAGAGGGGTTATCTTTCCTGACTGGACTGGCGATTAAGTTTTCGCGGACCCGGCGgttt                                                        >  2:711853/1‑139 (MQ=255)
                                                                                      ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTTCATGTGGAATGAGCGGTGAGAGGGGATAATCTGGCTGACTGAACATGCGTATAAATTTCCCCCAACCCGGCGGTTAg                                                       >  2:134378/1‑140 (MQ=255)
                                                                                          tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCCGAGGGGGGTTATCTTTCCTGGCTGGAATGGCGGTTAAAATTTCTCGTACCCGGCGCTTTGtttt                                                   >  2:169978/1‑138 (MQ=255)
                                                                                            aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCTTTGGAAGTGAGCGCGGAGGGGGGTTTTCTTGCCTGCCTGAACTGGCGGTTAAGTTTTCCCGCAACCCGCGTTTAGGTAttt                                                 >  2:141711/1‑139 (MQ=255)
                                                                                            aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGGCGAGAGGGGTTATCTTTCCTGACTGAACTGGGGATTAAGTTTTCGCGTACCCCGCGTTTTGTTaatt                                                 >  2:75519/1‑137 (MQ=255)
                                                                                            aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                 <  1:711853/140‑1 (MQ=255)
                                                                                            aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTCTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCCCGTACCCGGCGTTTAGTTATTg                                                 <  2:371291/140‑1 (MQ=255)
                                                                                             tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGGGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGGACCCGGCGTTTTGTTATTgt                                                >  2:239217/1‑139 (MQ=255)
                                                                                             tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCCGAGAGGGGGTATCTTGCCTGACTGAACTGGCGGTTTAGTTTCCCCGTACCCCGCGGTTTGTTATTgg                                                >  2:72790/1‑139 (MQ=255)
                                                                                              ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                               <  1:305821/140‑1 (MQ=255)
                                                                                              ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                               <  1:405981/140‑1 (MQ=255)
                                                                                                cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                             <  2:711246/140‑1 (MQ=255)
                                                                                                cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                             <  2:706310/140‑1 (MQ=255)
                                                                                                   gATGCCATGTCGCGGATCTCGTTTGCGGCGATCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga                                          <  1:600555/140‑1 (MQ=255)
                                                                                                    aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat                                         <  2:219054/140‑1 (MQ=255)
                                                                                                    aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat                                         <  1:169978/140‑1 (MQ=255)
                                                                                                                 gATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGGGAGGGGGTATTTTGCCTGACTGAAATGGGGGGTAAGTTTTCGCGGACCCGGCGTTTTGTTATTTTGGGTGATGATATTTAAGCTc                            >  2:623604/1‑140 (MQ=255)
                                                                                                                   tctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg                          <  1:280075/140‑1 (MQ=255)
                                                                                                                     tcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt                        <  1:72790/140‑1 (MQ=255)
                                                                                                                            cggcggTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGGGGGGGTTTTTTTTCTGACTGAACTGGCGGTTTAGTTTTCGCGTACCCCGCGGTTTGTTATTTTGGATGATGATATTGAAGATCGGCTGATTTGt                 >  2:397559/1‑140 (MQ=255)
                                                                                                                            cggcggTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATATTTCCTGACTGAACTGGCGGATAAATTTCCGCGGACCCGGCGTTTAGTTATTTTGGATGGTGAAATTTAAGCTCCGCTTTTTGGt                 >  2:81860/1‑140 (MQ=255)
                                                                                                                               cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa              <  1:562953/140‑1 (MQ=255)
                                                                                                                                 gataatttttcttttttTTCCATGAGAAGTTAGCGCAGCGAGGGATTATCTTTTCTGACTGAACCGGCGATTAAGATTTCGCGCTCCCCGCGGTTTGTTATTGCGGATGATGATATTTAAGCTCGGCTGATTGGTTCAtt            >  2:527103/6‑140 (MQ=255)
                                                                                                                                 gTCATTTTTTTTTTTTTTTCCCTGAGGAATGAACGCAGAGAGGGATTATCTTGCCTGTCTTAACCGGCGAATAAGTTTTCGCGTACCCCGCGGTTTGTTATTGCGGGTGATGATATTGAAGGTCGGCTGATTGGTTCAtt            >  2:667637/1‑140 (MQ=255)
                                                                                                                                 gTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGGGGGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTTGTTATTGCGGATGATGATATTGAAGCTCGGGTGATTTGTTCAtt            >  2:448745/1‑140 (MQ=255)
                                                                                                                                 gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAtt            >  1:28192/1‑140 (MQ=255)
                                                                                                                                         ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt    <  1:255325/140‑1 (MQ=255)
                                                                                                                                         ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt    <  1:112500/140‑1 (MQ=255)
                                                                                                                                         ttctttttctGCCACGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt    <  2:561875/140‑1 (MQ=255)
                                                                                                                                           ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCa  <  1:213738/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
AAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCA  >  NC_000913/4604205‑4604483

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: